AD 


r.TR.77-025 


A  DYNAMIC  ANALYSIS  OF  THE 
SELF-PROPELLED,  8-INCH,  M110A1, 
HOWITZER  UTILIZING  FORMAC 


TI^OMAS  D.  STREETER 
ROBERT  H.  COBERLY 

PTIC  QUALITY  mSPECTED  3 

JUNE  1977 

TECHNICAL  REPORT 


LARGE  CALIBER  WEAPONS  SYSTEMS  LABORATORY 


Ur  AILMY  ARMAAAENT  RESEARCH 


AND  DEVELOPMENT  COMMAND 
ROCK  ISLAND,  ILLINOIS  61201 


19970930  114 


UNCFASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Data  BnteretQ 


Streeter  5292 


4.  TITLE  (and  Subtitta) 

A  Dynamic  Analysis  of  the  Self-Propelled  8"  MllOAl 
Howitzer  Utilizing  FORMAC 


9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Weapons  Division 

Large  Caliber  Weapons  System  Laboratory 
Dover,  NJ  07801  _ 


11.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

[weapons  Division 

[Large  Caliber  Weapons  Systems  Laboratory 
[Dover,  NJ  07801 


REPORT  DOCUMENTATION  PAGE 


.  report  number 

R-TR-77-025 


2.  GOVT  ACCESSION  NO 


.  AUTHOR^*; 


THOMAS  D.  STREETER 
Robert  H.  COBERLY 


*4.  MONITORING  AGENCY  NAME  tk  AOORESSfl/  dtHarant  from  ControWni  OWce) 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


5.  TYPE  OF  REPORT  &  PERIOD  COVERED 

FINAL  REPORT 
Sep  76  -  Apr  77 


6.  PERFORMING  ORG.  REPORT  NUMBER 


8.  CONTRACT  OR  GRANT  NUMBERfsJ 

N/A 


10.  PROGRAM  ELEMENT,  PROJECT,  TASK 
AREA  A  WORK  UNIT  NUMBERS 

6. 26. 03. A 
1W662603AH78 


12.  REPORT  DATE 

June  1977 


13.  NUMBER  OF  PAGES 

150 


15.  SECURITY  CLASS,  fo/ fWa  reporO 

UNCLASSIFIED  ! 


15a.  DECLASSI  FI  CATION/ DOWN  GRADING 
SCHEDULE 


16.  DISTRIBUTION  STATEMENT  (of  ihia  Raport) 

Distribution  Unlimited 


17.  DISTRIBUTION  STATEMENT  (of  the  abatract  antarad  In  Bfock  30,  if  different  from  Report) 

Unlimited 


10.  SUPPLEMENTARY  NOTES 


19.  KEY  WORDS  (Continue  on  reverse  aide  if  neceaaary  and  identify  by  bfock  number) 


MllOAl 
FOKMAC 
DYNAMIC 
MATH  MODEL 


SIMULATION 

FORTRAN 

SELF-PROPELLED 

LAGRANGE 


20  ABSTRACT  (Continue  on  reverse  aide  it  neceaaary  and  identify  by  block  number) 

This  report  presents  a  new  method  for  solving  Lagrange's  equations  of 
motion  utilizing  FORMAC.  An  application  using  this  technique  is  given  with  an 
eleven  degree-of-freedom  problem  which  describes  the  motion  of  the  MllOAl, 
self-propelled  8"  howitzer  under  dynamic  conditions  of  firing. A  Computer  progranj 
has  been  written,  is  operational,  and  the  listing  is  contained  in  the  appendix. 
This  report  is  an  endeavor  to  automate  the  generation  and  solution  of  the 
equations  of  motion  for  dynamical  systems. 


DD 


FORM  1479 
AN  73  *473 


EDITION  OF  1  NOV  65  IS  OBSOLETE 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Data  Entered) 


CONTENTS 


Page 

ABSTRACT .  iv 

FORWARD .  iv 

1.0  INTRODUCTION .  1 

1.1  GENERAL  DISCUSSION  OF  MODELING .  1 

1.2  DESCRIPTION  OF  THE  MODEL . ^ 

2.0  SOLUTION  TECHNIQUE .  9 

2 . 1  OVERVIEW .  9 

2 . 2  LAGRANGES  EQUATION . 9 

2.3  FORMAC  APPROACH .  11 

2.4  KINETIC  ENERGY .  12 

2.5  OTHER  ENERGY  TERMS .  16 

2.6  NUMERICAL  INTEGRATION .  18 

3.0  APPLICATION  OF  SOLUTION  TECHNIQUE .  19 

3.1  SYMBOL  TABLE . 19 

3.2  TECHNIQUE  OF  GENERATING  ENERGY  EXPRESSIONS  AND .  29 

GENERALIZED  FORCES 

3.3  FORMAC  PROCEDURE .  35 

4.0  CONCLUSIONS  AND  RECOMMENDATIONS .  41 

REFERENCES . 

DISTRIBUTION .  ^5 

APPENDIX  A  -  DERIVATION  OF  ENERGY  EXPRESSIONS  AND  .  A-1 

GENERALIZED  FORCES 

APPENDIX  B  -  MATHEMATICS  OF  SOLUTION  TECHNIQUE .  B-1 

APPENDIX  C  -  FORMAC  PROGRAM .  C-1 

APPENDIX  D  -  THE  SQUARE  ROOT  METHOD .  D-1 

APPENDIX  E  -  REMOVAL  OF  IMBEDDED  TERMS .  E-1 

APPENDIX  F  -  FORTRAN  PROGRAM .  F-1 


iii 


ABSTRACT 


This  report  presents  a  new  method  for  solving  Lagrange’s  equations 
of  motion  utilizing  FORMAC*  An  application  using  this  technique  is 
given  with  an  eleven  degree-of^f reedom  problem  which  describes  the 
motion  of  the  MllOAl,  a  self-propelled  8”  howitzer  under  dynamic  con¬ 
ditions  of  firing.  A  computer  program  has  been  OTitten,  is  operational, 
and  the  listing  is  contained  in  the  appendix.  This  report  is  an  en¬ 
deavor  to  automate  the  generation  and  solution  of  the  equations  of 
motion  for  dynamical  systems. 


FORWARD 

In  October  1976,  it  was  requested  that  a  mathematical  model  be 
developed  of  the  MllOAl  to  study  the  dynamic  motion  of  the  system  dur¬ 
ing  firing  so  as  to  define  dynamic  loading  at  various  points  in  the 
weapon  system  for  a  subsequent  stress  analysis  of  weapon  components. 
This  project  was  funded  by  AMC  under  AMCMS  Code  664602,12.38900,  D,  A, 
Project  No.  1W664602D389 , 

The  authors  are  indebted  to  Miss  Ann  Marie  Lowder  and  Mrs,  Jan 
Keller  for  their  many  hours  of  typing  in  preparing  this  report. 

Also,  the  authors  appreciate  the  assistance  given  by  Mr.  A1  Smith  in 
helping  to  obtain  the  data  for  this  model. 
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1,0  INTRODUCTION 


1.1  General  Discussion  of  Modeling 

A  convenient  way  to  describe  a  physical  system  and  study  its  be¬ 
havior  is  to  represent  the  system  by  mathematical  equations.  Gener¬ 
ally,  the  resulting  equations  are  a  set  of  differential  equations 
based  on  the  laws  of  mechanics  and  the  geometry  of  the  physical  sys¬ 
tem.  The  use  of  a  computer  is  usually  necessary  to  solve  this  set  of 
equations  and  especially  if  the  system  under  study  is  complex  and  a 
Sufficient  amount  of  detail  is  required  to  characterize  it.  If  the 
use  of  a  computer  becomes  necessary,  the  mathematical  expressions, 
physical  constraints  and  logic  must  be  transformed  into  a  computer 
program.  The  equations  and  their  solution  technique  comprise  a  math¬ 
ematical  or  computer  model.  If  such  a  model  is  to  be  of  value,  this 
model  must  be  amenable  to  analysis  and  its  output  must  represent  the 
behavior  of  the  physical  system  with  sufficient  accuracy  so  that  useful 
information  can  be  obtained  concerning  the  actual  system. 

Since  real  world  conditions  cannot  be  modeled  exactly,  inherent 
inaccuracies  will  exist  in  the  model  because  the  modeler  must  resort 
to  simplifications  and  abstractions  of  the  actual,  physical  system. 

The  important  point  is  that  an  acceptable  level  of  confidence  be 
established  so  that  inferences  drawn  from  the  model  output  are  correct 
and  the  model  generates  the  same  behavior  characteristics  as  the 
actual  system.  Confidence  in  the  model  output  is  gained  through 
model  validation. 

Model  validation  is  a  check  of  the  agreement  between  the  behavior 
of  the  model  and  that  of  the  actual  system.  The  correctness  of  the 
model  can  only  be  measured  relative  to  the  physical  system.  Models 
may  be  validated  according  to  various  criteria.  One  such  criteria 
is  Bayes*  formula  for  conditional  probabiliities.  With  this  procedure, 
specific  model  outputs  are  compared  with  field  test  data,  and  the  con- 
probability  (that  the  model  is  valid  given  the  field  test 
results)  is  determined. 
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However,  many  times  the  system  to  be  modeled  is  still  on  the 
drawing  board  and  field  test  results  do  not  exist.  If  this  is  the 
case,  the  design  and  simulation  should  complement  each  other.  That 
is,  a  math  model  may  be  developed  while  the  design  of  the  system  is 
progressing  from  the  drawing  board  stage  to  system  integration.  Here 
the  model  may  be  used  to  provide  insight  into  how  the  system  will 
perform  under  dynamic  conditions.  Thus  the  model  may  be  used  to 
(1)  predict  dynamic  loading  on  critical  parts,  (2)  define  sensitivity 
of  weapon  performance  to  design  parameters,  (3)  evaluate  suggested 
weapon  modifications  and  (4)  provide  the  necessary  foundation  for 
design  optimization  studies. 

As  test  data  becomes  available,  it  is  necessary  to  "validate" 
and  "tune"  the  model.  If  a  correlation  between  test  results  and 
predicted  motion  is  poor  and  significant  motion  occurs  which  was 
neglected  in  the  model,  a  complete  revision  of  the  model  may  be  re¬ 
quired.  However,  if  all  significant  motion  has  been  accounted  for 
in  the  model  but  the  magnitude  and  timing  of  the  predicted  motion  is 
in  error,  a  "tuning"  of  the  model  is  in  order.  For  example,  the 
model  to  be  generated  in  this  report  requires  an  effective  spring 
rate  for  the  ground-weapon  interaction.  The  value  chosen  for  this 
parameter  may  be  poorly  estimated.  If  varying  the  value  of  this 
parameter  between  reasonable  limits  greatly  Improves  the  desired 
correlation  between  model  and  test,  the  model  is  said  to  have  been 
tuned . 

Xn  the  design,  development,  and  fabrication  of  any  given  system, 
the  proper  mixture  of  analysis,  simulation,  and  testing  is  necessary 
to  produce  the  payoffs  that  must  be  achieved  for  a  cost-effective 
product.  The  use  of  model  studies  may  reduce  costs  in  lieu  of  exten¬ 
sive  testing,  but  it  does  not  necessarily  preclude  the  need  for  per¬ 
forming  an  actual  field  test.  A  field  test  yields  Information  for 
the  evaluation  of  the  system  itself,  provides  a  data  base  to  be  used 
for  comparison  to  other  systems,  and  is  necessary  for  model  validation. 
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One  of  the  most  important  and  useful  tools  to  analyze  the  results 
of  a  given  design  is  through  mathematical  simulation*  Such  is  the 
case  for  the  development  of  the  model  contained  in  this  report*  Sim— 
ulation  is  a  powerful  tool  and  it  enables  the  analyst  to  become  famil¬ 
iar  with  the  behavior  or  performance  of  the  actual  system  under  study 
when  subjected  to  a  variety  of  different  conditions  or  parameter 
changes.  Generally,  experiments  with  the  actual  system  itself  are 
very  costly;  however,  they  can  be  performed  on  the  model  with  relative 
ease  and  at  low  cost.  Many  times,  experimentation  on  the  model  will 
provide  more  information  about  the  interaction  of  variables  than  test¬ 
ing  on  the  actual  system  because  of  the  controlled  environment  and  the 
ease  of  parameter  variation. 

Weapon  systems  of  today  are  becoming  more  complex  and  this  trend 
is  likely  to  continue  because  of  the  current  threats  being  proposed; 
as  a  result,  design  requirements  are  becoming  increasingly  more 
difficult  to  satisfy.  It  is  well  recognized  that  a  change  in  design 
or  a  change  in  even  one  aspect  of  the  weapon  system  may  very  well 
produce  changes  or  create  the  need  for  changes  in  other  parts  of  the 
system.  As  the  long  list  of  parameter  trade-offs  such  as  accuracy, 
caliber,  threat,  dispersion,  etc.,  are  considered,  the  engineer's 
intuition  and  experience  become  increasingly  more  difficult  to  apply 
and  it  becomes  more  important  to  define  the  design  procedure  mathemat¬ 
ically.  Design  requirements  specify  that  a  weapon  system  is  to  per¬ 
form  some  task  at  some  index  of  performance.  Thus,  the  design  of  a 
weapon  system  provides  a  natural  setting  for  an  optimization  problem, 
assuming  a  knowledge  of  all  environmental  factors  influencing  the 
design  process  are  known.  To  just  search  for  admissible  parameters 
so  that  the  system  is  enabled  to  perform  its  task  is  not  satisfactory. 
It  should  be  required  to  seek  those  parameters  so  that  performance  is 
optimized  (in  some  sense) . 
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In  the  past,  conventional  methods  of  analyzing  the  dynamic 
behavior  of  a  given  system  required  the  analyst  to  linearize  the 
generalized  coordinates  to  achieve  simplification  of  expressions. 

This  is  no  longer  necessary.  FORMAC,  a  language  developed  by  IBM, 
allows  the  analyst  to  proceed  with  a  nonlinear  analysis  of  the  system 
being  studied.  It  also  allows  mathematical  models  suitable  for  optim¬ 
ization  studies  to  be  formulated  with  relative  ease  (as  far  as  obtaining 
the  required  differential  expressions) . 

The  model  developed  in  this  report  is  completely  nonlinear.  The 
reasons  for  this  approach  and  an  in-depth  discussion  of  its  derivation 
are  given  in  the  sections  that  follow. 

1.2  Description  of  the  Model 

The  purpose  of  this  technical  report  is  to  document  the  work 
accomplished  to  date  on  the  development  of  a  mathematical  model  for 
the  Howitzer,  Heavy,  Self-Propelled,  Full-Tracked,  8-inch,  MllOAl, 
see  figure  1.1.  This  documentation  has  two  major  objectives;  (1)  the 
detailed  description  of  a  new  method  for  developing  and  solving 
Lagrange’s  equation  of  motion  and  (2)  the  generation  of  a  mathema¬ 
tical  simulation  to  describe  the  motion  of  the  system  and  to  define 
dynamic  loading  at  various  points  in  the  weapon  (which  are  to  be  used 
as  input  for  a  stress  analysis  of  weapon  components) . 

A  method  was  developed,  utilizing  FORMAC,  to  obtain  the  necessary 
symbolic  representations  for  the  differential  expressions  required  to 
solve  Lagrange’s  equations  of  motion  of  the  MllOAl.  FORMAC,  ^.s 
developed  by  IBM,  provides  for  the  symbolic  manipulation  of  mathema¬ 
tical  expressions,  i.e.,  the  expression  SIN(X)  can  be  differentiated, 
resulting  in  the  expression  COS(X).  Expressions  can  be  differentiated, 
evaluated,  replaced,  compared,  and  parsed.  After  differentiation, 
expressions  which  occur  repeatedly  can  be  replaced  by  new  variable 
names  and  thus  millions  of  arithmetic  operations  can  be  eliminated 
during  the  execution  of  the  computer  program.  This  is  illustrated 
in  appendix  E .  Since  PL/ I  is  a  subset  of  FORMAC,  all  of  the  facilities 
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of  PL/I  are  available  for  program  structure,  loop  control  and  I/O. 

The  FORMAC  output  consists  of  coded  differential  expressions  and  they 
are  automatically  punched  on  cards  (error  free)  in  FORTRAN  format.  The 
end  product  is  a  computer  program  written  in  FORTRAN.  Approximately 
95  percent  of  the  work  to  generate  the  simulation  model  was  accomplished 
on  the  computer.  This  procedure  essentially  reduced  the  11  degree-of- 
freedom  problem  described  below  to  a  rather  routine  operation. 

In  the  development  of  mathematical  models  of  weapon  systems  it 
is  desirable  to  limit  the  degrees  of  freedom  of  the  model  to  the  major 
gross  motions  of  the  system.  In  general,  the  larger  the  model  (in 
terms  of  degrees  of  freedom)  the  less  accurate  is  the  predicted  motion 
of  the  system.  This  is  primarily  due  to  the  difficulty  in  accurately 
defining  the  values  of  the  weapon  parameters.  However,  this  problem 
can  be  overcome  by  performing  carefully  instrumented  tests  (both  static 
and  dynamic) .  If  a  physical  system  does  exhibit  significant  multi- 
degree  of  freedom  motion  it  is  necessary  to  model  that  motion  and  accept 
the  complexity  of  the  resulting  system  of  equations.  In  this  report 
it  was  essential  that  all  major  motions  be  identified,  resulting  in  an 
eleven  degree-of-freedom  system. 

The  model  developed  for  the  MllOAl  has  five  distinct  masses: 

(1)  the  vehicle,  M^,  (2)  the  spade  assembly,  (3)  those  parts  that 

traverse  but  do  not  elevate,  M^,  (4)  those  parts  that  elevate  but 
do  not  recoil,  M^,  and  (5)  the  recoiling  parts. 

The  generalized  coordinates  were  defined  based  on  the  following 
logic:  the  entire  system  would  translate  laterally  (x) ,  roll  (0), 

and  yaw  (\^)  as  a  single  rigid  mass;  the  vehicle  would  translate  fore 
and  aft  (y) ,  translate  in  a  vertical  direction  (z),  and  pitch  ((()) 
about  its  own  mass  center;  those  parts  that  traverse  (M^  +  + 

would  yaw  (t)  relative  to  the  vehicle;  those  parts  that  elevate 
(Mg  +  Mg)  would  pitch  (y)  relative  to  the  vehicle;  the  recoiling  parts 
would  translate  (n)  relative  to  the  Mg  parts;  finally  the  spade  as¬ 
sembly  would  translate  fore  and  aft  (v)  and  pitch  (v)  relative  to  the 
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ground.  Note  the  initial  value  of  y  (Yq)  is  the  angle  of  elevation, 
the  initial  value  of  x  (x^)  is  the  angle  of  traverse,  and  the  initial 
value  of  n  (rijj)  locates  the  in  battery  position  of  the  mass  center 
of  the  recoiling  parts  relative  to  the  trunnion. 

A  brief  description  of  the  operation  of  the  physical  system  is 
given  below.  When  the  weapon  is  fired,  the  pressure  generated  by  the 
burning  propellant  drives  the  projectile  out  of  the  tube  and  forces 
the  recoiling  parts  (tube,  recoil  rods,  recoil  pistons,  etc.)  rearward. 
The  motion  of  the  recoiling  parts  is  primarily  resisted  by  oil  pressure 
on  the  face  of  the  recoil  piston.  The  piston  is  pulled  through  a  cyl¬ 
inder  of  oil,  the  oil  being  throttled  around  the  piston  via  orifice 
areas  which  are  designed  as  a  function  of  recoil  displacement  to  min¬ 
imize  peak  recoil  forces.  In  parallel,  to  the  recoil  cylinder  is  the 
recuperator  cylinder  wherein  the  recuperator  piston,  during  recoil, 
compresses  a  gas  and  hence  stores  enough  energy  to  return  the  recoil¬ 
ing  parts  to  the  firing  (In-battery)  position.  The  total  recoil  force 
is  transferred  to  the  under structure  through  the  trunnions.  As  a 
result  of  this  system  of  forces,  a  large  torque  is  generated  around 
the  trunnions  and  the  elevating  parts  will  tend  to  pitch  relative  to 
the  vehicle.  The  restraint  to  this  pitching  motion  is  offered  by  a 
friction  brake.  Until  the  moment  tending  to  produce  pitch  exceeds 
the  restraining  moment  of  the  friction  brake,  no  pitch  motion  results. 
When  the  moment  tending  to  produce  pitch  exceeds  the  restraining 
moment,  pitch  motion  is  initiated;  and  will  continue  until  the  con¬ 
stant  restraining  moment  of  the  brake  brings  the  pitching  motion  to 
a  stop.  A  similar  action  occurs  in  the  yaw  motion  of  the  traversing 
parts . 

The  spades  are  designed  to  offer  ground  resistance  to  rearward 
motion  and  to  the  pitching  of  the  spade  assembly  but  little  re¬ 
sistance  to  roll,  yaw  and  translations  in  the  lateral  and  vertical 
directions.  Thus  horizontal  ground  springs  acting  on  the  spade  as¬ 
sembly  in  the  longitudinal  direction  are  sufficient  to  restrain  the 
longitudinal  translation  and  pitch  of  the  spade  assembly . 
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The  resistance  to  lateral  translation  and  yaw  of  the  weapon  system 
is  the  ground  friction  between  the  ground  and  the  vehicle  tread  (assumed 
to  be  at  the  ground  contact  point  of  the  four  corner  roller  wheels) . 
Vertical  ground  springs  are  located  at  these  four  ground  contact  points 
and  restrain  the  roll  and  vertical  translation  of  the  weapon  system. 

Note  the  effective  spring  rate  goes  to  zero  at  any  of  the  ground  contact 
points  when  contact  between  the  ground  and  the  vehicle  is  lost.  The 
braces  and  spade  cylinders  are  modeled  as  springs  and  hydraulic  damping 
and  act  primarily  in  the  longitudinal  direction  between  the  vehicle 
and  the  spade  assembly.  Dampers  are  associated  with  the  vertical  ground 
springs  and  the  braces. 

The  method  of  mathematically  describing  the  above  physical  phe~ 
nomina  is  contained  in  Appendix  A. 


8 


2 . 0  SOLUTION  TECHNIQUE 

2.1  Overview 

Experience  lias  demonstrated  that  modeling  multi“degree  of  freedom 
systems  requires  many  hours  of  tedious  manipulation  of  expressions  ’with 
the  risk  of  generating  numerous  errors.  The  result  is  generally  the 
implementation  of  a  model  which  is  a  linearized  version  of  the  system 
under  consideration.  That  is,  in  the  conventional  technique,  the  analyst 
is  forced  to  linearize  most  of  the  generalized  coordinates  to  achieve 
simplification  of  expressions.  These  shortcomings  led  to  the  development 
of  a  semi-automated  procedure. 

This  report  develops  a  new  technique  for  obtaining  the  dynamic 
equations  of  motion  of  any  system  resulting  from  Lagrange's  equation. 

The  equations  generated  by  the  procedure  given  here  are  completely 
nonlinear.  This  approach  is  taken  for  several  reasons.  First,  the 
accuracy  that  might  be  lost  by  linearizing  due  to  cross— coupling  effects 
is  not  known;  secondly,  the  nonlinear  model  is  much  easier  to  obtain 
than  the  linearized  version  using  the  methods  developed  in  this  report 
and  the  nonlinear  approach  does  not  require  much,  if  any,  additional 
core  storage;  and  thirdly,  a  more  accurate  model  would  produce  results 
which  are  closer  in  agreement  to  that  of  the  real  world. 

It  was  found  that  the  term  -9 (Kinetic  Energy) / 3 (Generalized  Co¬ 
ordinates)  does  not  have  to  be  calculated  as  the  positive  of  this  term 
appears  in  the  expression  3 (Kinetic  Energy)/ 3 (Generalized  Velocities) 

and  they  cancel.  This  result  is  always  the  case.  If  the  procedure  in 
this  report  is  used,  two,  three,  and  perhaps  even  four  degree-of -freedom 
systems  can  be  accomplished  quite  easily  by  hand  if  the  FOEMAC  software 
package  is  not  available  to  the  analyst. 

2.2  Lagrange's  Equation 

The  expression  for  the  Lagrangian  method  which  yields  the  equations 
of  motion  for  a  dynamical  system  is 
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(2.1) 


_d  3(KE) 
dt 


3(KE) 


3(DE) 


3q 


34 


3  (PE) 
3q. 


where  j  =  1,2, . • . . ,k 

KE  =  Total  kinetic  energy 
DE  =  Total  dissipative  energy 
PE  =  Total  potential  energy 
Fj  =  Generalized  external  force 
q^  =  Generalized  coordinate 


q^  =  Generalized  velocity 
t  =  Independent  variable,  time 
k  =  Number  of  generalized  coordinates 

Equation  (2.1)  in  matrix  form  can  be  written  as 


A(q)q  =  B(q,q,t)  (2.2) 

The  first  term  of  equation  (2.1)  will  generate  the  A  matrix  plus  addition¬ 
al  terms  which  contribute  to  the  B  vector.  The  remaining  four  terms  of 
equation  (2.1)  make  up  the  rest  of  B.  It  will  be  seen  in  the  sequel  that 
equation  (2.1)  can  be  expressed  conveniently  in  matrix  form. 

In  generating  the  desired  equations  of  motion  by  utilizing  FORMAC, 
it  becomes  necessary  to  examine  Lagrange's  equation  in  detail,  see 
Appendix  B.  This  is  required  because  FORMAC  performs  only  partial 
differentiation.  The  results  of  the  analysis  in  the  appendix  determine 
which  derivatives  are  to  be  taken  and  how  they  are  ultimately  combined 
to  yield  a  set  of  2nd  order  nonlinear  differential  equations  describing 
the  dynamic  motion  of  the  system  under  investigation. 

In  the  conventional  method  of  actually  carrying  out  the  analysis 
of  the  Lagrangian,  an  imbedding  of  terms  occurs  within  a  given  expres¬ 
sion  (terms  occur  repeatedly)  due  to  the  matrix  operations,  dot  pro¬ 
ducts,  and  the  way  in  which  each  succeeding  mass  center  is  located 
from  the  proceeding  one.  To  avoid  the  imbedding  problem  as  much  as 
possible,  the  philosophy  of  the  technique  reported  here  is  to  operate 
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on  the  position  vectors  of  each  mass  before  dot  products  are  actually 
taken.  In  fact,  because  of  the  way  each  succeeding  mass  center  is 
located,  the  position  vectors  themselves  can  be  broken  up  into  smaller 
expressions. 

However,  imbedding  still  occurs  to  some  extent,  but  the  expressions 
are  small  enough  such  that  the  imbedded  terms  can  be  identified  rapidly 
and  replaced  by  new  variable  names  quite  easily  with  FORMAC.  Thus, 
millions  of  arithmetic  operations  are  saved  in  the  course  of  actually 
solving  the  differential  equations  during  the  running  of  the  computer 
program.  To  see  this,  the  reader  should  refer  to  Appendix  E  to  get  some 
idea  of  how  much  of  a  reduction  in  the  number  of  arithmetic  operations  can 
actually  be  accomplished  by  the  replacement  of  repeated  terms  with  new 
variable  names.  All  of  the  partial  differentiation  as  required  by 
Appendix  B  is  performed  on  the  smaller  expressions  and  then  a  reduction 
in  the  size  of  these  differential  expressions  is  obtained  by  removing 
the  imbeddedness.  The  dot  product  of  the  position  vectors  is  taken  at  a 
later  time  with  numbers  as  elements  instead  of  with  large  expressions  as 
is  normally  done  in  the  conventional  method.  FORMAC  necessitated  the 
analysis  covered  in  Appendix  B  since  only  partial  derivatives  could  be 
obtained  by  using  FORMAC.  The  outcome  of  this  analysis  provides 
a  technique  that  facilitates  the  derivation  of  the  equations  of 
motion  whether  done  by  hand  or  with  the  aid  of  a  computer. 

2,3  FORMAC  Approach 

An  objective  of  the  work  reported  here  was  to  develop  an  organized 
and  efficient  computational  scheme  that  would  handle  multl-degree  of 
freedom  problems  with  relative  ease.  FORMAC  provides  this  capability. 

It  is  an  IBM  software  package  which  allows  the  manipulation  of  mathema¬ 
tical  expressions.  Using  FORMAC,  expressions  can  be  differentiated, 
evaluated,  replaced,  compared,  and  parsed.  Since  PL/l  is  a  subset  of 
FORMAC,  all  of  the  facilities  of  PL/l  are  available  for  program  struc¬ 
ture,  loop  control  and  I/O.  The  computer  performs  the  necessary  matrix 
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multiplication  and  differentiation.  In  addition,  angles  are  coded,  ex~ 
pressions  are  optimized  and  finally  punched  in  FORTRAN  format  ready  for 
numerical  integration.  This  process  virtually  eliminates  any  mathematical 
or  key  punching  errors  for  approximately  85  percent  of  the  FORTRAN  pro¬ 
gram  as  this  amount  of  the  computer  program  is  obtained  from  the  results 
of  the  FORMAC  output.  This  technique  also  allows  the  analyst  to  spend 
considerably  more  time  formulating  the  problem  and  not  be  worried  about 
the  enormous  amount  of  mathematics  required  to  obtain  the  desired  equa¬ 
tions  of  motion.  In  short,  it  reduces  problems  which  seemingly  appear 
hopeless  to  a  rather  routine  operation. 

2.4  Kinetic  Energy 

This  section  covers  the  kinetic  energy  portion  of  the  Lagrangian. 

To  be  more  specific,  it  discusses  the  use  of  equation  B-11.  An  ex¬ 
ample  is  given  on  how  to  obtain  the  A  matrix  of  equation  (2.2)  by  use 
of  B-11  and  also  the  corresponding  FORTRAN  coding  generated  by  the 
Einstein  summation  notation  is  shown. 

Since  nearly  all  of  the  required  matrix  operations  for  a  particular 
problem  are  contained  in  the  definition  of  the  kinetic  energy ,  this 
becomes  the  obvious  starting  point  for  obtaining  the  necessary  differ¬ 
ential  expressions.  It  will  be  seen  in  later  sections  that  many  deriv¬ 
atives  for  other  energy  terms  do  not  have  to  be  calculated  as  this  will 
already  have  been  accomplished  in  the  kinetic  energy  portion. 

As  mentioned  previously,  it  is  more  efficient  (computationally)  to 
operate  on  the  position  vectors  before  dot  products  are  actually  taken. 

And  also  mentioned  before  was  that  the  position  vectors  themselves  are 
not  really  operated  upon  since  they  are  in  fact  reduced  to  sums  of 
smaller  quantities  (to  reduce  the  imbedding  problem)  and  it  is  these 
smaller  quantities  which  are  of  interest.  The  actual  break-up  of  the 
position  vectors  is  covered  in  section  3.3.  Because  of  the  ease  in 
actually  combining  the  smaller  terms  to  obtain  the  position  vectors, 
the  example  given  here  will  be  concerned  with  only  the  position  vectors . 
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Equation  B~ll  is  rewritten  for  convenience. 


dt  3q 


(2.3) 


where  all  quantities  are  defined  in  Appendix  B.  mass  asso¬ 

ciated  with  the  position  vector  Q,  where  Q  takes  on  the  various  letters 
Q,  R,  S,  etc.,  to  distinguish  each  of  the  n  position  vectors  and  its 
corresponding  mass. 

It  is  seen  from  equation  (2.3)  that  first  and  second  partial  deriv¬ 
atives  with  respect  to  the  generalized  coordinates  are  required.  The 
subscripted  variable  PKE  (i,j,k,l)  provides  a  convenient  bookkeeping 
notation  for  representing  these  derivatives,  the  mass  from  which  they 
originate  and  the  direction  (x,y,  or  z)  of  the  position  vector  with 
respect  to  the  inertial  coordinate  system.  From  past  experience  this 
bookkeeping  procedure  has  proven  to  be  very  adequate;  it  is  organized, 
efficient  and  yet  simple  and  it  also  conveys  a  maximum  amount  of 
information  such  that  a  particular  quantity  can  be  easily  identified. 

For  example,  the  i  runs  from  1  to  n  and  refers  to  the  particular  mass 
and  associated  position  vector  that  is  being  delt  with;  the  1  runs  from 
1  to  3  and  refers  to  the  x,y,  or  z  direction;  the  j  and  k  refer  to  the 
partial  derivatives  which  have  been  taken  with  respect  to  the  q^  gen¬ 
eralized  coordinate  and  the  q^  generalized  coordinate  if  second  partials 
have  been  taken.  If  only  one  derivative  has  been  taken,  then  k  takes 
on  the  value  1  plus  the  number  of  generalized  coordinates. 

Suppose  the  number  of  generalized  coordinates  is  eleven  and  the 
position  vectors  are  Q,  R,  and  S,  then 

^  '  =  PKE  (1,5,8,L),  where  L  =  1,2,3 

9q53q8 

9^5  =  PKE  (3 ,5,8,L) ,  where  L  =  1,2,3 

sqjaqj 
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PKE  (2,9,12,L),  where  L 


1,2,3 


An  example  of  using  equation  (2.3)  with  the  number  of  masses  n  =  2 
and  the  number  of  equations  k  =  3  and  considering  only  the  acceleration 
terms  of  (2.3)  yields  the  following  result. 


First  equation: 


+  M. 


R 


Q  •  Q  +  Q  •  Q  ^9  +  Q.  * 

qi  1  q2  2  q^  q3  3 


^  S  q,  +  R  •  ^  q„  +  1  ^  <i 

qi  q^  1  q^  q2  2  q^^  q3  3 


Second  equation 


M, 


Q  •  Q  *4.  +  Q  -^40  + 

q2  qi  1  q2  q2  2 


^2  ^3 


+  M 


R 


R  •  R  d-  +  R  •  R  q«  +  R  •  R  q« 

q2  q^  1  q2  q2  ^  q2  q3  3 


Third  equation: 


M, 


Q  *  Q  q,  +  Q  •  ^  4  4  <io 

^^3  ^1  ^  ^3  ^2  ^  ^3  ^3  ^ 


+  M 


R 


-y  -y 


R  •  R  q_,  +  R  •  R  q„  +  R  •  R  q. 
qS  ^  ^  *^3  *^3  ^ 
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These  three  equations  can  be  rewritten  in  the  form 


First  equation: 


M  Q  •  $  +  •  R  q,  +  M  ‘ 

Q  q^  R  \  1  Q  'll  'I2 

\\  •  \  «2  “q\  •  +  \\  •  *,3  ‘ 


Second  equation: 


-y 


Q  "q2  ^  ^2  ^1 


R  q2  ^2  ^ 


2  '‘a 


->■ 

M.  Q 

Q 

Q  q2 

q 

+ 

t 

R  q2 

q 

Third  equation: 


M  ^  t  ^  111  + 

Q  q3  q^  R  ‘I3  'll  ^  ^  ^3  ^2 

MR  •  R  a+MQ  *0+  MR  •  R  <i 

”r  q  q  ^2  Q  q-i  qo  R  q-j  'i-i 

32  3  3  J 

J  L- 

The  first  row  of  the  A  matrix  of  equation  (2.2)  is 

A(1  1)=MQ  *0  +M„^ 

AU,i;  '^q^  ^  q^  q^ 


A(1.2)  =  •  t 


—  L  Q^qi  ^^2  ^  ^2 

A(l,3)  -  M^Q  •  %  +Wk\  •  *,3 


and  so  forth  for  the  remaining  six  elements  of  A  of  this  example 


15 


The  FORTRAN  Coding  for  obtaining  the  elements  of  the  A  matrix  for 
the  translational  part  of  the  kinetic  energy  is  now  given.  The 
rotational  part  must  be  added  to  these  terms  and  this  can  be  found  in 
Appendix  F  which  contains  the  FORTRAN  program.  Similarly  the  remain~ 
ing  coding  for  the  kinetic  energy  is  in  appendix  F. 

Since  A  is  symmetric,  the  upper  triangular  terms  are  calculated 
to  reduce  the  number  of  arithmetic  operations.  Recalling  that  lEQS  =  3 
and  IMASS  =  2  yields 

DO  1  j  =  l.IEQS 

DO  1  k  =  j.IEQS 


SUM  =  0. 

DO  2  i  =  1, IMASS 

DO  2  1  =  1,3 

2  SUM  -  SUM  +  XMASS(i)*  PKE(i, j , 12 , 1) *  PKE(i,K,i2,l) 


A(j,k)  =  SUM 
1  A(k,j)  =  SUM 

2.5  Other  Energy  Terms 

The  differential  expressions  required  for  other  energy  terms  are 
either  obtained  from  those  which  have  already  been  calculated  or  are 
derived  separately.  As  an  example  of  how  to  utilize  expressions 
already  calculated,  the  reader  can  refer  to  equations  A-7  through  A-11 
in  Appendix  A.  These  five  equations  are  very  similar  to  equation  A-12 
except  that  A-12 considers  only  the  z  direction  and  involves  both 
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multiplication  by  the  gravitational  constant,  g  and  by  the  five  masses. 
Therefore,  3Uj^/9q.  can  be  obtained  very  quickly  from  the  following 
definitions  (which  are  explained  in  section  3.3).  For  equations 


The  3Uj^/3q.  makes  up  a  part  of  the  B  vector  in  equation  (2.2)  which 
forms  the  right  hand  sides  of  the  differential  equations.  The  con¬ 
tribution  of  9Uj^/9qj  (from  the  potential  energy)  to  the  right  hand 
sides  (utilizing  calculations  from  the  kinetic  energy)  is  evaluated 
as  follows 

“  do  14  J  =  1,IEQS 

14  RHS(J)  =  RHS(J)  +  GRAV*(XMASS(1)*PKE(1,J,12,3)  +  XMASS(2)* 


1  PKE(2,J,12,3)  +  XMASS(3)*PKE(3,J,12,3)  +  XMASS(4)*PKE(4,J,12,3) 


2  XMASS(5)*PKE(5,J,12,3)) 


Note  that  the  fourth  subscript  of  PKE  is  equal  to  3  which  signifies  that 
only  those  terms  in  the  z  direction  are  used.  The  third  subscript  is 
equal  to  twelve  since  only  one  partial  derivative  is  taken  (the  number 
of  generalized  coordinates  is  eleven) . 

As  much  use  as  possible  is  made  of  all  previously  defined  matrix 
operations  or  whatever  caculations  have  been  performed  to  aid  in  obtain¬ 
ing  new  partial  derivatives.  This  can  readily  be  observed  when  examin¬ 
ing  the  FORMAC  program. 

2.6  Numerical  Integration 

A  fourth  order  Runge-Kutta  integration  scheme  is  used  to  integrate 
the  matrix  differential  equation 

A(q) q  =  B(q,q,t) 

To  solve  for  the  q,  i.e,, 

^2  ° 


advantage  is  taken  of  the  symmetry  in  the  mass  matrix  A.  The  sub¬ 
routine  which  decouples  the  acceleration  terms  into  the  form  of  equations 
(  2.4  )  is  called  SOLVE  and  the  algorithm  used  in  the  Square  Root  Method. 
Appendix  D  describes  the  modifications  that  were  made  to  that  method  to 
eliminate  the  concern  of  pure  imaginary  numbers.  A  description  of  each 
subroutine  and  its  purpose  is  given  in  the  FORTRAN  listing  in  Appendix  F. 
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3*0  APPLICATION 
3. 1  Symbol 

Computer  Symbol 
Al,  A2,  A3 

AAA  (11,  11) 

AISUB,  (AIBAR), 
A2SUB,  A3SUB 

ALPHAl,  ALPHA2, 
ALPHA3 

ASTAR 

BETA 


BBETAE,  BETAE 

BBETAT,  BETAT 

BO  FT 


OF  SOLUTION  TECHNIQUE 
Table 


Corresponding 

Analysis  Symbol  Definition  Unit 

a-,  a«,  a^  Coordinates  of  the  mass  inches 

center  of  vehicle  in  0^- 
A**B*’C*' 


(A^),  A^, 


Coefficients  matrix  of 
acceleration  (mae*  Matrix) 
terms 


pound  second 
inches  ^ 
or  lb. sec.  in. 


Coordinates  of  the  attach¬ 
ments  of  the  brace  to  the 
vehicle  in  0,  - 


inches 


^1’  ^2*  ^3 


6. 


Coordinates  of  the  horizon-^-*  Inches 

tal  attachment  points  of 
ground  springs  to  spade  as¬ 
sembly  in  0^  -  ABC 

Angle  between  of  brace  degrees 

and  a  plumb  line  dropped  from (changed  to 
the  attachment  of  the  brace  radians  in 
to  the  vehicle  program) 

a  A^  2 

3  =  —t:  where  j  *  specie  pound  sec 

^  Aq  inches^ 

fic  weight  of  oil,  A  *  area 
of  piston  and  Aq  »  effective 

area  of  orifice.  6  is  a 
hydraulic  constant  for  the 
spade  cylinder 


^E’  ^EMAX 


torque  provided  by  friction 
brake  in  elevating  mech¬ 


anism. 


is  maximum. 

tMAA 


inch  pounds 


^tMAX  torque  provided  by  friction  inch  pounds 

brake  in  traversing  mecht^ 
anism.  i®  niaximum. 


B<t) 


Breech  force 


pounds 
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Computer  Symbol 


BRCHX  (105) 


BRCHY  (105) 


Bl,  (BIBAR) 
B2,  B3 


CBRCE 


cell,  CC12, 
CC21,  CC22 


COFT 


CU5 


CIP,  C2P,  C3P 


DA  (1,  1,  1,  L) 
L=  1,2, 3, 4 
l-(->(l,l) 
2<->(l,2) 

3^(2, 1) 

4^(2, 2) 


DP  (1,  1,  1,  4) 


Corresponding 


Analysis  Symbol 

Definition 

Unit 

Abcissa  of  breech  force 
curve (time) 

seconds 

Ordinate  of  breech  force 

curve 

pounds 

(B^),  B3 

Coordinates  of  the  attach¬ 
ments  of  the  spade  cylin¬ 
ders  to  the  vehicle  in 

0  “  X"Y"Z" 

4 

inches 

C, 

1 

Damping  coefficient  asso¬ 
ciated  with  K. 

1 

pound  sec 
inch 

C.  . 

iJ 

Damping  coefficients  asso¬ 
ciated  with  K. . 

ij 

pound  sec 
inch 

C(t) 

Recuperator  force 

pounds 

Damping  coefficient  for 
use  in  determining  initial 
conditions 

in,  lbs.  sec 

pt  pt  pT 

Cl,  C2,  c^ 

Coordinates  of  the  attach¬ 
ments  of  the  spade  cylin¬ 
ders  to  the  spade  assembly 
in  0^  -  U’V’W’ 

inches 

Extension/contraction  of 
vertical  ground  spring 

inches 

^'0 

iV\ 

Spring  deflections ) 

pounds 

inch 
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Computer  Symbol 
Dl,  D2,  D3 


DIP,  D2P,  D3P 


EB 


EC 


ER 


ELANG 


El,  E2,  E3 


FOFG 


Corresponding 
Analysis  Symbol 


Definition 


Coordinates  of  the  center 
of  the  trunnion  in 
0.  - 


Unit 

inches 


Coordinates  of  the  attach¬ 
ments  of  the  braces  to  the 
spade  assembly  in  0^  -  U'V’W’ 


inches 


Distance  from  mass  center  inches 
of  recoiling  parts  to 
point  of  application  of 
breech  force  in  n  direction 
only,  t  is  negative. 


Distance  from  mass  center  inches 
of  recoiling  parts  to  point 
of  application  of  recuperation 
force  in  n  direction  only, 
is  negative. 


Distance  from  mass  center  inches 
of  recoiling  parts  to  point 
of  application  of  recoil 
force  in  n  direction  only. 
e  is  negative. 


Angle  of  elevation  for  use  radians 
in  determining  initial  con¬ 
ditions  of  Q(7) 


®1’  ®2’ 


Coordinates  of  the  center  of  inches 
the  traverse  bearings  in 


0  .  x"Y”Z" 

4 


f(y) 


Equilibrator  force  as  a  pounds 

function  of  y* 
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Corresponding 

Computer  Symbol  Analysis  Symbol  Definition 


Unit 


FFl,  FF2,  FF3 


GAMMAX  (105) 


GAMMAY  (105)  F(y) 

GKST  K 

St 


Coordinates  of  the  mass  inches 

center  of  the  traversing 


Abscissa  of  equllibrator  radians 

force  curve  data  in  degrees 
but  converted  to  radians  in 
program. 

Ordinate  of  equllibrator  pounds 

force  curve. 


Effective  torsional  spring 
about  trunnion  used  in  de~ 
terminiag  initial  cond It Idas 


inch  pounds 
radians 


GRAV  g 


Acceleration  due  to  gravity  inches/sec 


Gl,  G2,  +  Q(2),  g-  g2  +  V, 
G3  g3 


Coordinates  of  the  center  inches 
of  pressure  of  the  spades  in 
0  -  XYZ.  Note  v  is  a  gener¬ 

alized  coordinate. 


HHl,  HH2,  HH3  h  ,  h^,  h  Coordinates  of  the  mass  cen-  inches 

ter  of  the  spades  in 
0^  -  U’V’W’. 


HOORl 


Delta  function  used  in 
defining  initial  conditions 


IBOFT 


table  lookup  counter  on 
breech  force 


IBPTS 


No.  of  ordered  pairs  in  breech 
force  table 


lEQS 


Number  of  generalized  coordinates 


IGOFT 


Table  lookup  counter  on 
equllibrator 


Corresponding 

Computer  Symbol  Analysis  Symbol  Definition 


IGPTS 


No*  of  ordered  pairs  of 
points  in  equilibrator  table 


IMASS 


Number  of  masses 


IROFT 


Table  lookup  counter 
on  recoil  force. 


IRPTS 


No.  of  ordered  pairs  of 
points  in  recoil  force  table. 


02,  03 


02. 03 


Coordinates  of  the 
equilibrator  attachment 
points  to  the  elevating 
but  non-recoiling  parts 
in  0,  -  H’H'Z’ 


inches 


PD  (3,  11,  12,  4) 


PG  (6,  12,  12,  3) 


Derivative  functions  in 
dissipative  energy 

Derivative  function  in 
generalized  forces 


PT  (8,  11,  12,  3) 


Derivative  function  in 
translational  part  of 
kinetic  energy 


PU  (4,  11,  12,  4) 


PW  (5,  12,  12,  3) 


Derivative  functions  in 
potential  energy 

Derivative  functions  in  ro¬ 
tational  part  of  kinetic  energy 


QDD  (I) 


Accelerations  of  the 
generalized  coordi¬ 
nates  I  =  1, ... ,11 


Units  of  Q(i) 
second^ 


QD(I) 


Velocities  of  the  gen¬ 
eralized  coordinates, 
I  =  1,...,11 


Units  of  Q(i) 
second 


QSAVE  (11) 


Dimension  variable  which 
saves  the  generalized  coor¬ 
dinates  in  Runge-Kutta  integration 


Corresponding 

Computer  Synibol  Analysis  Symbol  B e f ini t i Oji 


QDSAVE  (11) 


Q(l) 


n 


Q(2) 


V 


Q(3)  X 

Q(4)  y 


Dimension  variable  which 
saves  generalized  velocities 
in  Runge-^Kutta  integration. 

Translation  of  re¬ 
coiling  parts.  Has 
initial  value. 


Fore  (+)  and  aft  (-) 
motion  of  spade  as¬ 
sembly. 


Lateral  displacement 
of  total  weapon  system 


Fore  (+)  and  aft(-) 
motion  of  vehicle  and 
traversing  parts 


Q(5) 


z 


Up  (+)  and  down  (-)  motion 
of  vehicle  and  tra-? 
versing  parts 


Q(6) 


Pitch  of  vehicle  and 

traversing  p^arts 


Q(7)  Y 


Q(8) 


V 


Pitch  of 
relative 
May  have 
angle  of 


elevating  parts 
to  the  vehicle, 
initial  value 
elevation  - 

o 


Pitch  of  spade  assembly 


Q(9)  e 


Roll  of  total  weapon  system 


Q(10)  If; 


Yaw  of  total  weapon  system 


Unit 


inches 


inches 


inches 


inches 


inches 


radians 


radians 


radians 


radians 


radians 
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Corresponding 

Computer  Symbol  Analysis  Symbol  Definition 

Q(ll)  T  Yaw  of  traversing 

parts  relative  to  the 
vehicle.  May  have 
Initial  value  (angle 
of  traverse  “  t) . 


RHS  (11) 

RODY  (105) 

RODX  (105) 

ROFT  R(t) 

TIME  t 

TIMER 

TIMEH2 


Right  hand  side  of  equa¬ 
tions  of  motion. 

Ordinate  of  recoil  force 

Abclssa  of  recoil  force 
(time) 

Recoil  force 

Time  measured  from  Ini¬ 
tiation  of  solution. 

Integration  step  size 

Defined  as  TIMER/ 2. 


TIMERS 


Defined  as  TIMEH/S 


XI,  Q(l),  ZETA  Coordinates  of  the  mass 

center  of  the  recoil¬ 
ing  parts  In  0^  -  E*H’Z’ 
Note  n  is  generalized  co¬ 
ordinate. 


XIB,  Q(l)  +  EB,  H  ,  rig,  Cg  Coordinates  of  point  of 

ZETAB  application  of  breech 

face  In  0,  -  E^R’Z^ 

D 


Unit 

radians 


pounds 

seconds 

pounds 

sec 

sec 

sec 

sec 

inches 

inches 
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Computer  Symbol 

XIC,  Q(l)  +  EC, 
ZETAC 


XIR,  Q(l)  +  ER, 
ZETAR 


XIXX(I) 


XIXY(I) 

XIXZ(I) 

XIYY(I) 


XIYZ(I) 

XIZZ(I) 


Corresponding 

Analysis  Symbol  Definition 


^C’  "^C*  ^C 


’^R 


XII,  ETAl,  ZETAl 


Moments  of  inertia  of 
masses  about  their  prin¬ 
cipal  "Z"  axis  where  Q,  R, 
S,  T,  U  correspond  to 
I  =  1,2, 3, 4, 5 

Coordinates  of  the  mass 
center  of  the  elevated 
but  non-recoiling  parts 
in  Og  - 


Unit 


inches 


inches 


lbs  sec  in 


Coordinates  of  point  of 
application  of  recuperator 
force  in  0^  - 


Coordinates  of  point  of 
application  of  recoil 
force  in  0^  -  E'H’Z’ 

Moments  of  inertia  of 
masses  about  their  prin¬ 
cipal  ”X**  axis  where  Q,  R, 
S,  T,  U  correspond  to 
I  =  1,2, 3, 4, 5 


Cross  products  of  inertia  lbs  sec  ifi 


Cross  products  of  inertia  lbs  sec  1ft 


Moments  of  inertia  of  lbs  sec  in 

masses  about  their  prin¬ 
cipal  ”Y**  axis  where  Q,  R, 

S,  T,  U  correspond  to 
I  -  1,2, 3, 4, 5 

2 

Cross  products  of  inertia  lbs  sec  in 


lbs  sec  in 


inches 


XKKl,  XKK2 


K 


Effective  spring  rate  of  pounds/lnch 

one  of  the  braces 
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Computer  Symbol 


XKPll,  XKP12, 
XKP21,  XKP22 


XKYl,  XKY2 


XK(1,1),  XK(1,2) 
XK(2,1),  XK(2,2) 


XL  (1,1).  XM(1,1 
XN(1,1),  XL(1,2> 
XM(1,2),  XN(1,2) 
XL(2,1),  XM(2,1) 
XN(2,1),  XL(2,2) 
XM(2,2),  XN(2,2) 


XLENGH 


XMASS(l) 


XMASS(2) 


Corresponding 

Analysis  Symbol  Definition  Unit 


Horizontal  ground  spring  pounds/inch 

at  four  roller  wheels  to 
prevent  lateral  translation 


Spring  rate  of  horizon-  pounds/inch 

^  tal  ground  springs  asso¬ 

ciated  with  spade  assem¬ 
bly. 


Spring  rates  associated  pounds/ inch 

with  ground  springs  at 
front  and  rear  roller 
wheels 


’  hj’  ”‘ij’  “ij 

Coordinates  of  ground  inches 

springs  at  front  and  rear 

roller  wheels 

i  =  l(right)  i*2(left) 

j  «=  l(front  j»2(rear) 

in  0,  -  X"Y’*Z” 

4 

Length  of  equilibrator  inches 


M 

V 


Mass  of  the  vehicle 


pound  second 
inch 


M 

s 


Mass  of  the  spade  assem¬ 
bly 


^2 

pound  second 
inch 
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Computer  Symbol 


Corresponding 

Analysis  Symbol  Definition 


XMASS(3) 


Mass  of  the  traversing 
but  non-elevating  parts 


iQund  second 
inch 


XMASS(4) 


Mass  of  elevating  but  non-  pound  second 
recoiling  parts  inch 


XMASS(5) 


Mass  of  recoiling  parts 


>ound  second' 
inch 


Coefficient  of  friction 
between  vehicle  trdad  and 
ground 


XNN2,XNN3  N^, 


Coordinates  of  the  equili- 
brator  attachment  points 
to  the  traversing  but  non¬ 
elevating  parts  in 

0,  - 


inches 


ZZ(I)  I  «  1,2,,. .30 


Factored  algebraic  expressions 
to  reduce  the  number  of  arith¬ 
metic  operations 


28 


3.2  Technique  of  Generating  Energy  Expressions  and  Generalized 
Forces • 


The  techniques  used  in  the  development  of  Appendix  A  are  illus¬ 
trated  in  the  following  paragraphs. 

Assuming  a  vector  to  be  defined  in  one  coordinate  system,  it  is 
necessary  to  determine  a  coordinate  transformation  which  will  define 
that  vector  in  a  different  coordinate  system.  These  coordinate  trans¬ 
formations  are  needed  to  refer  velocity  and  displacement  vectors  to 
a  fixed  coordinate  system  and  angular  velocity  vectors  to  the  appro¬ 
priate  body  axes.  As  an  example,  three  coordinate  systems  are  defined. 
They  are:  (1)  Oj^-ABC,  an  inertial  coordinate  system  (fixed  in  space) 
having  its  origin  at  the  mass  center  of  the  weapon  system  in  the  em¬ 
placed  position;  (2)  02“A^B^C^,  a  coordinate  system  initially  coinci¬ 
dent  with  Oj^-ABC  and  remaining  parallel  to  Oj^~ABC  at  all  times  where 
0^  is  the  center  of  mass  of  the  system;  and  (3)  O^^A^B'^C’^  a  coordinate 
system  fixed  in  the  weapon  system  and  moving  with  it 


or 


A 

’ai' 

X 

B 

C 

s 

_ci_ 

+  where  X^  = 

0 

0 

A  *=  A^  +  X 
B  ==  b1 
C  = 


A^  =  A^  Cos  0  +  C  Sin  0 
fil  =  B^ 

=  “A^  Sin  0  +  C  Cos  6 


Fig.  3.2 
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thus 


and 


> 

_ j 

Cos  0  0  Sin  0 

s  Xa 

®i 

Xa  = 

0  10 

-Sin  0  0  Cos  0 

Hi*  ^ 

Fig,  3.3 


=  A"  Cos  -  B"  Sin 
■=  A"  Sin  ilJ  +  B"  Cos  ij) 


A 

‘  A"' 

Cos  ^  -Sin  i|J  0 

1 

=  ^b 

B" 

where  X^  = 

Sin  ^  Cos  ^  0 

_"l_ 

C" 

0  0  1 

- 

and  hence  by  successive  substitution 


"a” 

’a"* 

B 

=  Xi  +  Xz 

B" 

c 

C" 

where 


X  • 

a 


X 


b 


The  matricies  X^,  and  X2  are  orthogonal  and  hence  their  in¬ 
verses  are  their  transposes.  Now  if  the  coordinates  of  a  point  in 
the  weapon  system  are  known  in  02-A''B"C"  they  can  be  determined  in 
Oj^-ABC.  All  transformations  in  Appendix  A  were  determined  in  this 
fashion. 

The  absolute  angular  velocity  of  the  system  about  the  B  C 

coordinate  system  is  (see  Figure  3.1  and  3.2) 

0)  *  9  sin 

X 

0)  =  0  cos  ij; 

y 
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or 


0) 

■  o' 

X 

-1 

0) 

y 

II 

0 

• 

ca 

z 

A"B"C" 

X 

All  angular  velocities  were  obtained  by  similar  transformations. 
Defining  an  element  of  mass  in  the  weapon  system  to  be  dm  and  the 
vector  from  the  origin  of  the  fixed  coordinate  system  to  that  element 
of  mass  to  be  p,  then  the  kinetic  energy  (dt)  of  that  element  of  mass 
is 


dX  =  1/2  p  ‘p  dM 


Define 


p  =  P  +  P 

where  ?  is  the  vector  from  the  origin  of  the  fixed  coordinate  system 
to  the  mass  center  of  the  weapon  system  and  p  is  the  vector  from  the 
mass  center  of  the  weapon  system  to  the  element  of  mass  dm.  Now 

"  -V  •  • 

dX  =  1/2  (?  +  p)  •(?  +  p)dm 

Since  the  weapon  system  has  a  rigid  body  rotation,  then 
A- 

P  =  0)  X  p 


thus 


dX 


1/2  (^  +  ux^)*  (^  +  wx^)dm 


expanding  and  integrating 
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T=l/2MP*P  +  P»w  x  pdin+  (cjxp)  •  (a)Xp)dm 

M  M 


The  sum  of  the  first  moments  of  a  mass  about  its  own  mass  center  is 
zero.  Thus 


p  dM  =  0 


and  oj  X  p  can  be  written  as; 


and 


^  ^  (o)  X  u)  z)  +  k(a)  y  -  0)  x) 

(a3Xp)=:i(co2^a)y)+j  z  x^  y  ' 

y  2 


I  (oj  X  p  )  •  (o)  X  p )  dM  =0)^ 

J  X 

M 


M 


(y^  +  z^)  dM  +  0)^ 

y  J 


(x^  +  z^)  dM 


M 


+  0)' 


(x^  +  y^)  dM  -  2(jo  ^co  ^y  dM 

X  y  J 


M 


M 


“2co  CO  xzdM-2co  o) 

V  y  z  J 


y  z  d  M 


M  M 

By  definition  of  moment  of  inertia  and  cross  products  of  inertia 


1/2 


(co  X  ^)  .  (w  X  p )  dM  =  1/2  (I  0)^  +  I  cj^  +  I  0)^  ) 
^  XX  yy  zz 


M 


-(I  CO.  co+I  coco+I  coco) 
xy  X  y  xz  X  z  yz  y  z 


Neglecting  cross  products  of  inertia  the  kinetic  energy  of  the  weapon 
is 

T  =  1/2MP  •  P  +  1/2  (I  0)^  +  I  0)2  +  I  0)2) 

XX  y  y  z  z 

where 


and  the  o).  are  as  previously  defined.  I  ,  I  ,  I  are  the  principal  moments 
1  X  y  z 

of  inertia  of  the  system  about  the  axes. 

The  potential  energy  of  a  mass  is  simply  the  weight  of  that  mass 

times  the  Z-coordinate  of  its  mass  center  as  measured  in  the  fixed 

coordinate  system.  Define 

Xq  “  [001],  then  for  our  example 


Energy  stored  in  a  spring  having  a  spring  constant  K  is 
U2  =  1/2  K  (AL)2 

where  AL  is  the  extension/contraction  of  the  spring.  It  is  only  necessary 
to  define  L  as  a  function  of  the  generalized  coordinates. 

Associated  with  springs  in  a  system  is  damping.  Thus  for  the 
Spring  K  assume  a  damping  coefficient  C.  Then  the  dissipative  function 
for  damping  corresponding  to  the  energy  function  for  the  spring  is 

U2  =  1/2  C  (A'L)2 
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Assuming  a  forcing  function,  F(t)  to  be  applied  at  a^^,  a^, 
©2  -A"B"C"  then  in  0^  -  ABC  the  point  of  application  is 


The  direction  of  the  force  in  0^  -  ABC  is 


X 


2 


F^(t) 

FjCt) 

Fj(t) 


The  generalized  force  is  then  written  as 


F^(t) 

F2(t) 

F3(t) 


3qi 


X1  +  X2 


where  i  refers  to  the  particular  generalized  coordinate  q. 

There  are  four  basic  moments  acting  on  the  elevating  parts  around 
the  trunnion:  (1)  the  moment  due  to  the  weight  of  the  elevating  parts, 
(2)  the  moment  due  to  the  equilibrator ,  (3)  the  moment  due  to  the  breech 
force,  and  (4)  a  moment  due  to  the  friction  brake*  The  friction  brake 
will  offer  sufficient  torque  to  cancel  out  the  other  three  moments  and 
no  pitch  motion  occurs  until  the  algebraic  sum  of  the  other  three 
moments  exceeds  the  maximum  torque  that  can  be  generated  by  the  friction 
brake.  Then  the  brake  will  offer  a  constant  (maximum)  resistance  until 
the  pitch  velocity  becomes  zero.  This  is  the  argument  used  to  generate 
the  logic  for  the  mathematical  simulation  (in  Appendix  A)  of  the  friction 
brake. 
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3.3  FORMAC  Procedure 


This  section  describes  the  procedure  that  was  used  to  obtain  the 
necessary  partial  derivatives  to  solve  Lagrange's  equation  of  motion 
for  those  energy  expressions  derived  in  Appendix  A  by  utilizing  the 
results  of  Appendix  B.  The  translational  part  of  the  kinetic  energy 
will  be  given  first,  followed  by  the  remaining  energy  expressions. 

Equations  (A-7)  through  (A-11)  of  Appendix  A  can  be  obtained 
from  the  definitions  given  below  by  summing  combinations  of  the 
appropriate  expressions. 


PT(1,L)  =  ^3^  +  ^2^3 


PT(2,L)  = 


PT(3,L)  =  +  X^X^^  hi 


PT(4,L)  = 


PT(5,L)  =  A2A5X^  f^ 


PT(6,L)  = 


PT(7,1)  =^2^5^7^9  1 


PT(8,L)  =  ^2^5^7^9 


where  L  =  1,2,3 
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It  is  more  efficient  to  differentiate  these  eight  quantities 
rather  than  differentiate  the  five  original  position  vectors  Q,R,S,T 
and  U  as  repeated  terms  will  expand  the  derivatives  into  very  large 
expressions.  By  working  with  the  smaller  quantities,  the  analyst  is 
able  to  quickly  reduce  the  size  of  the  partial  derivatives  utilizing 
the  replace  statement  in  the  FORMA.C  program,  see  Appendix  E.  Also, 
computational  costs  are  higher  if  differentiation  is  performed  on  the 
original  expressions  as  considerably  more  core  storage  is  necessary 
and,  in  general,  the  expressions  just  become  too  large  to  handle 
efficiently.  The  entire  FORMAC  program  for  this  problem  required 
230K  bytes  of  core  using  TSO.  However,  in  order  to  keep  the  core 
size  at  a  minimum,  quantities  that  weren't  needed  for  certain  com¬ 
putations  were  commented  out. 

The  break-up  of  the  five  position  vectors  into  the  eight  smaller 
vectors  appear  exactly  the  same  here  as  in  the  FORMAC  program,  see 
Appendix  C.  After  the  definition  of  these  quantities  has  been  made, 
the  differentiation  can  proceed.  As  required  by  equation  B-11,  first 
and  second  partlals  are  to  be  taken.  Once  this  has  been  accomplished 
the  expressions  are  punched  in  FORTRAN  format .  The  first  and  second 
partials  of  these  eight  expressions  are  defined  by  the  subscripted 
variable  PT(1,J,K,L),  ...  ,  PT(8,J,K,L)  where  J  =  1,2,  ...  ,  11  and 
represents  the  partial  derivative  with  respect  to  the  J  generalized 
coordinate;  K  =  1,2,  ...  ,  12  and  represents  the  partial  derivative 
with  respect  to  the  generalized  coordinate  except  when  K  =  12. 

This  signifies  only  one  derivative  has  been  taken;  L  =  1»2,3  for  the 
x,y,  or  z  direction. 

Near  the  beginning  of  subroutine  NAME  in  the  FORTRAN  program,  the 
differential  expressions  PT(I,J,K,L)  are  combined  to  give  the  partial 
derivatives  of  the  original  position  vectors.  These  derivatives  have 
the  variable  name  PKE(I,J,K,L)  where  I  =  1,2, 3, 4, 5,  for  the  five 
position  vectors  and  the  J,K,  and  L  run  to  11,  12,  and  3  respectively. 
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The  next  energy  term  to  be  discussed  is  the  kinetic  energy 

(angular  part).  Equation  (B-14)  determines  the  derivatives  that  are 

needed  for  the  rotational  part  of  the  kinetic  energy.  Since  the 

inertia  terms  are  constant  in  the  A  and  B  vectors  of  equation  (B~14), 

all  that  is  necessary  at  the  moment  is  to  define  the  quantities 

W  (i) ,  W  (i) ,  and  W  (i)  where  i  =  1, 2,3,4, 5  for  the  five  masses  and 
X  y  z 

then  differentiate  these  terms.  This  procedure  begins  on  line  2390 
in  the  FORMAC  program  of  Appendix  C.  Because  equations  (A-5)  and 
(A«”6)  are  equal,  only  four  angular  vector  quantities  need  to  be 
differentiated.  The  definition  of  the  four  quantities  begins  on 
line  2650,  however,  the  development  starts  on  line  2390  and  because 
of  the  ease  in  following  the  program,  it  is  not  necessary  to  discuss 
the  matrix  operations  defining  the  angular  velocity  expressions. 

The  inertia  terms  are  combined  with  the  partial  derivatives  of 
the  angular  terms  in  the  FORTRAN  program.  Also,  the  coefficients 
of  the  acceleration  terms  are  added  to  the  acceleration  coefficients 
of  the  translational  part.  The  remaining  terms  are  added  into  the 
right  hand  sides.  Because  equations  (A-5)  and  (A-6)  were  equal, 
it  was  not  necessary  to  calculate  the  partials  in  (A-^b)  ,  However, 
their  corresponding  inertia  terms  have  different  numerical  values 
and  must  be  combined  appropriately.  Therefore,  equation  (A-6)  is 
set  equal  to  equation  (A-5)  in  subroutine  DER2  so  that  this  can  be 
accomplished.  Thus  the  subscripted  variable  for  the  masses  now 
runs  from  1  to  5. 

Since  differentiation  with  respect  to  the  generalized  velocities 
is  required,  the  definition  of  the  subscripted  variable  changes  for 
the  angular  quantities.  These  are  as  follows  for  the  variable 
PW(I,J,K,L).  I  =  1,2, 3, 4, 5  and  defines  the  mass  that  is  being 
delt  with.  J  =  1,2,..., 12  and  defines  which  partial  derivative  has 
been  taken  with  respect  to  the  generalized  coordinate.  K  = 
1,2,..., 12  and  defines  which  partial  derivative  has  been  taken  with 
respect  to  the  generalized  velocity,  L  =  1,2,3  for  the  x,y  or  z 
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direction.  When  J  =  12  and/or  K  =  12,  this  indicates  no  partial  has 

t  h  th 

been  taken  with  respect  to  the  J  generalized  coordinate  or  the  K 

generalized  velocity  or  both.  An  examination  of  equation  (Br*14)  shows 
both  the  A  and  B  vectors  are  required  with  no  derivatives.  The  deriv-' 
atives  for  equation  (B-IS)  already  exist  from  those  taken  in  equation 
(B'-IA) .  Since  there  are  not  any  similar  terms  in  these  two  equations, 
the  negative  of  (B-15)  is  required  as  was  not  the  case  in  the  trans¬ 
lational  analysis. 

The  partial  derivatives  of  the  subscripted  variable  OMEGA  (I,L), 
which  begins  on  line  2650,  are  defined  to  be  PW(I,J,K,L). 

A  discussion  of  the  derivatives  for  the  potential  energy  will 
now  take  place.  The  potential  energy  function,  PE,  is  defined  to  be 
PE  =  £  U^,  see  equation  (A-15) ,  As  required  by  equation  (2.1),  the 
partial  derivative  of  each  of  the  six  components  of  PE  with  respect 
to  the  generalized  coordinates  must  be  taken.  The  9U^/8q,  did  not 
actually  babe  to  be  performed  as  all  of  the  required  derivatives  have 
already  been  calculated  from  the  kinetic  energy.  This  is  discussed 
in  section  2.5 

The  derivatives  of  and  can  be  followed  in  the  FORMAC 

program  starting  at  statement  number  3710.  In  statement  number  3860, 
LAM25  (1,KKK)  was  LAM25  (3,KKK)  when  aU^/Sq.  was  taken.  The  3  was 
changed  to  a  1  for  the  partial  differentiation  of  in  the  dissipative 
energy.  Working  with  TSO,  this  change  was  very  simple  as  compared  to 
coding  additional  statements.  This  short  cut  is  rather  unfortunate 
as  the  step  by  step  procedure  in  the  FORMAC  program  is  not  entirely 
in  sequence.  Only  230K  bytes  of  core  was  allowed  and  thus  every 
possible  use  was  made  of  every  previous  operation. 

The  subscripted  variable  PU(I,J,K,L)  for  the  potential  energy 

has  the  following  definitions  for  the  subscripts.  I  =  1,2, 3, 4  and 

refers  to  U.  through  U. .  U  and  U  are  added  to  the  right  hand  sides 
1  4  5b 

in  subroutine  NAME.  J  =  1,2,.,. ,11  and  represents  differentiation 
with  respect  to  the  generalized  coordinate.  K  =  12  since  second 
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partials  are  not  taken.  L  =  1, 2,3,4  for  the  as  this  term  was 
broken  up  into  4  terms  due  to  the  4  ground  springs,  see  definition 
of  variable  DA(I,J,K,L)  in  section  3.1.  and  were  each  performed 

in  two  parts  as  indicated  by  equations  (A~13)  and  (A-14) . 

The  dissipative  energy  function,  DE,  is  defined  as  DE  =  + 

where  and  are  given  in  equations  (A-16, 17 , 18) ,  The  DE 

function  requires  differentiation  with  respect  to  the  generalized 
velocities.  Thus,  for  in  equation  (A~13),  let  EXP  stand  for  the 
expression  involving  the  sums  and  products  of  the  matrices.  Then, 


U  =  i  C.  .  ~  EXP 
1  2  ij  dt 


Consider  now 


d  aEXP  •  +  3EXP  •  .  .  3EXP  • 

i  z  n 


9  r  1  o  r  9EXP  •  .  .  9EXP  •  T 

9q.%q.  9q^‘^l+*-  Sq^^^n 

J  3 


„  9EXP  •  . 

C  .  .  ^  q.  + 

ij  1 


9 EXP  •  9EXP 

9q  3a . 

n  -j 


Therefore,  only  the  8EXP/8q^  is  calculated  and  the  suras  and  products 
of  terms  are  performed  in  subroutine  NAME,  The  utilization  of  the 
calculations  performed  in  the  potential  energy  were  taken  advantage 
of  in  obtaining  the  expression  EXP, 

1  *2 

The  necessary  derivatives  for  the  dissipative  function  ^2  ”  2  ^i  ^^2^ 
may  be  obtained  by  the  following  analysis. 


■  3L 

d  ^ _ 2  ;  .  _ 2  • 

dt  2  9q^  ^1  3q2  ^2 


^n 
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is  assembled 


9U  .  3L  .  3L 

--  =  C.  (L„)—  =  C.  (L  ) 

■i*  1  1  2  3q 

9qk  k 


9L. 


Thus,  only  SL^/Sq^.  needs  to  be  determined, 
in  subroutine  NAME, 

The  same  analysis  is  performed  for  U^,  That  is. 


3U 


3q 


3  .  2 

-  =  6(L.)^  ^ 

3  9q,. 


where 


.  9L3  .  9L3  . 

^ —  q.  +  ^ —  q  + 
3  1  9qo  2 


3L, 


3q  n 
n 


The  required  derivatives  for  the  generalized  forces  are  defined 
in  section  A6  of  Appendix  A.  The  reader  can  easily  follow  through 
the  FORMAC  program  by  ref ering  to  section  A6  and  subroutine  NAME, 
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4.0  CONCLUSIONS  AND  RECOMMENDATIONS 

This  report  exhibits  the  techniques  of  utilizing  FORMAC  for  the 
generation  of  the  equations  of  motion  of  a  complex  weapon  system. 

The  development  of  this  technique  was  a  major  objective  of  this  effort. 
It  allows  the  analyst  to  spend  considerably  more  time  formulating  the 
problem  and  not  be  too  concerned  about  the  enormous  amount  of  math¬ 
ematics  to  be  performed.  It  also  allows  a  completely  nonlinear  model 
to  be  developed  and  thus  it  is  no  longer  necessary  to  linearize  the 
generalized  coordinates  to  achieve  simplification  of  expressions.  The 
procedure  of  linearizing  is  very  time  consuming,  nearly  impossible 
for  large  degree  of  freedom  systems,  and  the  risk  of  generating 
numerous  errors  is  quite  high.  Utilizing  FORMAC,  the  differential 
expressions  are  punched  on  cards  in  FORTRAN  format  which  eliminates 
any  potential  key  punching  errors. 

A  computer  model  has  been  developed  for  the  simulation  and  is 
operational.  The  model  output  appears  reasonable  based  on  current 
data  available  and  a  knowledge  of  how  the  system  performs  under 
dynamic  conditions  of  firing.  Due  to  the  impending  closure  of  Rodman 
Laboratory  it  was  necessary  that  this  report  be  written  in  a  limited 
time  frame  since  both  authors  were  leaving  the  laboratory.  The  work 
presented  in  this  report  was  initiated  during  the  month  of  October, 
1976  and  was  terminated  in  January,  1977.  Therefore,  some  refinements 
and  corrections  to  the  model  which  would  have  been  made  if  time  per- 
mitted  are  discussed  below. 

A  better  understanding  of  the  friction  brakes  in  both  the 
elevating  and  traversing  mechanisms  and  an  improved  logic  criteria 
is  desired.  The  value  may  be  a  "break  away"  torque  and  under 

dynamic  conditions  is  inaccurate.  A  more  desirable  logic  might  be 

(IY)|  --  theny  =  Y=  0soAy  =  0 
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if 


_d^ 

dt 


(Signium  y) 


^  ^EMAX  ^E  ~  ^EMAX 


where 


'  -  "tx  'bx  "-"r  <"  ■"  "e('1  + 


A  more  accurate  value  of  3„  be  obtained  when  test  data  becomes 

E 

available. 

At  the  time  of  this  analysis,  only  the  total  resisting  force, 

R(t)  +  C(t),  was  available  and  was  defined  to  be  R(t),  These  forces, 
R(t)  and  C(t),  should  be  Individually  determined  as  functions  of  the 
generalized  coordinates  to  permit  parametric  variation  studies  on 
orifice  area  effects  and  to  justify  optimization  design  studies  as 
meaningful.  The  computer  program  has  been  written  to  accommodate 
these  changes. 

Determination  of  the  hydraulic  force  coefficient  (3)  resulted 
in  such  a  large  force  (due  to  the  exceedingly  small  orifice  area) 
that  it  exceeded  the  force  expected  from  deformation  of  the  cylinder. 
Therefore,  the  spade  cylinder  was  ignored  in  the  computer  runs  by 
setting  3=0.  However,  if  modification  to  the  spade  cylinder  are 
made  to  allow  the  cylinder  to  act  as  a  shock  absorber,  the  design  of 
the  cylinder  may  be  evaluated  by  this  model, 

Initial  conditions  may  be  determined  in  two  different  ways; 

3P  E 

(1)  set  ^  “  0  and  solve  the  resulting  system  of  algebraic 

equations  for  initial  values  of  the  generalized  coordinates  or  (2) 
assign  a  dashpot  with  each  spring  and  without  applying  B(t),  let  the 
computer  solve  the  system  of  equations  until  equilibrium  is  obtained. 
The  resulting  values  of  the  generalized  coordinates  are  then  used  as 
initial  conditions.  Due  to  time  limitations,  a  simple  static  analysis 
was  used  to  define  Initial  conditions  for  the  sample  output  of 
Appendix  F. 

A  secondary  output  from  the  program  is  the  dynamic  loading  on 
the  various  components.  Spring  loads  are  directly  determined  from 
the  program.  However,  to  determine  interaction  loads  on  various 
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components,  free  body  diagrams  will  have  to  be  developed  and  equations 
of  motion  written  for  the  components.  Since  accelerations  of  the 
components  are  known  from  the  program,  it  only  remains  to  solve  for 
the  forces  which  should  then  be  compared  with  the  results  from 
experimental  test  data.  Finally,  the  model  should  be  validated  against 
data  obtained  from  field  tests  which  are  to  take  place  during  FY  77, 
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APPENDIX  A 


Derivation  of  Energy  Expressions  and  Generalized 
Forces  for  MllO  Self-Propelled  8”  Howitzer 


A-1 


APPENDIX  A 


Derivation  of  Energy  Expressions  and  Generalized 
Forces  for  MllO  Self-Propelled  8"  Howitzer 

Al,  Coordinate  Systems 


0^  -  ABC 


0^  -  A'B'C’ 


0^  -  A"B”C" 


0^  -  XYZ 


0,  -  X*Y’Z» 
4 


an  Inertial  coordinate  system  having  its  origin 
at  the  mass  center  of  the  weapons  system  combi¬ 
nation  in  the  firing  position.  The  0^  -  AB  plane 
is  horizontal  and  parallel  to  the  earth’s  surface. 

The  0^  -  BC  plane  is  vertical. 

a  coordinate  system  having  its  origin  fixed  at  the 
center  of  mass  of  the  weapons  system  combination. 

This  coordinate  system  is  at  all  times  parallel  to 
0^  -  ABC.  The  coordinates  of  0^  in  0^  -  ABC  are 
(x,  0,  0). 

a  coordinate  system  reflecting  the  roll  (0)  and 
yaw  (ip)  of  the  weapons  system  combination.  The 
coordinates  of  the  mass  center  of  the  vehicle  (0^) 
are  located  in  this  coordinate  system  by  (a^,  a^,  ^3)* 

a  coordinate  system  having  its  origin  at  the  mass 
center  of  the  vehicle  prior  to  the  y  and  z  transla¬ 
tions  of  the  vehicle  and  is  at  all  times  parallel 
to  0^  -  A”B”C”. 

a  coordinate  system  having  its  origin  fixed  at  the 
mass  center  of  the  vehicle.  This  coordinate  system 
is  at  all  times  parallel  to  0^  -  XYZ.  The  coordi¬ 
nates  of  0^  in  0^  -  XYZ  are  (0,  y,  z). 


0  -  X”Y**Z'*  a  coordinate  system  fixed  in  the  vehicle  and  re- 

4 

fleeting  the  pitch  (^)  of  the  vehicle.  The  end 
points  of  the  braces  and  spade  cylinders  are  located 
in  this  coordinate  system  by  (A^,  A2,  ^3) >  ^2* 

A^)  and  (B^,  B^) ,  (B^^,  B^,  B^)  respectively. 


a  coordinate  system  having  its  origin  at  the  center 
of  the  traverse  bearings;  0^  -  Z**’  being  the  center 
of  traverse.  This  coordinate  system  is  at  all 
times  parallel  to  0^  -  X*’Y"Z*’.  The  coordinates 
of  0^  are  e^,  62  >  in  the  0^  -  X”Y’*Z^*  coordinate 
system. 

a  coordinate  system  that  reflects  the  traverse  an¬ 
gle,  T,  The  coordinates  of  the  mass  center  of  the 
traversing  but  non-elevating  parts  are  f^,  f2,  f^ 
in  this  coordinate  system. 


0.  -  EHZ 
6 


0^  -  H'H*Z’ 


a  coordinate  system  having  its  origin  fixed  at  the 
midpoint  between  the  trunnions.  This  coordinate 
system  is  at  all  times  parallel  to  0^  -  X^^Y^^Z^^. 

a  coordinate  system  that  reflects  the  pitch  of  the 
gun  (y)  relative  to  the  understructure.  Note  that 
the  initial  value  of  y  is  the  angle  of  elevation. 
The  coordinates  of  the  mass  center  of  the  recoiling 
parts  are  in  this  coordinate  system.  The 

coordinates  of  the  mass  center  of  the  elevated 
but  non- recoiling  parts  are 
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0^  uvw 


a  coordinate  system  having  its  origin  at  the  center 
of  pressure  of  the  spades.  The  coordinates  of  0^ 
in  0^  -  XYZ  are  (g^,  ,  g^) . 


0^  -  U’V’W’ 


a  coordinate  system  fixed  in  the  spade  and  reflecting 
the  pitch  (v)  of  the  spade.  The  coordinates  of 
the  mass  center  of  the  spade  are  (h^,  h2 ,  h^)  in 
this  coordinate  system.  The  coordinates  of  the 
end  point  of  the  spade  cylinders  and  braces  attached 
to  the  spade  are  (+  ^3)  (it  ^2^  ^3) 

respectively . 


A-4 


A3.  Kinetic  Energy 


Define 


Mass  of  vehicle/spade  combination 

Mass  of  vehicle 

Mass  of  spade 

Mass  of  traversing  but  non-elevating  parts 

Mass  of  elevating  but  non-recoiling  parts 
Mass  of  recoiling  parts 


Define  the  coordinates  of  the  mass  center  of  M 


M  :  X,  0,  0  in  0-  -  ABC 

vs  1 

M  :  a, ,  a^,  a^  in  0.  -  A"B"C’* 

V  12  3  2 

M  :  h, ,  h.,  h.  in  0^  - 

S  i^j  / 

Mt  :  i 2'  ^3 

\  n^,  In  Og  -  E’H'Z’ 

M  :  5,  n,  C,  in  0,  -  5’H’Z’ 

r  o 

Define  the  vectors  in  the  0^  -  ABC  coordinate  system  from  0^  to 
the  mass  center  of  M^  as 


and  the  vectors  to  an  element  in  M^  as 


Define  ?  as  the  vector  from  the  mass  center  of  to  the  element  of 

e.g. 

-y  -y 

Q  +  Pp  -  q 


The  angular  velocities  are  determined  as  follows 


px 

-1 

'0* 

• 

0) 

e 

py 

b 

• 

(JL) 

pz^ 

A'*B"C" 

Mi* 

and 


"U3 

”0) 

T 

qx 

-1 

px 

0) 

to 

•f 

0 

qy 

5 

py 

0) 

0 

qz 

X”Y"Z” 

- 

p* 

Also 


0) 


rx 


0) 


ry 


0) 


rz 


-1 

11 


u’v’w’ 


"  (0 

V  1 

px 

0) 

py 

0 

(0 

pz 

0 

1.  ^ 

(A-1) 


(A-2) 


(A-3) 


A-IO 


Now 


A-11 


so 


Now 


Q  +  cOq  X  Pq 


M 


q  dM^  -  is  [  Q  •  Q  +  2Q  •  Up  X  *  PQ  >  •  ^  *  PQ^  ‘^v 

M 


or 


q  dM^  =  %  Q  Q  +  Q 


M 


‘^Q  |pQ  ^  [(“q  Pq) 


M 


M 


(Up  X  Pp)  dM^ 


but 


1 


Pq  dM^  •  0 


M 


and  neglecting  cross  products  of  inertia 

I5  |(^p  ^  Pp)  •  ^  Pq)  d^v  “  “qX  ■"  V  “QY  ^QZ  "qZ^ 

M 


So 


KE  =  %M  Q^  +  %M  R^  +  %M  S^  +  i^M  T^  +  i$M 
V  s  e  r 


+  h 

<Iqx 

^QX 

+ 

^QY 

^QY 

+ 

^QZ 

“QZ> 

+  J5 

<'rx 

+ 

3 

^RY 

•"ry 

^RZ 

“rz 

+  % 

“sx 

+ 

^SY 

“SY 

+ 

^SZ 

“sz^ 

+  h 

^^TX 

“^x-^ 

Ity 

+ 

^TZ 

2 

TZ 

+  %  0)^  +  0)2^) 


t  b 

where  is  the  absolute  angular  velocity  of  the  i  body  around  the  j  axis 
of  that  coordinate  system  fixed  in  the  body. 
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Now 


(A-7) 


(A-8) 


Xi  + 


^4  ■"  S 


Xe  + 


(A-9) 


Xi  +X2 


4*  X 


(A- 10) 


+  X, 


(x^  +  X4  +  Xj 


+  X. 


(A-11) 


A4.  Potential  Energy 
Define 

Xo-  [0  0  1] 

Then  the  potential  energy  of  the  component  weights  is  given  by 
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"l  ■  fv  *0 


+  ^2  ^z,)!  +  M  A  A  A  +  1  j.  1 

J  a  0  I  1  2  [^3  ^  ^10  ^11 


+  A  +  A- 

3  4  5 


^6  +  S 


+  M  A- 
e  0 

Ai  +  A2 

A.  +  A,  +  A^ 

345 

■'i 

*• 

u 

^1 

•0  j 

• 

i 

4 

(A-12) 


Ai  +  A2 


A3  +  A^  +  A5 


Ag  +  A^ 


A8+  Ag 


Assume  vertical  ground  springs  at  the  front  and  rear  roller  wheels.  In 
0^  -  X*’Y”Z**,  define  the  coordinates  of  the  point  of  contact  between  ground 
and  roller  wheels  as  j  ’  ^ij  ^  ^  =  1>2  where  i  =  1  is 

right,  i  «  2  is  left,  j  =  1  is  front  and  j  -  2  is  rear. 


The  coordinates  of  the  roller  wheels  ground  contact  in  0^  -  ABC  are 


ri  ' 

•1 

ij 

A_  +  A,  +  A_ 

m. , 

3  4  5 

ij 

L  ij 

and  the  extension/contraction  of  the  springs  is 


- 

P 

1. , 

- 

ij 

ij 

A^ 

A^  4”  Art 

Art  4"  A,  4*  A- 

m. . 

A^  + 

m. . 

0 

1  2 

3  4  5 

ij 

3 

ij 

n 

n 

ij 

ij 

and  the  energy  stored  in  the  springs  is 
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if  6  - .  <  0  then  K . .  *  0 
ij 


The  coordinates  of  the  end  points  of  the  braces  prior  to  motion  are 
(in  0^  -  U'V'W) 


Thus  the  original  length  of  the  brace,  is 


Thus  the  length  of  the  brace,  L2,  is 
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Ai* 


-  A 


"2 


°3J 


The  energy  stored  in  the  braces  Is 


U.  =  J5  K 


2  (L^  -  L2). 


i  =  1,2 


(A-13) 


The  coordinates  of  the  attachment  of  the  two  ground  springs  (always  per¬ 
pendicular  to  the  spade)  in  ©2  -  ABC  are  (+  02,  a^) .  In  the  0^  -  U'V'W* 

system  the  coordinates  are 


-1  !  ,-i 

11 


-  “1 


"  ^3  ■  ^10 


Prior  to  motion,  the  coordinates  of  the  spring  attachments  in  0^ 
are 


-  U’V'W’ 


r 

—  ““ 

ai  + 

“2 

- 

^2  ^2 

“3 

^3  % 

_ 

1—  — 

The  displacement  of  the  springs,  6^,  is  the  "y"  component  in 


-1 

11 


-1 


•± 

\  N 

/ 

±“l 

■  +  g^- 

“2 

-  h  -  S  -  ho  , 

!■( 

“2 

- 

^2  H 

-  “3 

/  J 

V 

“3 

^3  ^  H 

and  so 


i  =  1,2 


(A-14) 
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The  energy  associated  with  the  tipping  j^arts  is  defined  as 


E  L 


f^2  [~  ‘''3  ”  "  *^2^  Y  +  (0^  -  d^)  Cos  y] 


+  (O2  -  d^)  Cos  Y  +  (O3  -  d^)  Sin  y 


L={  N2  -Y-M  O2  -X3 


2  ^  1/2 


where 


=  B(t)  •  C„until  B(t)  •  C  >  then  Bg  = 


when  Y  <  0>  Pg  *  0 


Also,  the  energy  associated  with  the  traversing  parts  is  defined  as 


where 


U6  =  e,  •  T 


=  B(t)  (-Sg)  until  B(t)  Ug  I  > 


then  B^  -  when  t  <  0,  =  0 


The  potential  energy  function,  PE,  is  defined  as 


PE  =  lU. 
1  ^ 


(A-15) 


A5.  Dissipative  Energy 

The  dissipative  function  associated  with  the  roller  wheel  springs  is 
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The  dissipative  function  associated  with  the  brace  is 


or  since  -  0 

The  coordinates  of  the  end  points  of  the  spade  cylinders  during  motion 
(in  0^  -  U’V’W’)  are 


Thus 

=  1/3  e  (Lg)^  (A-18) 

The  dissipative  energy  function,  DE,  is  defined  as 
3_ 

D.E.  =  lU 
1 


-16) 


A-18 


Note  that  there  is  no  dissipative  function  associated  with  the  ya\v  and 
pitch  springs  of  the  traversing  and  elevating  parts,  respectively.  Also 


where  o  =  specific  weight  of  oil,  A  -  area  of  piston,  and  A^  effective 
are  of  orifice. 


A6.  Generalized  Forces 


The  breech  force  acts  at  in  the  0^  -  coordinate 

system  and  in  the  ~n  direction.  The  recoil  force,  R(t) ,  acts  at 
(^A,  in  the  0^  -  H’H’Z*  system  and  in  the  n  direction  but  is 

an  external  force  in  the  n  equation  only.  The  components  of  the  breach 
force  in  0^  -  ABC  are 


B(t)^ 

=  X2  Xg 

J 


In  -  ABC  the  coordinates  of  the  points  of  application  are  Ag,  Bg,  Cg 


The  generalized  force  is 


,  where  =  n  +  e 

D 


Similarly 


and 


and 


Also, 


where 


and 


RCt)/ 

r 

Q"n  = 

R(t)„ 

0 

n 

"rJ 

,  where  \  ^ 


M‘)a 

ll(t)B 

-  »2  •  Xj  •  X7  • 

1 - 

0 

1 _ 

Lo  J 

»R 


Xi.X2(X3-.X,.X3(X6-HX7(X^-.X^ 


I  recuperator  force,  C(t),  acts  at  C^-  Thus 

c(t). 


n 


c(t) 


RJ 


where  n  -  r\  +  i 
c 


Lc(t)cJ  [c. 

C(t)A 

"  0 

C(t)g 

=  X2  •  X5  •  X?  •  Xg 

C(t) 

Lc(t)J 

0 

B 


Xi  +  X2  A3  4-  /  Xe  +  x?  Xg  +  X^ 
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Define 

X'.  [l  0  o] 

then  the  displacement  of  the  horizontal  ground  springs  is 


and  the  horizontal  force  is  -  K' V The  generalized  force 
such  that 


if  I  -  K-,,  6%  J  <  y  K^.  6^. 


ij  ij 


ij  ij 


-y  K,,  (signum  if  j  -  6'^^)  > 


ij  ij 


Note  when  6^^  £  0,  “  ^'ij  ~  ^ 
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APPENDIX  B 


Mathematics  of  Solution  Technique 

Since  FORMAC  takes  only  partial  derivatives  the  following  analysis 
was  required  to  determine  which  derivatives  are  actually  needed  and 
once  they  are  obtained  in  what  manner  are  they  combined  to  produce  the 
equations  of  motion. 

To  begin  the  analysis  the  expression  for  the  Lagrangian  which 
yields  the  equations  of  motion  is 


JL  9(KE)  _  a(KE)  3(DE)  a (PE)  ^  p 

dt  aq^  aq^  Bq^  Sq^  j 


(B-1) 


where 


j  ”  2,  •••>  k 

KE  =  Total  kinetic  energy 

DE  =  Total  dissipative  energy 

PE  =  Total  potential  energy 

F.  =  Generalized  Force 
J 

q .  =  Generalized  coordinate 

q.  =  Generalized  velocity 

t  =  Independent  variable,  time 

k  =  Number  of  generalized  coordinates 


and  write 


Consider  at  first  only  the  kinetic  energy  terms  of  equation  (B-1) 


d  a (KE)  a(KE)  ^  d  a (t+w)  a (t+w) 


dt  3q. 


dt  8q. 


(B-2) 


where  T  =  T  is  the  translational  part  of  the  kinetic 


i=l  energy 
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w  =  EW.,  W  is  the  rotational  part  of  the  kinetic  energy 
i=l^ 

n  =  Number  of  masses 

•  • 

=  1/2 

Wi  =  I  p^dM^+1/2  [  (m^xp  .).(w^Xp^)dM^ 

1  1 
,  th 

M  =  1-  mass 
i 

^  .  th 

Q,  =  Position  vector  of  the  i-  mass 

Q  =  Vector  from  the  mass  center  of  M  to  an  element  of  M 
=  Angular  velocity  of  mass 

Now  consider  only  the  and  dropping  the  subscript  on  the  T  for 
convenience  as  though  only  one  mass  is  being  studied  at  this  time*  An 
examination  of  the  expression  (B-3)  will  proceed. 


_d_  3Q-Q 
dt  9q 


(B-3) 


dq  q  q  q 


where 


Since  Q  is  a  function  of  the  generalized  coordinates  only,  i.e., 


Q  -  Q(q2’^2’  ’  *  ’ 


dt  ^  ^  3q2^2 


+  A  =  ^  q 


(B-4) 


B-3 


where  the  double  subscript  on  i  refers  to  the  Einstein  summation 
notation  and  runs  to  the  value  k.  From  (B~4) 


=  Q  .q^  +  Q  ,  +  ^  -q,  +  ^  q^.  +  ...  +  Q  +  Qq 

q,q  1  qi  iq  q2‘i  ^  q2  2q  qj^q  «•  \ 


Since  Q  is  a  function  of  only  the  coordinates 


■5r*9i2  ■55'"^  •••  H  »1 


Q-  -  ? 
q  q 


(B-5) 


Therefore, 


i|2 . 24.4  -  24,4 


4^  r 

in.n  I  A  I 


dt  3q  J  dt  ^q 


For  one  mass,  say  T  =  1/2MQ*Q,  then 


—  Ill  =  A  r ^1/2 

dt  9q  dt  9q 


MQ  -Q 

q 


Examination  of  the  term 


yields 


+  [Q^.q 


*  '\'^t52^2  +  •••  + 


di-  %-Ql  -  [q„^-3  -e  q,-5  1^1  +  •  •  •  +  (5  '5  + 


•  •  •  • 

+  [3  .  +  3-3.  ]q,  +  ...  +  [3_.  ^  -Q.  ]q 


•q  ’1 


q  "q,,^"k 


Since  Q  is  a  function  of  only  the  coordinates,  terms  like  Q  .  go 

qq^ 

to  zero  and  using  the  results  of  equation  (B-5)  the  above  expression 
reduces  to 


^  =  1%/^  +  +  •••  +  +  %\  ]q 


q 


'*"Q*Q  Qi  +  •••  +  Q*  Q  q 


q  X-'k 


(B-6) 


B-5 


Since 


iQ..^  +  +iS_5 

3q2^2  +  •••  +  aq^\ 


a^Q  .  ,  32q  .  ^  .  92$  . 

3q^3q*^l  +  3q23q  *12  3qj^3q 


\q’l  <“-’> 


And  using  the  results  of  equations  B-4,  7,  equation  B-6  becomes 

q..  + 


_d_ 

dt 


- 1 

1 

_ 1 

JJ. 

^  -1 

Qno  *^n  5  -Q  4,. 

q-  +  ,  .  •  + 

L’  J 

qqj^  q^  1  q  i 

qqk  qi  1 


^  q . 


•••  *1 


q  'q^’k 


qi  Q  ‘  Q  q,-  +  • 
1  qq^^  q^  1 


.  +  q  5  •  Q  q . 

k^qqj^  qi  ^ 


->  -v 


+  Q  -Q  q.qq  +  .. .  +  Q  *0  q.q, 

q  qfqi  1  1  q  qiqk  ^  ^ 


+  Q  *Q  q,  +  .. .  +  Q  ‘Q  8, 
^q  qq  ^q  V  ^k 


qi.5  '*‘•••+4,  ^  *Q  4- 

1  qqi  qqj  qi  ^ 


+  Q. 


•Fq  q. 
LVi  " 


q-  +  .  . .  +  Q 


A 'a] 


"l  ■"  •••  « 


Ok’*' 


q.Q  ‘Q  q.  +  Q*$  q.q.  +  ^*Q  ?1. 

j  qqj  qi  1  q  q^q^  i  j  \ 
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Thus,  for  one  mass  the  following  result  is  obtained 


.t,  _  •  A  q.  +  Q  'Q  q.q-i 

dt  3q  dt  \  ^qq .  \  ^1  \  ^  ^ 


(B-8) 


When  more  than  one  mass  is  involved  a  summation  over  all  masses 
must  be  performed.  Equation  (B-8)  would  take  the  form 


,  31  T. 
_d _ 1 

dt  3d 


(B-9) 


Before  proceeding  with  the  rotational  part  of  equation  (B-2) ,  an 
examination  of  the  translational  part  of  the  second  term  of  that  equation  is 
given  below.  For  ease  of  notation  consider  only  one  mass.  Then, 


3q  9q 


1/2MQ-QJ  ;  q  =  J^q^,  ^2*  ] 


^4  M  4  4 


1/2MQ  -  Q  =  “  Q^-Q  +  Q  -Qq  =2  ^ 


M  r.4  4 


aT  ^  .  T  . 

—  =  MQ  •Q=MQ  q.*Q  q. 
3q  q  ^^q  1  qj  J 


(B-IO) 


from  B-4  and  B-7.  The  results  of  equation  (B-10)  also  appear  in 
equation  (B-8);  and  since  (B-10)  is  preceded  by  a  minus  sign,  these 
terms  cancel  and  equation  (B-2)  becomes 


(B-11) 


A  ^  ^ 

dt  8q  dq 


=  EM 


Q  •Q  4.4.“^Q  ‘Q 

q  ^^q^  ^q  ^q^ 


n 

It  was  seen  earlier  that  W  =  E  W.  ,  where 

i=l  ^ 


W.  = 
1 


‘  ^iJm  ^  a  |m 

i  '  1 


.)  dM 


=  1/2  I  (i)  0)2(1)  +  1/2  I  (1)  0)2(1)  +  1/2  I_(i)  0)2(1)  (B-12) 

A  j  j  y  Z  Z  Z 


fl  (i)  0)  (i)  0)  (i)  +  I  (i)  0)  (i)  o)  (i)  +  I  (i)  (jj  (i)  o)  (i) 
L  xy  X  '  y  '  xz  ^  X  z  yz  y  z^ 


I 

ince  ~  Equation  (B-12)  can  be  written  in  matrix  form 


as 


M 


W.  =  A(i)  B(i) 

Where  A(i)  and  B(i)  are  defined  as 


A(i)  = 


1/2  I  (i)  0)  (i)  -  I  (i)  0)  (i) 
XX  X  xy  y 


1/2  I  (i)  00  (i)  -  I  (i)  0)  (i) 
yy  y  yz  z 


1/2  I^^(l)  0)^(1)  -  I,, (i)  0)^(1) 


IT 


;  B(l)  = 


0)  (1) 
X 


uiy(i) 


u)^(i) 


The  subscripts  are  dropped  again  as  a  matter  of  rotational  convenience. 


9W 

^  =  W.  =  (AB)  .  =  A.B  +  AB. 
dq  q  q  q  q 


Since  W  is  a  function  of  the  generalized  coordinates  and  generalized 
velocities  i.e.,  W  —  Wiq^^,  ^2*  ****  ^1*  ^2*  ***^  then 
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^  ‘"4’-  %,  “’i  * \\ ’1  ^  \\\ 


+  ^  \‘k\* 


(B-13) 


^  Wl'^  '”  \^q^k  ■" 


However,  terms  like  A.,  and  B..  go  to  zero  because  each  of  the 

qq^  qqi 

vectors  is  linear  with  respect  to  the  generalized  velocities.  Equation 
(B-13)  reduces  to 


fw.  1  =  q^A.  B  +  A.B  q,  +  q.A  B.  +  AB.  ^ 

dt  L  q  J  ^1  qq^  q  qj  j  i  q^  q  qq^  ^ 


+  A.B  q  +  q.A.  B. 
q  qH  q^  q 


(B-14) 


The  remainder  of  the  rotational  terms  are 


—  =  A  B  +  AB 

9q  q  q 


(B-15) 


The  above  analysis  describes  which  partial  derivatives  must  be 
taken  and  how  they  are  combined  for  the  kinetic  energy.  The  derivatives 
of  the  other  energy  expressions  are  either  obtained  from  the  kinetic  en¬ 
ergy  or  are  completed  separately.  For  the  most  part,  those  differential 
expressions  which  cannot  be  obtained  from  the  kinetic  energy  are  easily 
calculated  by  the  use  of  FORMAC,  In  order  for  the  reader  to  see  more 
clearly  the  use  of  equations  (B-11) ,  (B-14) ,  and  (B“15) ,  a  discussion  is 


B-9 


given  in  section  2.0  with  FORTRAN  to  see  exactly  how  terms  are  combined. 
Also,  in  section  3.3  this  is  discussed  further  with  actual  application 
to  an  artillery  problem. 
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APPENDIX  C 


FORMAC  Program 


AIM'ENDIX  C 


FORMAC  Trograni 

This  appendix  contains  a  listing  of  the  FORMAC  program  which  was 
used  to  generate  the  differential  expressions  to  solve  Lagrange’s  Equa¬ 
tion.  This  program  was  developed  on  the  IBM  360/65  operating  system 
Time  Sharing  Option  (TSO) .  A  maximum  of  230K  bytes  of  core  storage  was 
available  and  thus  segments  of  the  program  were  commented  out  when  not 
utilized.  Because  of  the  ease  in  which  TSO  converses  with  the  user, 
the  values  of  some  variables  were  changed  to  accommodate  a  quick  form¬ 
ulation  for  the  next  energy  expression  to  be  evaluated.  At  all  times 
maximum  use  was  made  of  all  previous  coding  and  so  some  expressions  may 
appear,  at  first  glance, out  of  place. 
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ee32e  LET(CHANll-CHAIN(ZZ(54)««a  *  ZZ(4S)ttB  *  ZZ(46 )tt2, ZZ(S6  ) . 

00330  ZZ(54)XX2  *  ZZ(5S)««2  ZZ(49)»«2.ZZ(57} ) ). 

00340  LET(CHAN12«CHAIN(-ZZ(S0)XX  S  *  ZZ(S7)tt  5.ZZ(58))). 

00350  LET(CHAN13-CHAIN(XKK8»ZZ(58)»Z2(5S  )«*(»  §  ).ZZ(59))). 
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APPENDIX  D 


The  Square-Root  Method 


This  appendix  presents  the  method  that  was  used  to  decouple  the 

•  • 

acceleration  terms,  i.e.,  solve  the  matrix  equation  AX  =  B.  The  method 
can  also  be  used  to  solve  any  system  of  linear  equations  of  the  form 
AX  =  B  if  the  A  matrix  is  symmetric,  which  is  the  case  when  solving 
Lagrange's  equations  of  motion. 

The  algorithm  used  is  called  the  square-root  method.  In  case  the 
coefficient  matrix  of  a  system  is  symmetric,  finding  the  solution  is 
made  more  expedient  by  taking  advantage  of  the  symmetry.  Using  the 
equations  presented  in  Reference  1  produces  pure  imaginary  numbers  in 
the  computational  scheme  and  so  a  modification  to  that  algorithm  is 
given  here  to  alleviate  the  problem  of  imaginary  numbers. 

A  FORMAC  program  was  developed  using  the  modified  square-toot 
method  to  solve  a  system  of  equations  either  numerically  or  symbolically. 
The  program  is  listed  here  along  with  results,  both  numeric  and  sym¬ 
bolic.  Cards  can  be  automatically  punched  in  FORTRAN  format. 

FORMAC  can  be  used  to  obtain  S3nnbolic  solutions  in  problem  areas 
which  heretofore  could  only  be  approached  numerically.  This  is  well 
exemplified  in  this  appendix. 

The  solution  of  a  system  using  the  square-root  method  reduces  to 
the  solution  of  two  triangular  systems.  For  the  equation 
AX  -  F 

the  algorithm  according  to  Reference  1  is 
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The  final  solution  is  found  by  the  formulas 


1  k.  -  Z  s.  X 

k  1  . ,  -  im  m 

n  m=i+l  .  ^ 

X  =  —  ,  X,  =  -  ;  1  <  n 

n  s  1  s.  . 

nn  11 


In  case  the  elements  of  the  matrix  are  such  that  radicands  of  the  ex¬ 
pression  s^^  are  negative,  pure  imaginary  numbers  appear  in  the  r  w 
for  which  s^^  <0.  To  alleviate  this  problem,  the  following  defini¬ 
tions  are  made 


^11 


s  .  =  a  . 
13  ij 


Cl  =  1/Sii 


c .  =  1/s  . . :  i  >  1 

1  11 


s..=a..-Z  s.c,i>l 
11  11  1  mi  m 

m=l 


s..=a..-Z  s.s.c,i>i 
13  13  mi  mj  m’ 


i  >  3 


ki  =  fi 


Z  s  .  k  c 
m=l  mi  m  m 


X  =  k  c 
n  n  n 


X.  =  c.k. 
1  11 


s .  X  c .  :  i  <  n 
im  m  1 


D-3 


n 

Z 

m==i+l 


With  this  method,  approximately  n  /2  elements  of  the  S  matrix  and  2n 
components  of  the  vectors  K  and  X  are  recorded.  The  square-root  method 
is  widely  employed  where  the  solution  of  symmetric  systems  is  called 
for  and  it  is  recommended  as  one  of  the  most  efficient  methods. 

A  listing  of  the  FORMAC  program  is  given  below.  To  solve  the  system 
of  equations  AX  =  F,  only  the  value  N  representing  the  size  of  the  square 
matrix  A  and  the  variable  SOLVE  must  be  changed  to  produce  either  a 
numeric  or  symbolic  solution. 

For  N  =  11  and  SOLVE  =  NO,  the  output  is  contained  in  the  FORTRAN 
listing  under  the  subroutine  named  SOLVE  in  Appendix  F. 
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00750  CHAREX(LINE»XK( I ) ); 

00760  CALL  PUNCHl ( LINE  ) ) 

007-0  IF  IDENT(SOLOE;yES)  THEN  GO  TO  LOOPS* 

00780  ATONIZECXKd  ) ), 
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For  N  =  3  and  SOLVE  =  NO  the  following  calculations  as  computed  by 
FORMAC  in  symbolic  form  solve  the  system  of  equations  with  symmetric  matrix 

A. 


^3x3  Sxl  ^3x1 


POTENT  CARDS  DATA 


READV 


S(l. 1 ) 
s(i.a) 
S(l,3) 
C(l)  • 
S(2.3) 
CO)  • 
C(3)  • 
XKd  ) 
XK(2) 
XIC(3) 
X(3)  ■ 
X(2)  • 
X(l)  ■ 


•  Add) 

-  Ad. 2) 

-  Ad. 3) 

Sd.l)»»(-1  ) 

•  A(2.3)-Sd.2)»S(1.3)tCd  ) 
(A(2.2)-Sd.2)*»2*C(l )  )»*(-!  ) 
(A(3,3)-S(2.3)t*2tC(2)-Sd.3)*»2*Cd  )  )»*(-! 

-  Fd) 

-  F(2)-X)Cd  )»Sd  .2)tC(l ) 

-  F(3)-XICd)»Sd.3)»C(l)-XK(2)*S(2.3)*C(2) 
XIC(3)»C(3) 

XK(2)»C(2)-X(3)»S(2.3)*C(2) 

XICd)<C(l)-X(3)*S(l.3)tCd)-X(2)*Sd.2)tCd  ) 


) 


For  N  =  15  and  SOLVE  =  YES  (with  the  appropriate  A  matrix  now 
defined) ,  the  following  numeric  solution  is  given.  The  solution  of 
course  is  the  vector  x(l)  through  x(15). 
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APPENDIX  E 


Removal  of  Imbedded  Terms 

This  appendix  displays  the  reduction  in  the  number  of  arithmetic  op¬ 
erations  that  were  accomplished  by  using  FORMAC.  CHAIN  and  REPLACE  op¬ 
erations  enabled  character  strings  involving  algebraic  expressions  to  be 
replaced  by  new  variable  names  and  thus  eliminate  millions  of  arithme¬ 
tic  operations  during  the  execution  of  the  program. 

The  first  five  equations  (53  lines  of  FORTRAN  formated  output)  rep¬ 
resent  the  partial  derivatives  of  (  part  2  of  two  parts  )  in  the  po¬ 
tential  energy  with  only  the  sines  and  cosines  replaced,  i.e.,  S6  = 

SIN  (Q(6)),  C8  =  C0S(Q(8)),  etc.  The  last  five  equations  (7  lines)  are 
the  result  of  removing  imbedded  terms. 

A  similar  type  of  reduction  in  the  number  of  operations  was  performed 
on  all  expressions  before  the  FORMAC  output  was  punched  on  cards. 
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APPENDIX  F 


FORTRAN  Program  and  Sample  Data 

This  appendix  contains  a  listing  of  the  FORTRAN  program  and  a 
sample  of  the  output  data.  The  correspondence  between  the  column 
headings  of  the  output  and  the  Q(i),  QD(i),  QDD(i)  along  with  the 
units  can  be  obtained  in  section  3.1.  All  units  are  in  inch*  pound, 
seconds,  and  radians. 


FOBtRAN  iV  G  LEVEL 


MAIN 


DATE  •  711*0 


ll/S*/5* 


RA*E  0001 


0001 

0002 

0003 

0004 

0005 
0006 
0007 
0008 
0009 
0010 
00  U 
0012 
0013 
0014 
0015 
0016 
0017 
0018 
0019 
0020 


0021 
0022 
0023 
0024 
0025 
0026 
0027 
0028 
0029 
00  30 
0031 
0032 


C  ELEVEN  OOF  NLOEL  OF  MHO  TQM  STREETER  JAN  1977 

C 

C 

C  THE  kaIK  program  ESSENTIALLY  READS  AND  NAMES  DATA.  CONSTANTS 

C  ARE  1NITULI2E0  IN  DERIV  AND  SUBROUTINE  AUTTA  INTEGRATES  THE 

C  CIFFERENTUL  EQUATIONS. 

C 

COMMON  /OATAl/T IME  tTIMEH,TlMEH2  .TINcH8 
COMMON  /OERI  Vl/Q{m,QDlll),QDDIll  » 
common  /KUT  r i3/l tPTS  ,  IRPTS.IGPTS 

COMM  jN  /XR  =  ADI /BRCHX (105)  ,BR CmY ( 10 5  ) . RQDX ( 105  I  ,ROOV(105I  . 

I  CAMMAXI1051  .(.AMMAVHOS) 

COMMON  /NAME9/H0LR1 

MEAD  I  ,  IBPTS  .1 BRCHXI  1  I  ,eRCHY4  I  )  .iMl.IfiPTS) 

READ  l,lfiPTS,(ROCX(l),RDDY()),|«l.lRPTS) 

READ  1  .!GPTS,(GAMMAX( 1  I  ,G AMMAY ( I ) , I •! t ICPTS  ) 

1  FORMAT H 10/ief 10  .0)  I 
DO  15  Ml.IGPTS 

15  GAMMAXHI-GAMMAXd  1*3.14159/180. 

F1NI*.005 

HOORl-1. 

HOORl-0. 

TIHEH2.TIMEH/2 . 

TIMEMe«TIMEH/e. 

PRINT  2 

2  FQRNATI6X,4HTIME/) 

PRINT  3 

3  FORMAT  (6XMnETAH1XMHY,1  IX, 1KX,11X,1HY,11X,1HZ»  9X«3HPNI, 

1  9X,3HGAM,  9X,2HNU,  9X,4MTHET,  9X,3HPSI.  9X,3MTAU» 

C 

C  INITIALIZE  CONSTANTS  IN  SUBROUTINE  OERIV 
C 

CALL  OERIV 
20  call  MUTTA 
PRINT  10, TINE 
PRINT  10,(0(11  , ft  ,111 
PRINT  10,IQDIM,fl,ll) 

PRINT  10, (000(1), f 1,11) 

10  F0RMAT(HP12.4) 

PRINT  11 
It  FORMAT)//) 

IF(TIME  .LE.  FINl  )  CO  TO  20 

STOP 

END 


0001 


0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 

0011 

0012 

0013 

0014 

0015 
0016 
0017 
00  IB 
0019 
0020 
0021 
00  22 
0023 
0024 
0025 
0026 
0027 
0028 
0029 
0030 
0031 
00  32 
0033 


00  34 
00  35 
0036 
0037 
00  38 


0039 


SUBROUTINE  OERIV 

C  this  SUBROUTINF  initializes  constants  and  defines  the  IMBEDDED 

C  TERMS,  ZZ(I),  NHKH  HERE  REMOVED  FROM  THE  PARTIAL  DCRIVMIVES. 

C  terms  are  used  in  DERI  AND  DER2. 

C 

COMMON  /DER1V1/Q(11).Q0(UI  ,000(11  ) 

COMMON  /DcRlV2/Al ,A2,A3,ASTAR,XMY1 ,XKY2 

COMMON  /OERI V3/A1SUB,A2SU6  ,A3SU8,XKK1 ,XKK2 , A1 BAR 

COMMON  /DER1V4/XL(2.2),XM(2,2)  .XN(2,2),XK(2,2I 

COMMON  /DERIV5/tl,E2,E3,fFl,fF2,FF3,01,02,03,XU,ETAl,ZETA| 

COMMON  /DERI V6/B1 .B2,83,61BAR,XI ,ZETA 

COMMON  /OERI V7/PC (6,12,12  ,3) 

COMMON  /OERlVe/XIB,EB,Z6TA6,XlR,£R,ZETAR,XlC,EC,Z£TAC 
COMMON  /DERIV9/8LFT  ,COFT,RDFT  ,FOFG 

COMMON  /NAHEI/PT (8,11 ,12,31 ,PM (5,12,12,3) ,PU(4,U ,12,41 
1  ,PD( 3  .1 1 .12,4) 

COMMON  /NAME  7/DA (1  ,I  ,1 ,4) ,DP (1,1,1 ,4) 
common  /XZ290/ZZ  190) 

COMMON  /DERCDN/  G 1  ,G2  ,C3 ,  OIP  ,  02P  ,D3P , C 1 P , C2P ,C3P , ALPH A I , 

I  ALPHA2 ,ALPHA3,HH1 ♦MM2,MH3 
COMMON  /TRIG/C6,S6,C7,S7,CB,SB,C9,S9,C10,S10,C11.S11 
AST AR> AST AR*3. 14159/ 180. 

GfO. 

G2--125.1  -  54.193fSIN)ASTAR  ) 

&3--7.15  -  54.193fCaSl  ASTAR  ) 

01PO9.5 

D2P.54.193f  SIN(  ASTAR)  -  54 .0327 ‘S  I  N(  AST  AR  ♦  0.2626) 
D3P*54.193fC0S(  ASTAR)  -  54. 0327*COS(  ASTAR  ■*  0.2*26) 

CIP-DIP 

C2P-02P 

C3P»D3P 

ALPHAl.39.5 

ALPMA2«-125.1  -  54  .0327«>S1N(ASTAR  •  0.2626) 

ALPHA3*-7.15  -  54.0327*C0S(ASTAR  *  0.2626) 

HHl *0. 

hH2-54.l93fSlX(  ASTAR)  -  2  .*54 .0327/3  .‘SIN)  ASTAR  ♦  0.2626) 
HH3«54.1931*C0S( ASTAR)  -  2  .  *54 .032 7/3 . *COS I  AS T AR  ♦  0.2626) 
RETURN 

entry  oerfuc 

C 

C  ZERO  ALL  DERIVATIVE  FUNCTIONS  W/R  TO  KINETIC  ENERGY 

C  (NOT  INCLUDING  THE  ANGULAR  TERMSI 

C 

DO  I  f 1 ,8 
DO  1  j>i,n 
00  1  K«1.12 
CO  I  L»l,3 
1  PT ( I , J,R,L )*0. 

C 

C  ZERO  KINETIC  ENERGY  ANGULAR  TERMS 

C 

DO  12  1-1. 5 


F-3 


THESE 


0040 

0041 

0042 

0043 


0044 

0045 

0046 

0047' 

0046 


0044 

OCSO 

0051 

0052 

0053 


00  54 
0055 
00  56 
0057 
0051 
0054 
0060 
0061 
0062 
0063 
0064 
0065 
0066 
0067 
0068 
0069 
00  70 
00  71 
00  72 
0073 
00  74 
0075 
00  76 
0077 
0076 
00  79 
0080 
0061 
0082 
0083 
0084 
0085 
00  66 
0087 
0088 
0089 
0090 
0041 
0092 
0093 
0094 
0095 
0096 
0097 
0098 
0099 
0100 
OlCl 
0102 
0103 
0104 
0105 
0106 
0107 
0108 
0109 
0110 
oin 
0112 
0113 
0114 
0115 
0116 
0117 
0116 
0119 
0120 
0121 
0122 
0123 
0124 
01*5 
0126 
0127 
0128 
0129 
0130 


01  31 


0132 


C 

c 

c 


c 

c 

c 


DO  12  J«1 .12 
DO  12  K-1 .12 
UO  12  l.«l,3 
12  1  .J.K.L )-0. 

Z£RO  DERIVATIVE  rUKCTlOHS  M/*  TO  POTENTIAL  ENERGY  U2 

OG  4  1-1.4 
OU  4  J*1«I1 
DU  4  H<  1  .12 
011  4  L*lt4 

4  PUI I . J.K.L) >0. 

£ERD  derivative  functions  n/r  to  dissipative  energy 

DC  5  1-1.3 
OU  5  J-l.U 
DO  5  K-1  .12 
DO  5  L-1.4 

5  POn  ,  J,K.L)-0. 

ZERO  DERIVATIVE  FUNCTIONS  N/R  TD  GENERALIZED  FORCES 


DO  6  1-1.6 

DO  6  J-lt}2 
DO  6  K-1, 12 
00  6  L«1  .3 
6  POM  .  J.K.L  1-0. 

C6*C0SIQt6n 

56- SIN(Q(6»I 
C7-C0S<QI7M 

57- S1N(Q(7M 
C8>C0S(Q(8> I 

58- SIN<0(«n 
C9>CUS(0(9) I 
S9>S1NI0<9M 
ClO-COSIOtlOl > 
SlO*5INI0M0n 

cii-cosiQiin ) 

Sll-SINIQMin 
ZZM)-C8-S9  ♦  S8*C9»SI0 
ZZI21--C8-S9  -  S8*C9*S10 
2ZI3>-Ce*C9*S10  -  S8*S9 
ZI141  .-C8-C9*510  ♦  S«*S9 
ZZ1S1-C8OS96SI0  ♦  5B-C9 
ZZ{6)--C8*S9*S10  -  S8*C9 
ZZ<7)-C8*C9  -  S8-S94S10 
ZZI8)--Cfl*C9  ♦  SB*S9*S10 
ZZI91-C9-S10-S6  ♦  S9-C6 
ZZ  MO  1 --C9-S  104S6  -  S9-C6 
ZZM  n-C9«S10-C6  -  S9*S6 
ZZM2»--C94S10*C6  ♦  59*56 
ZZM3)-C9*S6  ♦  S9*S!0*C6 


ZZM4I--C9*S6  -  S9*S10*C6 
ZZM5>-C9*C6  -  S9*S10*S6 
22M6I--C9-C6  ♦  59'»S10*S6 
ZZ M 7  I  p-Cl 1 •C9*S10  -  S11*C9*C10*C6 
ZZM8>-CM*C9*S10  ♦  SI  1  *C9*C10*C6 
Z2M9|p-Cll*C9*C10*C6  ♦  S11*C9*S10 
ZZ120  )-Cn-C9«»ClC*C6  -  S11*C9*S10 
2ZI21  I-CII^CIO  -  SU*S10*C6 
Z2I22 >--Cll*C10  ♦  511*S10*C6 
ZZ(23»»-:il*S10*C6  -  Sll*C10 
ZZ(24)-C11*510*C6  ♦  Sll*C10 
ZZI25)«C11*C10*C6  -  S11*S10 
ZZ1261— Cll*C10*C6  ♦  SU*S10 
Z21271«-Cn*S10  -  S11*C10*C6 
ZZI28I-CU-S10  ♦  S11*C10*C6 
ZZ (291 -Cll -S9*S1C  ♦  S11*S9*C10*C6 
ZZI30  »--CU*S9*S10  -  SU*59*C10*C6 
ZZ  (  31  I  «CU  *S9*C10*C6  -  511*59*510 
ZZI32)'-Cll*S9*C10*C6  ♦  511*59*510 
ZZ133l--Cll*C9*ClO  ♦  Sll*C9*S10*C6 
ZZM4I-C11*C9'»C1D  -  S11*C9*510*C6 
ZZU5MC1 1  »C9*S10*C6  ♦  511*C9*C10 
Z2t36I--Cll*C9*S10*C6  -  S11*C9*C10 
ZZ137J-C11*59*C10  -  5M*S9*S10*C6 
ZZ  I  38  I  »-Cl 1 *59*0 10  ♦  51 1 *59*S10*C6 
ZZ(39J*-CU*S9*S10*C6  -  S11*S9*C10 
ZZ  140  I  *C U *59*510»C6  ♦  Sll*S9*C10 
21(41 ) »A2SU6*C6  -  A35UB*S6 
ZZ (42  ) - A2SUB-S6  ♦  »3SUB*C6 
2ZI43I--&3  ♦  0(51  ♦  22(421 
Z2(44»--C2  -  0(21  ♦  0(41  ♦  ZZ(4ll 
ZZI451--D2P  ♦  22(431*56  *  Z2(44l*C8 
ZZ(46l«-D3P  ♦  Z21431*C8  -  12(441*58 
ZZ(47I--D1P  -  G1  ♦  A1SU8 
ZZ(48)--02P  -  G2  ♦  A2SU8 
ZZ(49I*-D3P  ♦  A3SUB  -  &3 

2Z  1501-ZZ  (47  l'’*2  ♦  ZZ(45M*2  ♦  ZZ(461**2 
Z2(511«ZZ(471-*2  ♦  Z2(4B1*»2  *  ZZ(491**2 
ZZ(52)--ZZ(50}**C.5  ♦  ZZ(511**0.5 
22(53  1-Z2(52)‘‘XK81*ZZ(  50  !•*(-. 51 
ZZ(54)-D1P  ♦  AIBAR  -  Cl 
2Z(551--G2  -  02P  ♦  A2SUe 

Z2(561-ZZ(54 1 **2  ♦  ZZ1451**2  ♦  ZZ(461**2 
22(57I-ZZ(541-*2  ♦  22(551**2  ♦  ZZ(491**2 


ZZ158)--2Z(5tt**C.5  ♦  ZZ(571**0.5 

22(  59  1  >XKK2*ZZ (581*221 56  !••(-. 51 


22(60  1--XNM  .1  1 

1  XNM.l  l*2ZM3) 

2  *1*S9*C10 
Z2(61 1*-XN(1,2I 

1  XNM  .2  I  *ZZM5  1 

2  ♦  A2*S9*510 
22(621  --XN(2 ,1  1 


-  A3  ♦  XNM, 11*22(151  -  XLM,11*S9*C10  * 

♦  A2*59*510  ♦  A3*C9  ♦  0(41*59*510  ♦  Q(51*C9 

-A3  -  XLM  ,21*S9*C10  ♦  XM  M  ,  2  1  *22  (  U 1  ♦ 

♦  C(4}*S9*S10  ♦  Q(51*C9  -  A1*S9*C10  ♦  A3*C9 

-  A3  -  XL(2,.ll*59*C10  ♦  XM  ( 2 . 1  I  *ZZ  (1 3  1  ♦ 


F-4 


C(*.)«S9»S10  ♦  i;i5l*C9  -  *1*S9*C10  ♦  A3*C«  ♦ 


01  i3 


013<> 


ons 


0136 

0137 

0136 

0139 

01<i0 

0141 

0142 

0143 

0144 

014S 

0146 

0147 

0146 

0149 

0150 

OlSl 

0152 

0153 

0154 

0155 

0156 

0157 

0156 

0159 

0160 

0161 

0162 

0163 

0164 


1  X!<f  2,11-1205)  ♦ 

2  *2*S9»5I0 

Z2t63)*-xNI2,2l  -  i3  ♦  XH  1 2 , 2  1  *2 2 1 1 3 1  *  XMI 2 . 2  I *22  1 1 5 1  - 

1  XL12,2)='59*C10  ♦  C(4)»S9*S10  •»  ai5l«:9  -  61*S9*C10  ♦  A3»C»  ♦ 

2  *2*S9»5lO 

22 (64 1 1-ALPMA2  *02*  0(2)  ♦  42  ♦  Cl*l-C2  -  012)  -  *2  ♦ 

1  41PH62*C1C  ♦  ALFHi3*59<>SlO  -  C9*S1C»ULPH*1  -  0(3)1)  ♦ 

2  S8»(-C.3  -  hi  *  ALPH*3»C9  ♦  S9*(4LPhll  -  COII) 

22(65 )<-ALPHA2  *  02  *  0(2)  «  A2  «  CI*l-02  -  0(2)  -  hi  * 

1  ALPMA2<'Cn  ♦  ALPHA3»S9»S10  -  C9*$1C»  l-ALPMAl  -  0(3)))  ♦ 

2  S8-(-A3  -  C3  ♦  ALFHA3-C9  ♦  S9»(-ALPHil  -  3(3))) 

Z2(66)»“02  -  0(2)  ♦  0(4)  ♦  B2*C6  >  63*56 
22(67). -03  4  0(5)  *  62*56  *  B3*C6 

22(63). -CIP  4  01  -  01 

22(69). -C2P  4  Z2(67)*S8  ♦  2Z(66)*C6 

22(70). "CSP  ♦  22  167)*C«  -  22(661*56 

22(71 ).(22170)**2  ♦  22{6a)**2  ♦  22 (69 ) **2 ) ♦• ( - .5 ) *0 . 5 
22(72 ).-32*S6  -  fc3*C6 
22(73)-62*C6  -  03*56 

22(74). (Z2(7&)**2  *  22(69)**2  ♦  161013  4  CIP  -  61 ) **2 ) • • ( - .5 ) *0 .5 

22(75). -02  -  0(2)  4  C(4)  -  A35Ue*56  4  62SU|*C6 

22(76). -03  4  0(5)  4  A3SUe*C6  *  A25U6*56 

22  (77)— C2P  4  22  (76)*S8  4  Z2(7S)*C8 

22(78). -D3P  4  2Z(76)*C8  -  Z2(7S)*S6 

22<79).-A3SU8*S6  4  A25UB*C6 

22  (80  ).-A35ue»C6  -  A2SUB*56 

22(81 l«(Z2  (77)**2  4  22(7B1**2  4  |D1P  4  11866  -  61 ) **2 ) •• ( - .51 *0 .5 

2i(82).(22  (77)**2  4  22  {78)**2  ♦  (AISUS  -  OIP  -  61  I **2 ) •• ( - .5 ) *0 .5 

2Z(63)*A2*S9*S10  4  A3*C9  4  0(4)*S9*51C  4  0I5)*C9  -  AI*59*C10 

22(B4).-A2*C9*510  4  63*59  -  Q(4|*C9*510  4  Q(5)*59  ♦  61*C9*C10 

22(85I.2Z(12)*C11  -  C9*C10*S1I 

22(e6)"-C9*C10*Cll*C6  ♦  C9*SIO*511 

221871 .C10*C11*C6  -  510*511 

22 (88). -CIO *511*66-510*611 

22(89  I  .S9*C10*C11»C6  -  59*510*521 

2Z(90)*ZZ(13)*C11  4  59*00*511 

C6LL  DERI 

C6LL  0Eft2 

RETURN 

ENO 


0001 


0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 

0011 

0012 

0013 

0014 

0015 


0016 
0017 
0016 
0019 
0020 
0021 
0022 
0023 
0024 
0025 
0026 
0027 
0026 
0029 
0030 
0031 
00  32 
0033 
00  34 
0035 

0036 

0037 

0038 

0039 

0040 

0041 


5UBR0UTINE  DER)  5 

C 

C  TH15  SUBROUTINE  CEfINES  ThE  PARTIAL  DtRIVATIVES  PERTAINING  TO  THE 
C  TRAN5(.ATiaNAL  PART  OF  THE  KINETIC  ENERGY. 

C 

COMMON  /DER IVl/Q 1 11 )  ,00( 1 1) ,000111) 

COMMON  /DEP I V2/A1,A2  ,  63 , 6  STAR ,XHT1 .IRT2 

COMM  UN  /DERI  V3/A1  SUB  ,A2  sue,  A35U8,)(K«1  ,XKK2,A1  BAR 

CCMMDN  /0ERlV4/XL{2,2},XH(2,2>iXH(?.2>.RK(2,2) 

COMMON  /DER|V5/E1,E2,E3,FF1,FF2.FF3.:uC2.03,XI1,ETA1,2ETA) 

COMMON  /0ERIVb/31,B2,B3,BlBAR,)(|.2ct6 
COMMON  /DERIV7/P6(6,12,12,3) 

COMMON  /0ER]V6/X18,EB,ZETABtX(R.ER.2£TAR,XICiEC,2ETAC 
COMMON  /1JER1V9/BLFT,C0FT,R0FT,F0FG 

COMMON  /NAMEl/PT  l8,ll,12,3),PH(S.U.12»3I.PU(4,ll,12t4) 

1  ,PD( 3 ,11 ,12  .4) 

CUMMQN  /NAHE7/DA(1,1,1,4)  ,DP(l,l,l,4l 
COMMON  /X2290/Z2 (90) 

COMMON  /DERCUN/  C 1 ,02 ,6 3 ,0 1 P  , :.2P  , 33P  , C 1 P  , C2P  ,C3P  ,  At  PHAl  , 

1  AiPHA2  ,ALPHA3 ,HH1 ,HH2 ,HH3 
COMMON  /TRI6/C6, S6,C7,S7,C8,Se,C9,S9,ClO,S10,Cll,Sll 
C 

C  8ER1VATIVE  FUNCTICNS  N/R  TO  K.E.  ISCT  INCLUDING  ANGULAR  TERMS) 

C 

PTI  1 ,  3,12,  1  )  .  1 

PT(  1,  9,12,  1)  .  A2*S9*510  4  A3-C9-11 *C10*59 

PTI  1,10,12,  n  .  -A2*C10*C9-A1»C9»S1Q 

PTI  1,10,12,  2)  .  -A2*5lO  4  Al»C10 

PT)  1,10,12,  2)  .  -A2*S10  4  Al*tlO 

PTI  1,  9,12,  3)  «  A2*C9*S10-A3»59-il.C10*C9 

PTI  1,10,12,  3)  *  A2*C10*S9  4  Al»59.5'10 

PTI  2,  4,12,  II  «  -C9*S10 

PTI  2,  5,12,  1)  .  59 

PTI  2,  9,12,  11  .  0(4)«S9*S10  4  a(5)»C9 

PTI  2,10,12,  11  «  -Q14)*C10*C9 

PTI  2,  4,12,  21  .  CIO 

PTI  2,10,12,  2)  «  -Q(4)*S10 

PTI  2,  4,12,  31  .  S9*S10 

PTI  2,  5,12,  3)  .  C9 

PTI  2,  9,12,  3)  .  Q(4)*C9»S10-C(5)  »S9 

PTI  2,10,12,  31  •  0{4>*C10*59 

PTI  3,  2,12,  1  )  .  -C9*510 

PTI  3,  8,12,  1)  .  HM2*22(1)  ♦  hm3«22I)I 

PTI  3,  9,12,  1)  .  -HH1*C10*S9  4  •<f-2»Z2(5)  4  HH3*zz(7I  4  G3*C9-C1*C 

110*S9  4  S9*S10*IC2  4  0(2) ) 

PT(  3,10,12,  1)  .  -HHl*C9*S10-t-M2»Ci:*C2*C9  4  HH3 -Cl  0  <  58  »C  9-0  1  *C9* 

1S10-C10*C9*I02  4  CI2)) 

PTI  3,  2,12,  2)  .  CIO 

PTI  3,  8,12,  2)  *  -HH2’C10*S8-KH39Cl7=’Ca 

PTI  3,10,12i  2)  .  t  Hi  *C10-HH2''Ce»>i:  -  hH3*58='S10  ♦  0 1  *C  10- S 1 0  *  ( G2 

1  4  0(21) 

PTI  3,  2,12,  3)  .  S9*S10 

PTI  3,  8,12.  3)  «  HH2*22(7)  4  hh3*22(«>) 


F-5 


00^2 

PT( 

3.  9,12i 

31 

* 

-rtHl*C10*C9  ♦  NH2*ZZC3I  ♦  MH3*ZZC21-6l»5*-Gi*ClO 

1  *C9 

♦  C9»S10‘>CC2 

Qf2H 

00>.3 

PI  1 

3.10,12. 

3) 

■ 

hh1*S9*S10  ♦  HH2*C10*C8*S9-HH3*Cl0*S9*Si  •  tl*S» 

l*S10  ♦  C10'-S9*IC2 

♦  CC2n 

0044 

PT< 

1  .  9.  9, 

1  ) 

A2*C9»510-A3*S9-*l*ClO*C9 

0045 

PT( 

1  .  9,10, 

1  1 

A2*ClCi*59  «  81*59*510 

004b 

PH 

1,10,10, 

1  1 

82*C9*510-Al*C10*C9 

0047 

PT  ( 

1.10.10, 

2  ) 

-A2*C10-A1*S10 

0048 

PT  1 

1,  9,  9, 

3) 

-A2*S9*S10-A3*C9  ♦  Al*Ct0*59 

0049 

PH 

1,  9,10, 

3) 

A2*C10*C9  ♦  A1*C9*S10 

0050 

PH 

1  ,  10,10, 

3) 

-82*59*510  *  AI*C10*S9 

0051 

PH 

2.  4.  9. 

1  ) 

59*510 

0052 

PH 

2,  4,10, 

1  1 

-C10*C9 

0053 

PTl 

2,  5,  9, 

1  1 

C9 

0054 

PH 

2,  9,  9, 

1 ) 

CC41*C»*510-0I5I*S9 

0055 

PH 

2,  9,10, 

1  ) 

CC4>*C10*5« 

OOSb 

PH 

2,10,10, 

1  ) 

QI4I*C9*S10 

0057 

PT< 

2,  4.10, 

21 

-510 

0058 

PTC 

2,10,10, 

2> 

-0C4I*C10 

00  59 

PH 

2,  4.  9. 

31 

C9*S10 

OObO 

PH 

2.  4,10, 

3) 

C10*59 

0041 

PH 

2.  5.  9. 

3) 

-59 

0042 

PH 

2,  5,  9, 

31 

-59 

0063 

PTC 

2,  9,  9, 

3) 

-g(4>*S9*510-Q15l*C9 

0044 

PTC 

2.  9,10, 

3  1 

0(4)*CIO*C9 

0045 

PTC 

2.  9,10, 

3> 

ec4i*cio*c* 

0044 

PTC 

2,10,10, 

31 

-eC4)*59*510 

0047 

PTC 

3,  2,  9, 

1  1 

59*510 

0048 

PTC 

3.  2,10. 

1  ) 

-C10*C9 

0049 

PTC 

3.  8.  8, 

1  ) 

hH2*Z2C31  ♦  HH3*ZZI2) 

0070 

PTC 

3.  8,  9. 

1  > 

HH2*ZZC7)  *  HH3*ZZC41 

0071 

PTC 

3.  8,10, 

1  ) 

HH2*C10*Se*C9  ♦  HM3*C10*CI*C9 

0072 

PTC 

3.  9,  9, 

1  ) 

-HH«>C10*C9  ♦  HH2*ZZ(3J  ♦  HH3»ZZC2  >-C3*59-61*C10 

l*C9 

♦  C9»SI0'»CG2 

012  n 

0073 

PTC 

3,  9,10, 

1  ) 

- 

HMl*59*StO  *  HH2*CI0*C8*S9-HN3*C10*S«*5I  *  (1*59 

MSIO  ♦  C10*S9*C62 

♦  0C2H  ' 

00  74 

PTC 

3,10,10. 

1) 

■ 

-MM1*C10*C9  ♦  HH2*C«*C9*SIO-HM3*5I*C9*S10->61*CIO 

1PC9 

♦  C9»S10«(G2 

0121) 

0075 

PTC 

3,  2,10, 

2) 

-SIO 

0074 

PTC 

3.  8,  8, 

2) 

-MH2*C10*C8  ♦  HH3*C10*S8 

0077 

PTC 

3,  8.10, 

2) 

HH2*5e*S10  ♦  HH3*C8*$10 

0078 

PTC 

3,10,10, 

2  ) 

-Mhl*51C-HH2*C10*C8  ♦  HH3*C10*S8-G1 •510-C10*CG2 

i«  oc2n 

0079 

PTC 

3.  2.  9, 

3) 

C9*510 

0080 

PTC 

3,  2,10, 

31 

C10*59 

0081 

PTC 

3,  8,  8, 

31 

HK2»ZZ16)  *  HH3*2ZC«t 

0082 

PTC 

3,  8,  9, 

3) 

hH2*Z2C2)  *  KH3*ZZC4I 

00^3 

PTC 

3,  8,10, 

3  > 

-mh2 •Cl  0*59*51 -HH3*C10*C8*S9 

0084 

PTC 

3,  9,  9, 

3  I 

KHl*C10*S9  ♦  NH2*ZZC4}  ♦  HH3*Z2Ce>-C3*Cf  ♦  H*CI 

10*S9-S9*510*CG2 

QC21 1 

0085 

PTC 

3,  9,10, 

3  1 

NH1*C9*S10  ♦  MH2*CIO*C8*C9-MH3*C10*S8*C9  ♦  G1*C9 

1*S10  ♦  C10‘>C9*CC2 

*  0C2)> 

0084 

PTC 

3,10.10, 

3) 

• 

HHI*CI0*S9-HH2*C8*S9*510  ♦  HH3*S9*S8*S10  ♦  «l*Cl 

10*$»“S9*S10*C&2 

ccati 

0087 

PTC 

4,  4,12, 

1  ) 

E2*ZZC9»  ♦  E3*ZZCU) 

0088 

PTC 

4,  9,12, 

1  T 

E2*2ZC13l  ♦  E3*ZZC15I-E1*C10*S9 

0089 

PTC 

4.10.12, 

1  1 

-£2*C10*C4*C9  ♦  E3*C10*54*C9-EI*C9*S10 

0090 

PTC 

4,  4.12, 

2  ) 

-E2*C10*56-E3*C10*C4 

0091 

PTC 

4.10,12. 

2) 

-£2*C4*S10  ♦  E3*54*S10  ♦  E1*C10 

0092 

PTC 

4.  4,12, 

3) 

E2*2ZC15)  ♦  E3*ZZC14I 

0093 

PTC 

4,  9,12, 

3  ) 

E2*ZZCin  ♦  E3*2ZC  10I-E1*C10*C9 

0094 

PTC 

4,10,12. 

3  1 

E2*C10*C4*59-E3*C10*S4*59  ♦  E1*59*S10 

0095 

PTC 

4,  4,  4. 

1  1 

f2*ZZClU  ♦  I3*22C10» 

00  94 

PTC 

4,  4.  9, 

1  1 

E2*ZZC15I  *  Ea*ZZC14) 

0097 

PTC 

4,  4,10, 

1  ) 

E2»C10*S4*C9  ♦  E3*C10*C6*C9 

0098 

PTC 

4,  4,10, 

1  J 

e2*C10*56*C9  ♦  E3*C10*C4*C9 

0099 

PTC 

4,  9,  9, 

1  > 

E2*ZZCin  ♦  E3*Z2C10I-E1*C10*C9 

0100 

PTC 

4,  9,10, 

1) 

E2*C10*C4*S9-E3*C10*S4*59  ♦  E1*$9*S10 

0101 

PTC 

4,10,10, 

1  1 

E2«C6»C9*510-£3*54*C9*S10-E1*C10*C9 

0102 

PTC 

4,  6,  6, 

21 

-e2*C10*C4  ♦  E3*C10*56 

0103 

PT  ( 

4.  6,10, 

2  ) 
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PTC 
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3> 

-£2*C10*54*S9-E3*C10*C4*59 

0108 

PTC 

4,  9,  9, 
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PT  C 

5,  9,12. 

3  ) 
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1*C- 

CI1*C10*C9  ♦ 

$lloZ2<ll  ) ) 

0120 

PTC 

5,10,12, 

3) 

FF2*ZZC 31  1-FF3*C10*54*59  ♦  FF1*ZZ(291 

0121 

PTC 

5,11,12, 

31 

FF2*(Cll*C10*S9-Sll*ZZCl3n  ♦  f  FI  *(CU  *22  (13  1  ♦ 

IS11*C10*S9) 
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PTC 
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PTC 
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FF2*CU •C10*$6*C9  ♦  FF3*C10*C4*C9  ♦  FFI*S11*C10* 

154*C9 

0125 

PT  C 

5,  4.11, 

1  > 
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1S1MC104S9) 
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PTC 

5,10,10, 

1  > 

B 

FF2*22C351-FF3*54*C9*S10  *  FF1*Z2C33I 

0130 

PTC 

5,10,11, 

1  ) 

B 

FF2*22(ia}  ♦  FF1*ZZ(191 
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19-SlloZH12n 

PT(  5,  6,  6,  ■  -f  F2<-C1  1 '>C10«C6  ♦  f  F  3*C10*S6-8F1  *511  *C10*C6 

PT(  5,  b,10i  2»  •  FF2 •Cl  1 ‘Sb^SlO  •  FF3*C6*S10  ♦  FFl *51 1 4$**$10 
Pit  5,  6rll.  2)  •  FF2*Sn  •C10*S6-FF1*C1MC10*S6 
pn  5,10.10,  2>  -  FF2*Z2I26>  ♦  FF3*C1C»S6  ♦  FFl»Zt<271 
PTt  5,10,11,  21  •  FF2*ZZ(22)  ♦  FF1»ZZI231 

pn  5,11,11,  2)  *  FF2*ZZ«261  •  Ffl*ZZ(271 

PTI  5,  6,  6,  31  «  FF2*C11*ZZI14I  ♦  FF3>ZZ(16>  ♦  FF 1 *51 1 *22 116 1 

PT(  5,  6,  9,  3)  -  FF2*C11 ‘ZZIIOI  ♦  FF3»ZZI121  ♦  FFl •$! 1 *22 1 10 1 

PTt  5,  6,10,  31  •  -FF2®cmClC'>S6»S9~Ff3»C10*C6*S9-FFl*Sll*C10*S6* 
159 

PTt  5,  6,11,  31  »  *FF2*S1I«Z21151  ♦  FF 1 *0 1 1 *2 2 ( 15 1 

PTt  5,  9,  9,  3)  -  FF2*tCl  l*2Ztl4l-Sll  •C10*59)  ♦  FF3»ZZ(U1  •  Ffl*l 
K11‘>C10«>S9  ♦  Sll«2Zfl41l 

FT(  5,  9,10,  3)  *  FF2*22t201-FF3*C10*S6*C9  ♦  FF1*Z2(18) 

PTt  5,  9,11,  31  •  FF2*tCll*C10*C9-Sll*22111 11  •  FFI • 1  Cl  1 *22 1 II  I  ♦ 
1S11->C10*C91 

PTt  5,10,10,  3)  >  FF2*ZZ(39l  ♦  FF3''56*S9*S10  ♦  FFl*ZZt37| 

PTt  5,10,11,  3)  •  FF2*ZZ(301  ♦  FF|*ZZ<3n 

PTt  5,11,11,  31  ■  FF2*(-Cn*ZZn3l-Sll*ClD*S91  ♦  Ff  1  •  I C 1 1*CI0*S9-S 
lll*Z2tl3n 

PTt  6,  6,12,  It  «  C3*ZZ(m  ♦  01»Sll*ZI19l  ♦  t)2*C  1 1 'ZZ 1 9 1 

PTt  6,  9,12,  II  «  03*221151  ♦  01 • ( -C 1 1 *C 10*S9  ♦  $11»ZZ{13I1  ♦  02*4 

1C11»ZZ(13I  ♦  S11*C10*S9) 

PTt  6,10,12,  11  «  03*C10*S6*C9  ♦  01*224171  ♦  02*224191 
PTt  6,11,12,  II  •  01«4CU»ZZ412I-5U*C10«C91  ♦  02*4-Cll  •ClO*C9-Sll 
1*ZI1121) 

PTI  6,  6,12,  21  *  -D3*C10*C6*01*Sll»C10*56“D2*tll*ClO*S6 
PTt  6,10,12,  21  •  03*Sb*SiO  ♦  01*ZZt21l  ♦  02*211231 
PTt  6,11,12,  2)  »  Dn2Z1251  ♦  02*221271 

PTt  6,  6,12,  31  •  03»ZZtl4l  ♦  01 *51 1 *ZZ 4 15 1  ♦  D2«C 1 1 *22 4 15 1 

PTI  6,  9,12,  3)  •  C3*ZZ4101  ♦  0 1 • I -C 1 1 •€ 1 0*C9  ♦  511*22(111)  ♦  02*4 

ICll *22  4  1  1  1  ♦  511 •C10*C9  1 

PTt  6,10,12,  3)  •  -03*C10*$6*S9  ♦  01*22(29)  ♦  02*22131) 

PTt  6,11,12,  31  •  Dl*4Cll*22tl3)  ♦  S11*C10*59)  ♦  D2*4CH *010*59-51 
11*22(13)) 

PTI  6.  6,  6.  n  •  03*22  1  10  1  *  01  *5 1  1  *2 2 1  1  1  )  *  02*01  1  *22  4  11  ) 

PTt  6,  6,  9,  11  .  03*22114)  ♦  D1*S11*:21151  ♦  02*011*22415) 

PTt  6,  6,10,  1)  •  03*010*06*09  ♦  01*511  *0 10*56*09  ♦  02 *0 1 1 ♦C10*$** 
109 

PTt  6,  6,11,  1)  ■  D1*C11*22(9)-D2*511*22I91 

PTt  6,  9,  9,  1)  •  D3*22tl01  ♦  01*1-011*010*09  ♦  511*22(1111  *  02*( 
1011*22111)  ♦  511*010*09) 

PTt  6,  9,10,  II  •  -03*010*56*59  ♦  01*22(291  ♦  02*22(311 
PTt  6,  9,11,  1)  •  01 • (Oil *22( 131  ♦  511*010*591  ♦  02* ( 0 1 1 *0 10*59-51 
11*22(1311 

PTt  6,10,10,  n  •  -03*56*09*510  ♦  01*22(33)  ♦  02*22(351 
PTI  6,10,1),  1)  •  Cl*22{19)  •  02*22(16) 

Pn  6,11,11.  1)  •  Dl*(-Cll»ClO»C9-5U*22(12l  1  ♦  02*  ( -C 1 1  *22  ( 1 2  )  ♦ 
1511*010*09) 

PTt  6,  6,  6,  2)  •  03*0  10*  56-01*511  •Cl0*C6-02*Cn  *010*06 
PTt  6,  6,10,  2)  •  03*06*510  ♦  01*511*56*510  ♦  02*011*56*510 
PT(  6,  6,11,  2)  •  -01*0  11  *0  10*56  +  02  “511 *010 *56 

PTt  6,10,10,  21  •  03*010*56  ♦  01*22127)  •  02*22(26) 

PTt  6,10,11,  2)  •  01*22(23)  ♦  02*22(22) 

PTI  6,11,11,  2)  >  01*22127)  •  02*22(261 

PTI  6,  6,  6,  3)  •  03*ZZ(ltl  *  D1*S11*:2I14)  ♦  02 *0 1 1  *2 2 (  14 1 

Pit  6,  6,  9,  31  «  03*22112)  ♦  01*511*22(10)  ♦  02*011*22(10) 

pn  6,  6,10,  3)  «  -D3*Cl0*C6*i9-Dl''511*C10*56*59-02*Cll*Cl0*S6*S9 
pn  6,  6,1),  3)  «  Dl*Cn*ZZ{lS)-D2*5ll*2Z(151 

PT(  6,  9,  9,  3)  «  03*22(  16)  *  01  *(01  1  *010*59  ♦  511*22(141)  ♦  02*4C 
111*22(141-511*010*591 

PTt  6,  9,10,  3)  -  -D36C10»56»C9  ♦  01*22118)  *  02*22120) 

PTt  6,  9,11,  3)  *  01*1011*22111)  *  511*010*09)  ♦  02*4011*010*09-51 
11*22(11)} 

PTt  6,10,10,  31  -  03*56*59*510  ♦  01*221371  *  02*22(39) 
pn  6,10,11,  31  *  01*22(31)  ♦  02*2  21  30  1 

Pn  6,11,11,  3)  *  01*(C  1  1  *0  1  0*59-51  1  *22  (  13  1  )  ♦  02  »  ( -C 1  1  *2  2  4  1 3 1  -  511 
1*010*59) 

PTI  7,  6,12,  1)  :■  XI 1  *511  *22  (9 )  +  n*l*(C7*Cll*22(9l  ♦  57*22(1  1  1  ) 
1*  Z6TAl«(C7*2Z(in-S7»Cll*22l9)) 

PTI  7,  7,12,  1)  •  ETAl»tC7»2Z(9)-57*(Cll*2Z(12l-511*C10*C9)l  ♦  2ET 
1A1*(-C7*(C11«22(12 1-511 *C10*C9)-S7e22(9)) 

PTt  7,  9,12,  1)  -  X (1*1-011*010*59  ♦  511*22(13)1  ♦  ETfil • (C7*( 01 1 *2 
121131  ♦  511*010*59)  ♦  57*22  (  15  ))  ♦  2 H »1  * ( 07*22 1  1 5 1 -S 7* I  0 1 1 *2 2 ( 1 3 1 

I  ♦  511*010*59)1 

PTt  7,10,12,  1)  «  X11*2Z117)  ♦  Er41*(C7»22(19)  ♦  57*010*56*091  ♦  2 
1£TA1*(C7 *010*56*09-57 *22(19)1 

PTt  7,11,12,  1)  «  Xll*(Cll*2Z(12)-5n*C10*C9)  ♦  ET  Al  *C  7*  (-Oil  *010* 
109-51  1  *22  1  1  2  )  )-2  ETAl  *57*1 -Cm«“  10*19-5  11*221  12)1 
PTt  7,  6,12,  21  »  -XI1*511*C1C*$6  ♦  E T A  I  *  I -07*0  1 1 *C 1 0  *  56-57*0 1 0*06 
1)  ♦  ZETAl *(-C7*C 10*06  ♦  57*011*010*56) 

PTI  7,  7,12,  2)  •  tTAl*(-C7»C10*56-S7-Z2t25l>  ♦  2 E TA 1 • ( -07*2 2  I  25  I 
1*  57*010*56) 

PTt  7,10,12,  2)  .  Xn*ZZI2U  ♦  EIAl  1  07*22  (23  )  ♦  57*56*5101  ♦  Z£T* 
1 1*(C7*56*S10-S7*Z2I23) 1 

PTt  7,11,12,  21  »  xn*2Zr25l  ♦  Em  :  7*22  1  27l-2E7ai*S7*2ZI21l 
PTt  7,  6,12,  3)  »  xn*511*22tl5l  »  E  n  1  *  (  C  7  *C  1 1 '2  2  ( 1 5  )  ♦  57*221  14) 

II  ♦  ZeTAl*(C7*2Z(l4)-57*Cn*Z2(15)  i 

pn  7,  7,12,  31  *  ETA1»(C7»22I1SI-57*(C11*2Z(13)  ♦  511*010*59)1  ♦ 
12ETA1  *  l-C?”  (Cl  1*22  ( 13  1  ♦  511*C10*S9!-57-'2Z(15»  1 
PTt  7,  9,12,  3)  *  xn*(-Cll«C10*C9  ♦  511*22(11))  ♦  E  T  Al  >  (C  7*  ( C  U  *2 
12(11)  ♦  511*010*09)  ♦  S7*:2(lon  ♦  iE7Al*tC7»22<  1C1-S7*(C11*2ZI11) 
1*511 *010*09) ) 

PTI  7,10,12,  3)  »  xn*Z2(29)  ♦  E  T  A  I  - 1  0  7  *2  2  (  31  )  -  5  7  >C  10  »  56  *59 1  ♦  2ET 
1A1*( -07*010*56*55-57*22(31) I 

PTt  7,11,12,  1)  ,  XmtCll  22113)  ♦  511*010*591  ♦  ETAl  *0  7*(  01 1  *010 
1  *59-51  1  *22  (  1  3  )  )-2£TAl *57* toil *0  10*59-  51  1  *22  1  1  3  )) 

PTt  7,  6,  6,  1)  .  XI 1*511 *22(11  )  ♦  ETai*ic7»C11»22(11)  •  57*22  (10) 
I)  ♦  2ETA1*IC7“  22  (  10  1  -57*0  1  1  *22  (  1  1  ) I 
PTI  7,  6,  7,  1)  «  ETAl*(C7*ZZin)-57*Cll*22(9))  ♦  2ETA1  *(-07*011*1 
12(9)-S7*2ZI11 )  I 

PTI  7,  6,  9,  1)  .  X(1*S11*22(15)  ♦  t T Al » | 07 *0 1  1  *2 2  ( 15  I  *  57*22(  14) 
1)  ♦  2£TAl*tC7*22(14|-57*Cll*22(15H 
pn  7,  6,10,  1)  •  xn  *si  I  *010*56*09  ♦  ETAI  *(07*0  U*C10*S6*C*  ♦  ST* 
1010*06*09)  ♦  2ETA1 *(C7*ClC«C6*C9-S7*Oll*C10*S6*C9) 
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nT)  7,  b,n .  11  * 

1  Z  (6  ) 

PT(  7,  7,  11  ■=  ET41»(-C7«(Cn‘-Z2(12l-Sll'>£;iC*C9J-S7*ZZ(9n  ♦  It 
n41“(-C7  'ZZIVI  ♦  S7«'(CU=*ZZ<  12)-S1I«’C10*C9H 
PM  7,  7,  9,  11  *  f  T41*{C7»ZZ(15I--57«>(CU»ZZM31  •*  S11»C10»S9)»  ♦ 
IZPTAl*'!-':? -(Cl  1»ZZ(13)  ♦  S11'»C10*S9)-S7«'ZZ(1SI1 
PM  7,  7,10,  11  ■  PTAl  *(C7*Cir!'>56*C9-S7*ZZ  (19)  »  ♦  Z  E  T  *1  •  | -C  7*ZZ  <  19 
1  )-S7*C10oS6<-C9  ) 

PT(  7,  7,11,  1)  »  -ETAl«-S?<-(-Cll*C10oC9-5n*ZZ(12))-ZETAl*C7*(-Cn 
1*C10^C9-SU‘‘ZZ(UI  ) 

PM  7,  9,  9,  1)  *  XI1«(-CH*C10»C9  ♦  S11»ZZ11U)  ♦  ETAI<»(C7*1C11*2 
1Z(U)  ♦  S11«>C10»C91  ♦  S7i>ZZ(10))  ♦  ZeT61*(C7*2Z(10)-S7*(Cn*ZZ(ll» 

I  ♦  Sll*C10»C9)  I 

PM  7,  9,10,  1)  •  *ll*ZZ(29J  ♦  ET4I  *(C7*-2Z  (31  )-S7*C10*S6*S91  *  ZET 
lAl*(-C76C10*Sb*S9-S79ZZ(3ll) 

PM  7,  9,11,  1)  •  Xll*(Cll«»Z2(13»  ♦  Sll*CiO*S9j  ♦  ETAl  •C7*(C1 1  *€10 
!•  S9-SU0ZZ  113  >)-ZETAl*S7o(Cll*C  10*  59-511  *2  2(1  31) 

PM  7,10,10,  1)  •  xn*ZZI33)  ♦  £TA1»(C7*ZZ(35)-S7*56*C9*S10)  •  ZIf 
lAl*(-C7#S5*C9oSl 0-57*22 (35)) 

PT(  7,10,11,  1)  •  XI1*ZZ(19)  ♦  £T*1*C7»2Z)18)-ZET*1*S7*ZZ(1I) 

PM  7,11,11,  1)  •  XI  1*(-C11*C10«C9-SU*2Z  ( 12)  )  ♦  ETAl  7*( -Cl  1  *11  ( 
112)  ♦  Sn*C10*C9)-Z£T41«57‘(-CU*ZZ(l2)  ♦  Sll*ClO*C9) 

PTI  7,  6,  6,  2)  •  -Xll*Sil'»C10*C6  ♦  ETA1*(-C7*C11*C10*C6  ♦  S7*C10* 
156)  +  ZETaI  •(C7*(,10*S6  ♦  S7*Cl  1  •C10*C6  ) 

PT(  7,  6,  7,  2)  •  ETA1*(-C7*C10*C6  ♦  S7*C1 1 10*56 1  ♦  ZETA1*(C7*CI 
ll*C10*S6  ♦  S7*C10*C6) 

PTI  7,  6,10,  2)  •  XIl*5ll*S6*S10  ♦  ET A1 * ( C 7 *C 11 *56*510  ♦  S7*C**S10 

II  ♦  ZEIAl*(C7*C6*SI0-S7«cn*S6*S10  ) 

PT(  7,  6,11,  2)  -  -XI1*C11*C10*S6  *  ETA1*C7*5U*C10*56-IETA1*S7*51 
11 *C10*S6 

PM  7,  7,  7,  2)  •  ETa1*{-C7*ZZI25)  ♦  S7*C10*S61  ♦  lEt A1 *  I C7*C 10*54 

I  ♦  57*22(251 ) 

PT(  7,  7,10,  2)  •  ETAl*(C7*56*S10-S7*ZZ(23n  ♦  2f TAl •(-C7*ZZ(23»-S 
17*56*510) 

PT(  7,  7,11,  2)  •  -ETA1*S7*Z2(271-Z£TA1 *C7*ZZ (27) 

PM  7,10,10,  2)  •  XIM2Z127)  ♦  £T  A 1  *  ( C  7*ZZ  ( 26  )  ♦  S7*C10*S61  ♦  ZETA 
11*IC7*C10*S6-S7*2Z(26I ) 

PT(  7,10,11,  2)  •  XI1*2Z(23)  *  ET A  1 «C 7*ZZ ( 22 ) -ZET A 1 *S7*ZZ ( 22 1 

PT(  7,11,11,  21  •  XI1*ZZI27)  «  ETA1*C7*ZZ(26I-ZETA1*57*1Z(26) 

PT(  7,11,11,  2)  •  XI1*ZZ(27)  *  ETA1*C7*ZZ(26)-ZETAI*57*2Z(26) 

PTI  7,  6,  6,  3)  •  xn*Sll*ZZ(16|  ♦  E  T  Al  •  (  C7*C  It  *22  ( 14  )  ♦  S7*ZZ<14) 

II  ♦  ZETA1*(C7*ZZ  (161-S7*Cn*ZZ  (1611 

PT(  7,  6,  7,  31  •  ETA1*<C7*ZZ( I61-S7*C11*ZZ(1511  ♦  ZET Al •( -C7*C I 1 • 
1ZZ(15)-57*ZZ(16)  ) 

PTI  T.  6,  9,  3)  •  xn*SIl  *ZZI  101  ♦  ET  Al  *  (  C7*Cl  I  *Z  Z  ( 10  I  ♦  $7*221121 
1)  ♦  ZETA1*(C7*Z2 (12  l-S7*C  ll*ZZ  110)  1 
PTI  7,  6,10,  3)  •  -X  ll*Sn*C10*56*59  ♦  E  T  A 1  *  ( -C7*C  1 1  •€  10*S4*$«-S7* 
1CI0*C6*59)  ♦  2£TA1*(-C7*C10*C6*S9  ♦  ST*C 1 1 •£ 10 *54* 59  I 
PTI  7,  6,11,  3)  •  Xll*CU  •2Z(15l-tTAl*C7*Sn*Zlll51  ♦  ZETA1*ST*SU 
1 »ZZ 1 15  > 

PTI  7,  7,  7,  3)  -  ETI1*(-C7*(C1I*Z2« 13)  ♦  SI  I *C10*&« }-57*22 ( 1 4)1  * 
1  ZETA1*I-C7*Z2(15)  ♦  57 • I  Cl  1 *22  1 1 3  )  ♦  Sll*C10*S9)l 
PTI  7,  7,  9,  3)  ■  ETA1*IC7*ZZ(101-57*(C11*ZZI11)  ♦  Sll*C10*C«n  * 


1ZETA1*(-C7*(C11*ZZ  111  )  ♦  51i*Ci0*C9)-S7*ZZH0)) 

PTI  7,  7,10,  31  ■  ET41*(-C7*C10*S6*S9-57*ZZI31  II  ♦  ZETAl *1 -(7*22 () 
11  )  ♦  S7*C 10*56*56  ) 

PT(  7,  7,11,  31  •  -tTAl*S7*(Cll*ClO*59-5U*ZZ  (131  )-ZET41*C7*(C11*C 
110*59-511*22(13)  ) 

PT(  7,  9,  9,  3)  .  xn»(Cn*C10*S9  ♦  511*22(16))  ♦  E  T  A I  •  (  C7  *  (C 1 1  *ZZ 
l(16)-sn*C10»59)  ♦  S7*ZZ(16))  ♦  2  E  T  A 1  •  (  C  7  *ZZ  (  16  I -S7*(  C 11  •  2  Z  1 14  ) -51 
11*C10*39)  ) 

PM  7,  9,10,  3)  *  XU*2Z(iei  *  £TA]*(C7*Z2(20I-S7*C10*56»C9I  ♦  ZET 
lAl*(-C7»C10*i6*C9- 5  7*22(201 

PT(  r,  9.11,  3)  >  X  I 1«IC11*ZZ(11  )  ♦  511*C10*C9)  ♦  E T A  1 *C 7* ( Cll •€ 10 
l*C9-5  U*2Z  (  1  1  1  )-2ETil*S7*  (Cl  1  9CKJ*C9-5n  *ZZ  (1 1  )  ) 

PM  7,10,10,  31  •  M1«ZZ(37)  ♦  £  T  A  1  *  (  £  7  »Z  Z  (  39  )  ♦  57*56*55*510  )  ♦  I 
ieTAi*(C7*S6*5v*SlC-57*ZZ(j9)) 


PT(  7,10,11,  3)  *  XU*ZZ{3n  ♦  :TA1*C7*ZZ(30)-2ETA1*S7*ZZ(30) 

PM  7,11,11,  31  ■  Xl  l*(Cn*C10*59-511*Z2(13))  ♦  £  TAl  *C  7*  ( -C 1 1  *22  1 1 
r3)-sii*cio*S9)-ZETAi‘S7*(-cu6ZZ(n)-5ii*cio*59i 
PTI  0,  1,12,  II  -  C7*(C11*ZZ(12)-S11*C10*C9)  ♦  57*22(9) 

PM  8,  6,12,  1)  «  Xl*sn*22(9)  ♦  3(1  )*(C7*Cll*22(9l  ♦  $7*ZZ(UI)  ♦ 
1  ZFTA*(C7*ZZII1  )-S7*Cn*ZZ(9)  ) 

PTI  0,  7,12,  II  •  C{n*(C7*ZZ(9)-S7*ICll*22(121-5n*C10*C9))  ♦  ZET 
l4*l-t7*(Cn*2Z(12)-511*ClC*C9)-S7*ZZ(9)) 

PM  0,  9.12,  1)  ■  X  1  *(-CU  *00*59  ♦  511*22(13)1  ♦  01 1  1  *(£7*  (Cl  1  *22 
1(13)  ♦  511*00*59)  ♦  57*22  (  1  5  ))  *  Z  E  T  A*  (  £  7  *ZZ  ( 1 5  1 -57  •  I C 1 1  *22  ( 13 )  ♦ 
1  51 1 *C 10*59  I  I 

PM  0,10,12,  II  •  XMZZI17)  ♦  £lll)*(C7*ZZ(19)  ♦  S7*C10*56*C9)  *  2E 
lTA*(C7*O0*56*£9-57*ZZil9)) 

PT(  3,11,12,  1)  •  XI  •(Ol*ZZ(12l-51l*O0*C9)  ♦  0(  1  )*C7*  ( -Cl  l*O0*C 
19-S11>ZZ(  12  ))-ZeTA*S7*(-CU*O0*C9-5U  *22(121) 

PT(  8,  1,12,  2)  •  £7*22(25  1-57*00*56 

PT(  e,  6,12,  2)  •  -XI*511*C10*56  ♦  0(  1  )•( -C7*Ol*O0*56-57*Cl0*C6» 
I  ♦  ZFT  A®( -C7*O0*C6  ♦  57  *  01*00*56) 

PM  a.  7,12,  2)  -  C)  1  )*<-C7*O0*S6-57*Zl(25)  )  ♦  Z  ETA*  ( -C  7*22  1 2$ )  ♦ 
1  57*00*56) 

PT(  0,10,12,  2)  «  XI*ZZ(21)  ♦  Q(I )*(C7*Z2(23)  ♦  57*56*510)  ♦  ZETA* 
I  (£7*56*510-57*22  (23) ) 

PM  8,11,12,  2)  *  XI*ZZ(25)  ♦  Q ( 1 )  *C 7*22  (27  )-2ETA*S7*ZZ (27) 

PT(  6,  1,12.  3)  *  £7*101*22(13)  ♦  511*00*59)  ♦  57*ZZ(15) 

PM  8,  6,12,  3)  *  X1*5U*2Z(15)  ♦  C  <  I )  *  ( C  7  •£  1 1  *22  ( 1 5  1  *  S7*2ZI14)» 

1  +  ZETA*(£7*ZZ116)-S7*01*2Z(151) 

PT{  0.  7,12,  3)  -  Cll  )»(£7»ZZ(15)-S7*(01»Z2(13)  ♦  S11*C10*591)  ♦ 
UETA*{-C7*(£n*ZZ(l3)  ♦  SU*CiO*59J-S7*ZZ115)  ) 

PM  8,  9,12,  3)  •  Xl*l-Ol*£10*£9  *  511*22(11)1  ♦  01 1  1  •(  C7*IC1 1  *22 

Mil)  *  Sn*£10*C9)  ♦  57*22(10)1  ♦  2E  TA*IC7*Zll  10  1 -57*  (C11*Z2  111  1  ♦ 

I  511*O0*C9)) 

PT(  8,10,12,  3)  .  XI*ZZ(29)  ♦  Q 1 1 1  •  (  C7*22  (3 1 1  -  57*00  •S4*S»  1  ♦  2CTA 
1*  (-C 7*0 0*56*59-57*22 (31  1  1 

PM  0,11,12,  31  •  Xl*(Cll*Z2(13)  ♦  S11*C10*S9)  ♦  0(1 1 •C7*(C11*C10* 
159-511*22(  131  )-2ETA*S7*(  Cl  1*00*59-51 1  *22  (13)  I 
PM  8,  1,  6.  II  •  £7*01*22(9)  *  57*22(111 

PTI  8,  1,  7,  1)  •  C7*2Z(9)-S7*(C11*2Z(12)-511*C10*C91 

PTI  8,  1.  9,  1)  •  C7*(C11*ZZ(I3)  *  S11*C10*S9I  *  $T*ZZ(191 

P7(  8,  1,10,  1)  •  C7*ZZ(191  ♦  S7*C10*S4*C9 


OTI  ri  1  11.  I)  •  t  ?♦  I -:t  1 1  ot  10  •C9-S1 1  Z  C 12  >  t 

PT(  d!  fc!  b,  I)  >  Ki*sn*22«nj  ♦  oil )*<c7*cii*iztm  ♦  s7*2nion 

I  ♦  zMftp(C7ozinci-s7*cii*22nin 
PTI  s,  b,  7,  n  -  tmMC7*ZllUl-S7»CU*ZZJ9»)  *  ZETAM-C7*cn*ll 

-  Xl.sn*zz(l5)  ♦  0(UMC7.C11*ZZ(15I  ♦  S7*2Zn*l» 

1  ♦  zi  T  A  =  ( c  7*/z  u  t n  »zz  1 1 5  n 

*PM  a,  6,10,  It  «  xMSn«>C10*S6*C9  ♦  on  J*K7*C11*C10*S6*C9  ♦  S7K 
110<»1.6*C9»  ♦  Zf  T4nC7«C10e-t6»C9-i7*Cn  •C1C*S6*C9) 

PH  a,  6,11,  IJ  ■=  XI  ni  1  “ZZ  (  9J-CU  t*C7*Sl  MZZ  19)  ♦  2ETA*S7*SJI*2Z1 

*pTi  8,  7,  7,  n  •  cn  )*i-c7*nii*zznz)-sii*cio*c9)-S7*zzi»n  *  ze 

ltAn-C70ZZ19l  ♦  S7MC1HZZ  nz)~sii*cio*C9))  * 

PM  0,  7,  9.  1)  •  CinMC?»22n5)-S7.(Cn*ZZn3l  ♦  S11*C10*S9I»  ♦ 

1  7 P T 4 ®  1 -f 7 * { C 1 1 ® Z Z ( I  3 t  *  S 1 1 *0 10*S9 ) - S7*Z Z ( 1 5)  I 
PTI  9,  7,10,  n  •  cn  )» 1C7»C10»S6*C9-S7*ZZ  C 19)  )  ♦  ZeTA*l-C7*ZZn9» 

^PH’el^niirn  «  -on  i«s7*i-cn®cio*c9-sii*2ziut  )-ZETA.C7M-cn* 

*PT*(*0!  9y9!Ml^«*xi®(-Cl  1»C10*C9  ♦  Sll-ZZlU))  ♦  0(1 )  •1C7*KU*ZI 

ini)  ♦  sii«cio*c9i  ♦  s79Zznon  ♦  zeta»(c 7*Z2 no)-57*(cii*iiiii i  * 

^Pirs^^JIlol’l  >  .  X1»ZZ(29)  ♦  on  )*(C7*ZZ(3l)-S7*C10*56*S»»  ♦  ZEIA 

'M('8r5;nrn''Milcn®zztn)  ®  5n.c10.s9)  .  on )*c7.kii.cxo. 

I S9-S 1 1 ®ZZ ( 1 3 ) I -Z t T A.S7. ( Cl  I .010*59-51 1 ‘Z Z (13) )  ,cf* 

PT)  a.lO.K,  1)  =  Xl.ZZ(33)  ♦  on).(C7.ZZ(  35  )-57.S6.C9.S10l  ♦  ZETJ 
PTI  a.n.io,  II  -  Xl.ZZ(33)  ♦  Q(l).{C7*zn35l-57.Sfc*C9.Sl0»  *  ZETA 

‘;);S:;Snn^ir^/nnnUt  .  on,.c7.zzni.-zETA.s7.zz,i.. 

PT(  8,11,11,  It  ®  *i«c-cii.cio.c9-5n.zzn2))  .  on).c7.i-cii*zin 


sii 'cio«c9) -zer4.57*(-cii*zzn2)  ♦  5n.cio.C9) 


PTI  8,  1,  b,  2)  »  -C 7 *0 1 1 "C I  0 • 56- 57 •€ 10 *06 
PTI  8,  1.  7,  2)  *  -C7.C10«S6-S7.2Z  1251 
PTI  8,  1,10,  2)  '  C7«ZZ(23)  ♦  57*56.510 

PTI  0!  ll'll  21  •  -XHS11.C10*C6  ♦  On).(-C7*Cn*C10.C6  ♦  57*010.5 
‘^^^:^:rn^^^‘':^'^l;™:sKr^^7.cn.c.o.56, .  zEtA..c7.cii 
511*56*510  ♦  iin)*(C7. 011*56*510  ♦  S7*C6*510» 
0.i,*c7.5n.cio*56-ZEiA*S7.sn* 
7,  7,  2)  »  on  |.|-C7.ZZI2S)  ♦  57*010*561  «  ZE 1  A*  I  C7*C10*S6 
2)  «  CI11.IC7*56.510-S7*ZZ(23))  ♦  2E  T  A*  I -C7.ZZ  U  J 1-57 

1  956‘'  sit'  I 

PTI  8,  7,11,  2)  .  -on  )«S7.ZZI27)-ZET»*C7*Z2(27) 

PTI  i.!lCno!  2)  .  X1.ZZ(  2  7  )  .  01  1  1  *  107*22  1  26)  *  $7.C10*S6I  ♦  lETA* 


-xl*cn  *010.56  ♦  01  l)*C7.Sll.C10*S6-ZE7A*S7*5n* 
on  )  .  I -C7.ZZ  I  2S  )  ♦  57.010*561  «  ZE 1  A* I C7*C10*56 
C I  I )  . I C 7 *56*510 - S7*Z Z 1 23 ) )  ♦  2ET A* I -C7*ZZ liJ 1-57 
-on  ).$7.ZZI27)-ZET»*C7*Z2(27) 

X1.ZZI27)  ♦  01 1 )*IC7*ZZ 126)  ♦  $7*010*561  ♦  lETA* 


*prrHllo!in^2r!^x!*ZZl23)  ♦  on  )*C7*Z2i22)-ZETA*S7*ZZ  (Zij 
Pl|  a, 11,11,  2)  »  XI*ZZ(27)  ♦  ,n I *C7*ZZI 26)-ZETA*S7*2Zi2*l 

(Ml  8,  1,  6,  n  «  07*0  1  1  .22  (  15  )  *  57*21  1  16) 

B,  j,  ],  3)  -  C7.ZZ(15)-57*ICll*ZZn3l  ♦  S11*CI0*S9) 

PTI  8,  1,  9,  J)  =  t7.icn*zzni)  ♦  511*010*09)  ♦  s7*zzno) 

PTI  6,  1,10,  3)  »  t7*2Z  131  1-57. 010*56*59 

VA  8;  l'Ai\  3I  n*iniIz(iir^^lin)McUcn*2znA)  ♦  S7.2zii6n 

.i.sii-i/dn  •  ei,i«n?»cii-:i<io,  . 

.  a(ii..-CT.cii.cio.s,.j,-5TC 


(Ml  8  ,  1  ,  6 

PTI  »,  1  ,  1 

PTI  8,  1 ,  9 
PTI  6 ,  1,10 
PTI  8  .  1.11 
PTI  8  ,  6,  6 


110.(6  591  ♦  ZrTi  1-07.0  10.06*  59  *  57*tli -tlU-Ob-oy 

PTI  t,  6.11,  5)  .  xJ*cn.ZZll5)-in)*C7*51l*ZZn5)  ♦  ZE7A*S7*5U*I 

'pTl^a,  7,  7,  3)  .  01  1  1.1-0  7*  1  0  1  1  *22  1  1  3)  *  51  1  *010 *59  ) -57*ZZ II S )!  ♦ 

1  i-I-O?  ZZt  IS)  ♦  S7.(Cll*ZZn3)  ♦  511*010*59))  _ 

^TI  8.  7,  9,  81  ^  01 1  i.(C7^z2no)-S7*icii.zzni)  ♦  5n*cio*c9M  * 
WE(a  M-O7.l0ll»zznn  ♦  Sll*C10«09)-S7*2ZnO)  ) 

PTI  8,  7.16,  3S  .  cn  I  *1 -07»CIO»56. 59-57*22131  ) )  ♦  ZE T A* (-07*221 31 

*PTI  A  nin  31  -  -on  )*S7. 101  1  *010*59-511  *22(13)  )-ZETA*C7*(Cll*Cl 

.  X1.|C11*C10*59  ♦  511*zzn4))  ♦  on  )*(C7*(Cll*Z2( 
lUl  Mi’OiO.sn  '  S/.ZZIUn  »  2ETA*(C7.ZZn6>-S7*ICU*2Zll6)-Sll* 

*PT(  n'i.l-.  3)  =  XI*2Zn8)  ♦  Oni*(C7*2ZI20)-57*C10*56*C9)  *  ZETA 

^pT(^i  -‘r-^n.icuizziii)  ♦  5n*cio*c9i  ♦  on  )*C7*(CU»cio* 

ir^-sii  zzci  n-7rTflts7ncii*cio*c9-sii*22nn» 

P((  ".rMM  n  =  X1*22I37)  ♦  01  1  1  *107*22(39)  57*56*59*510 )  ♦  2E 

1 1 M  7  ■ '1  S-; n-s7=-z2 1 39  n 

,11,  )\  --  xi*2ZI3n  *  on  )*C7*2ZI30)-2ETA*S7*IZ(30I 
PM  -.li,!!;  7  xincii.cio*S9-su*2zn3n  *  on)*C7*(-cii*z2n3 

1  l-Sl  ,  .r  :.  VM -2“:  14^57*1-0  11  *22  n3)-511*C10*S9) 


F-9 


S0B8JUT1NE  O'cft^ 

THIS  SUBROUTINE  LEFINES  THE  PSRTUL  DERIVATIVES  PERTAINJNC  TO  fMI 
ftNOULAR  VELUCniFS,  POTENTIAL  ENERGY,  DISSIPATIVE  ENERGY,  AND 
GENERALIZED  FURCES. 

COMMON  /[>ERlVl/0(ll>,t0111»,000(lI» 

COMMON  /i tRIV^/Al  ,a2.A3,ASTAR,KRYl,XKY2 

C-iHMlJN  /LERI  V3/41SUB,A2SU0,A3SUB,XKK1  ,XKK2,AlfA« 

COMMON  /l.ER  U<./XL  12 .2  I  .XH  (2  ,2  )  ,XN(  ?,2I  ,XM2,2  I 

C.MM JN  /J£RlvS/El.E2.E3,FFl,FF2,Ff3,01,C2,03,Kll,ETAl,2EIAl 
CDMHilN  /0ERIVt/81,B2,e3,BlBAR,Xl,ZET4 

COMMON  /OER  lv7/PG«fe  ,12,12 *3  ) 

COMMON  /L£RlVE/XII,cfl,ZETAB,X|R,ER,2ETAfc,.XlC,EC»ZETAC 
COMMON  /GERiV9/dLFT,CCFT,RDFT ,FOfG 

COMMON  /NAME  I/PT  (8  ,11  ,12  ,3  1  ,PM(  5,12 ,12 ,31  ,PU(  A,U  ,U,GI 
1  ,PD  (3  .1  I  •  12  I 

CCjMMON  /NAME7/CA(l,l,l,4),0P(l,l,ltNl 

1  ALPHA2 , ALPH A3 ,Hhl ,HH2 ,Hm3 
COMMON  /TRl&/C6,Sfe,C7,S7,Ce,S8,C9,S9,C10,S10,Cll,SU 

derivative  FUNCTIONS  N/R  TO  R.E.  lANGULAR  T£»NS  ONLY) 


),12,  11  -  CDI9I-CI0 
>.12,  2)  ■  -00(9)RC10*S6  ♦  ODaO>*C6 
1,12,  21  •  -00(9 1»S10*C6 
>.12,  31  ‘  -gD(9l*C10»C6-00ll0)*S6 
),12,  31  -  COl9)*SlO*Sfe 
),12.  11  •  CO(9)*C10 

1.12,  2)  •  -0D(91*S6*C10  ♦  00(10>*CB 

3.12,  21  •  -0D(9)*CB*S10 

3.12,  J)  «  -aC(91*C8*C10-ODI10)*SB 
Ph(  2,10,12,  31  •  CD(9l»S8*S10 

PW(  3,  b,12,  11  *  Sll*(-aC(9MClO*5t.  ♦  ODdOJRC*) 

PM(  3,10,12,  II  ■  -0D(9 1*511  *S10»Cfe  ♦  0C|9)*C10*CU 
PMl  3,11,12,  11  •  -Sll*lOC(<.)  ♦  OD(91*S101  ♦  CU*(80(91*C10*C»  ♦  Q 
10(101*561 

pg(  3,  6,12,  21  ■  Cll  *( -ClC(91  •C10*S6  ♦  001101*Cfcl 
Py(  3,10.12,  21  -  -00(91*Sll*C10-Q0(9l*S10*Cll*C6 

PW(  3,11.12.  21  »  -SU*(QC(9)*tlO*C6  *  OD  ( 10  1  •  56  I -C 1 1  •  ( 00  ( G »  ♦  0D( 
191*510) 

Pm{  3,  6.12,  31  ■  -QD(91*C10*C6“00<10)*S6 
Py(  3,  10,12,  31  •  00(91*510*56 

PW(  4,  6,12,  11  ■  5) 1 *('QC (9  )*C10*56  ♦  CD(101*C61 

PW<  4,10,12,  1)  *  -0D(9)*511*S10*C6  ♦  0D(9)*C10*Cll 

P„(  4,11,12,  II  •  -51l*(0C(6l  ♦  00(91*5101  ♦  Cll  *(00  (91  •C10K6  ♦  0 

*p1(°4,^6!i2,  21  •  CIl*C7*l-0D(f 1*C10*56  ♦  00(I01*C61  ♦  S7*l“Q0(fl* 

7a2°*2J^!  C7*(0D(11)-0D(9I*C10*S6  •  00  ( 10  » *04 ) -S7*( -SI  1  •<  0 
1l(6I  ♦  00(91*5101  ♦  Cll*{C0(9l*C10*C6  ♦  0D(10>*S6I» 

C,12,  21  -  C0(9I*S10*S6*S7  ♦  C 7* ( -001 9 1 *5 U •ClO-OOl 9 » •S10*C 

1.12,  2)  ■  C7*(-511»(00<9)*C10*C6  ♦  00 ( 10 1 *56 1 -Cll • ( Q0( 6 »  ♦ 


PW(  1,10.12.  11 
PN(  1  .  6.12,  21 
PH(  1 .10.12.  21 
Ph(  1 ,  6.12,  31 
PM(  1,10.12,  31 
PW(  2,10,12.  It 


QD(10)*C61  ♦  C7*(-0D(9) 


PW(  4,10,12,  21  -  00(91*510*56*57  ♦  C 7* I 9 l -a n -ttu-wui ▼ r - j*w-v 
'  PVf  (^^!  1 1 , 12  ,  2)  ■  C7*(-511»(00<9)*C10*C6  ♦  00  ( 10  1  *56  1 -Cll  •  (  Q0(  6  »  ♦ 

3)  •  -Cl  1  *S7*(-Q0(9l*C10*56  ♦  QD(lO)*C6l  ♦  C7*(-0D(9) 

1  »C10*C6-a[)(  10  1  *561 

PrtI  4,  7,12,  31  •  -C7»{-511*100(6)  ♦  OD(9)*S10)  ♦  Cll * (QD(9>*C10*C 
16  ♦  JCIKI^SEI  l-S7»(C0m  l-OD(9)*tlO*56  ♦  QD(10I*C61 
PP(  4,10,12,  31  *  CD(91*510*56*C7-S7*(-CD(9)*Sll*Cl0-00(91*S10*Cll 

^PwV4,ll,12,  31  •  'S7*(-5U*(0D(9)  *C10*C6  ♦  QD  ( 10  1  *56  1 -C 11  ♦  ( 00  <*  1 
1*  Q0(9)»S10) 1 


py(  1 .)?,12,  1 1 
Pr.(  1  .12,12,  21 
PW<  1,12.12,  31 
py (  2  , 12 , 1 2 ,  II 
PMt  2,12,12,  21 
P21  2 ,12,12,  3  1 
P  2 (  3.12.12,  11 
1(91*5101 
P  »  I  3  , 1 2 , 1 2 .  2  1 
10(101*561 
PW(  3,12,12,  31 
P-(  4,12,12,  1  I 


' S7* ( -5 U * ( 00(9 1 *C10*C6  ♦  QD ( 10 1 *56 )-Cll*(00<*l 

LD<6l  ♦  00(91*510 
GD<9)*C10*C6  ♦  00(10>*S6 
-QD (9 1 ‘C 10*56  ♦  00(10I*C6 
C0(81  *  00(91*510 
;CD(91*Ce*C10  ♦  00(101*59 
-CD (9  1  •S8*C10  *  00(10)*C8 

SI  1 •( 00(91 •ClO*Cb  ♦  00(101*561  ♦  Cll*lODI6)  ♦  QD 
-511*100(61  ♦  00(9>*S10)  ♦  C11*(0D(91*C10*C6  ♦  0 
00(11  l-CD(91*C10*S6  ♦  QD(10I«C6 

CD(71  ♦  SI  1  * ( OD ( 9  1  •C10*C6  ♦  00(101*561  ♦  Cll*IOD 


(  4,12,12,  21  *  C7*(-S11*(QD(61  ♦  00(91*5101  *  C 1 1 • ( OD 1 9 1 •C10*C6 
au(l0)*56)l  ♦  S7*(0D(ll  l-OD(91*C10*S6  ♦  Q0(101*C6I 
i(  4,12,12,  31  -  C7*((.D(  11  l-wD(91*C10*S6  ♦  00  ( 1 0  I  *06  1 -57*  ( -5  1 1  •  1 0 
6)  ♦  UD(9l«510)  ♦  Cll •( CD(9l*C10*Cb  ♦  00{10)*561) 


7 

Pi,  ( 

1,12,  6, 

1  1 

1 

PWl 

1.12,  9, 

1  1 

SIO 

PH{ 

1,12,  9, 

2) 

C10*C6 

- 

PW  ( 

1.12,  9, 

3  » 

-C10*S6 

1 

py  ( 

1  ,12,10, 

2  1 

56 

PH( 

1  ,12,10, 

3  1 

C6 

\ 

Pw  ( 

2,12,  6. 

1  1 

1 

Py  ( 

2,12,  9, 

1  1 

510 

PW( 

2,12.  9, 

2) 

C8*C10 

py  ( 

2,12,  9, 

31 

-5«*C10 

7 

Pm( 

2,12,10, 

2  1 

59 

py  ( 

2,  12,10, 

3) 

C8 

py  ( 

3,12,  6, 

1  ) 

Cll 

py  1 

3,12,  6, 

2) 

-511 

p«  ( 

3,12,  9, 

1  1 

51 1*C10*C6 

♦  S10*C11 

p„  ( 

3,12,  9, 

21 

-511*510  ♦ 

Ci0*Cll*C6 

py| 

3,12  ,  9, 

3  1 

-CIO* 56 

py( 

3,12,10, 

1  1 

511*56 

< 

PW  ( 

3,12.10, 

2  1 

Cll*56 

PM  ( 

3,12.10, 

31 

C6 

7 

PW  ( 

3,12.11. 

31 

1 

py( 

4,12,  6, 

1  1 

Cll 

PW  ( 

4,12,  6, 

21 

-SU*C7 

py( 

4.12.  6, 

31 

511*57 

F-IO 


OKS 

0146 

0147 

014« 


OK*) 
01  SO 
OlSl 
0152 

01  S3 

0154 

0155 

0156 

0157 

0150 

0159 

0160 

0X61 

0162 

0163 

0164 

0165 

0166 

0167 

Ol6r 


0169 


C  WUX.ifUY.WUZ  «  wtx»wty«htz 

c 

OU  2  J>6,12 
00  2  X»6,12 

DO  2  L-1  .3 

2  PMIS.J.X.L  i«PW<4.J,K,LI 
C 

C  DERIVATIVE  FDNCTIOSS  W/R  10  POTENTUl  EhERGY  U2 

C 

POI  2,  4,12,  1»  »  XK(l,n»ZZI60J*S9*S10 

PUI  2,  5,12,  It  ■  XKI  1  ,1  I #ZZ 160) *C9 

PUI  2,  6,12,  1)  *  XK I  1 , 1 ) *22 160) •( XN n , 1 ) *22 ( 14 1  ♦  XH( 1  , 1  I *22  1 1 5  I) 

PUl  2,  9,12,  1)  »  XKU  ,1  )*22ItO)*I0I4)*C9*SlO-0(5)*S9-*l*C9*C10  ♦ 

ixNd.u^zzdoi-iiss^-xLn  ,n*c9*cio  ♦  xKn,ii»z2nn  ♦  *2*C9<»5io» 
PU(  ^,10,12,  1)  »  XK  ( 1  , 1  )  1 22  (60 ) * ( Q 14 > *59*C10  ♦  *1 •59»S IC'XN ( 1 , 1 ) • 
1$9*C1C*S6  ♦  XL ( 1 ,1 1*59»S1C  ♦  XM(  1  ,  1) •S9*C 10*C6  ♦  A2*S9*C10) 

PUl  2,  4,12,  2)  «  X«  (1  ,2»#ZZ(61 l*S9*S10 

PUl  2,  5,12,  2)  =  *X(  1  ,2 ) *22 (61 ) *C9 

PUl  /,  6,12,  2)  •  XR  (1  ,2  i  *ZZ  161 )*(XNI 1 ,2)*Z2( 14)  ♦  XH ( 1  , 2 ) *22 ( 1 5 ) I 

PUl  2,  9,12,  2)  ■>  XM1,2)*ZZ(61)*|CI4)*C9'»S10-C(5)*59-AI*C9*C10  ♦ 

1  XN( 1 ,2 ) OZZ ( 10  ) -XL ( 1  ,2  )*t9*C10  ♦  Xm 1 , 2 ) *2 2 ( 1 1  ) - A3 *59  ♦  A2*C9*S10J 
PUl  2,10,12,  2)  '  XK I  1 ,2  I »ZZ <61 ) *  I  014 )*S9*C 10  ♦  41 *59*51 0-XHl 1 ,2 ) ♦ 
1S9*C10*S6  ♦  X L 1 1  ,2 ) *59*51  D  ♦  XM ( 1 , 2  )  *  59 *C  10»C 6  ♦  A2*S9»C10) 

PUl  2,  4,12,  31  •  XM2,1  )*22  162  )*59*510 

PUl  2,  5.12,  3)  •  XK (2 ,1) «22I62)*C9 

PUl  2,  6,12,  3)  *  Xh  (2  , 1  )  *22  (  ^2  )  *1  xm  2 , 1 )  *22  1 14  )  ♦  XH(2 ,1 1  *22  1 15  ) ) 

PUl  2,  9,12,  3)  *  XX (2,1  I *22 (62 »* (XNI2,n *22 ( 10)-XLI2 ,1  )  •C9*C10  ♦ 

IXM 12 ,l )*22 ( 1  n  ♦  CI4)9C9*51C-3(5)»59-A1*C9*C10-A3*S9  ♦  A2»C9*S10) 
PU(  2,10,12,  3)  -  XM2,1  )  *22  (62  )*{ -XN(2,1  )*59*C10*56  ♦  XII2, 11*59* 
1510  ♦  XK(2,1  )w.59*ClD*C6  ♦  C(4I*59*C10  ♦  A1*S9*S10  ♦  A2*59*C10) 

FU(  2,  4,12,  4)  .  XK|2,2)*22I63)*S9*S10 

PU(  2,  6,12,  4)  .  XK(2,2I *22(63)*C9 

PU(  2,  6,12,  4)  r  XK(2.2)«22tb3)*(XN(2,2)*2Z(14)  *  XH 1 2  ,2 ) *22 < 1 5 1 ) 

PU(  2,  9,12,  4)  *  XK(2  ,2  ) *22  (63  ) *( XNI2 ,i ) ‘ZZl 10 )  *  XH(2,2>«22<11)- 

1  XL  (2,2 )*C9*C10  ♦  0(4)*C9*510-0(5I*59-41«C9»C10-A3*S9  ♦  A2*C9*SJ0) 
PUl  2.10,12,  4)  .  XK  :2,2) *22  (63  )*(-xN(2,2)*S9*C10*56  ♦  XN(2,2)*S9* 
^  1C10*C6  ♦  XL (2  ,2 )  *59*510  ♦  C|4)»59*C10  ♦  A1*S9*510  ♦  A2*59»C10) 

C  DERIVATIVE  FUNCTIOKS  V/i  TO  POTENTIAL  ENERGY  U3  PARI  ONE 

PUl  3.  2,12,  1)  •  •2Z(53)*(-22  I45)*C6*2  ♦  22 1 46 )  *51*2  I • ,S 


F-11 


.  1  1G 

.  j  n 

:  ii 

:in 


'  1  75 
:i  7fc 
:j  77 

cue 

0179 

0180 

0161 

0162 

0U3 


0184 

01s5 
01  b6 

0187 


PO(  1.  ‘..w.  n  *  -ZZ(53I*(U(45>*C8‘>2-1Z  146U5««2)*.S 

pji  5,u.  n  *  -zz  u  i»  »<zz(45)  •»  zz(46»*ce*2i*.5 

PUI  j,  fe,l2,  n  =  -ZZ153l''UiK5)»(ZZ<4n*S«  -  CB*ZZ(42I»* 
u.-*  ZZ«4-,|  MZZ<41>‘>C8tS8*ZZ(42)>»2,>*.5 
PU<  3,  8,12,  IJ  ■  -ZZ<53)MZZ<45)MZU43»»C8-ZZ1441*S0J*2  ♦  ZIU6I 
Ut-ZZ(43l*i8-ZZ(44)*C8»*2)*.5 

r 

f  0EPlV*TtV6  PUNCT!Ct*S  ^(/«  TO  POTENTIAL  EKEBCY  U3  PART  TNO 

C 

Put  3,  2.12.  2>  •  'ZZ(59I  M-ZZt45»*C8*2  ♦  ZZ 146  J  »•  .5 

PUI  J.  4,12,  2>  •  -ZZ  159U  (ZZ  C45)*Ce*2-ZZ(461  »S8*2»*,5 
Put  i,  5,12,  2)  «  -ZZ(59>MZZ(45>*S8P2  ♦  Z  Z  1 46  1 ‘C  8*2  1  *.5 
Put  3,  6.12,  2)  »  -Z21S9I  *(Z2C45)*1ZZI4U»S8  -  Ce*IZt42»)p 
)2.*  ZZt46)  MZ2  (4n»Ce«Sfi*2Z  142  )  >»2,  l*.5 
PUt  3,  0,12.  21  •  -ZZ1591  MZZt45)*UZ143MC#-ZZ(44l*S»l*2  ♦  2Z1461 
l»l-ZZt43I»S6-ZZ144»*C0»*2»«.5 
C 

C  OERIVATIWE  FUNCTKNS  N/R  TO  POTENTIAL  ENERGY  04  PART  ONE 

C 

PUt  4,  2,12,  n  -  XKYl *Z2 I64»el-Cfl  ♦  U 
1  -012  »«  II  .-C6  I -< -1  ALPhAl -C  t  3  U ‘>C9*S10  •»  ALPHA2*C10  ♦  AlPMA3*59pS10 
1  -A2-G2-Ut2n‘'C8-t  lALPHAl-Ct  31  )*S9-ALPHA3*C9-l3-G3>*Si  ♦ALPHR2 
1  >A2-C2 

POl  4,  3,12,  11  -  )tliYl*ZZI641*lC6*C9*S10-SepS91 
1  -Ot2)»C9»51O*C0  ♦  0(2)«S9oS8 

Put  4,  0,12,  11  »  AAYl'>ZZI641»tCi4t-A3-C3  ♦  ALPHA3*C9  ♦  S»*!ALPHA1 

I- U13I n-S09(-G2-GI2J~A2  ♦  ALPhA2*C10  ♦  ALPHA3pS9*S10-C9»SI0*1  ALPHA 

II- 0(31111 

1 *01 2 1 »l - ( alpha  1 -C( 31 1 *C9»S10  ♦  ALPHA2*C10  ♦  ALPHA3*S9*S10-A2-4a 
I  -0(2 1 1*58-QI2 )» ( ( ALPmA1-C(3  1  1 *S9  ♦  AL PMA 3*C9-A3-G3 I PCi 
Put  4.  9,12,  1)  «  xky1»ZZ(64)=>1C0p|ALPHA3*C9*SIO  ♦  S9*S10*  1  AlPHAl- 
IC(3II1  ♦  S6»(-alPhA3«S9  ♦  C9'>1  ALPHAl-COm  1 

I  -3< 2 1 »1 ( ALPHAl-Ct 3  1  I  •S9»S10  ♦  ALPHA3pC9*S10 I *CI-0(2 I •! ULPHAI - 
I  0I3))*C9-ALPHA3*S91»SB 

PUt  4,10,12,  n  «  A4Y1  «ZZ ( 64 1 *€8 *1 •ALPHA2*S10  ♦  ALPHA34S9*Cie*C9*C 
llO*  t ALPHAl-OlS  1 )  I 

I  -Q(21*t -l6LPHai-C(3npC9*ClO-ALPHA2*SlO  +  AL  PHA3  *S9*C10  !•€• 

C 

C  DERIVATIVE  FUNCTICNS  N/R  TO  POTENTIAL  U4  PART  TMO 

C 

PUt  4,  2,12,  2)  •  *KY2*ZZ (65  IPI-CS  *  II 
1  -0(2l»(l.-C8)-t-l-AL PH A1 -0(311  *C9#510*ALPHA2*C10*ALPHA3*S*PS10 
1  -A2-G2-Q  (  2 1 IOC8 - ( I -ALPHA  1 -01 3 1 1 *59- ALPhA3*C9-A3-G3l*SI  ♦  ALPMA2 
I  -A2-02 

PUI  4,  3,12,  21  «  **tY2*ZU651*Ka*C9*S10-S8»59l 
I  -0(2  1  »C9'»S1G®C0  ♦  0121'>S9»S@ 

Put  4,  6,12,  21  »  XKY2'»ZZ  (65}®(C8*(-G3-A3  ♦  ALPMA3»C9  *  $9*I-ALPHA 
n-Q(  1  I  1  l*S8'>(-C2-C(2}-A2  ♦  ALPMA2<»tlO  ♦  AlPMA3*S9*S10-C9*S10*( -ALP 
IHAI -0<  3  1  I  I  1 

1  *012  I  -  1  -  I -ALPHAl -31  3  1  m9»S10  ♦  ALPHA2*(10  ♦  ALPHA3*S9*$10-A2-C2 
I  -0(2  )  l‘‘Sa-U2  I  M  (-ALPhA1-C13H*S9  ♦  ALPHAI*C9-A5-63  ••Cl 
PUI  4,  9,12,  21  •  *KY2's^ZZ  (651o(C8AtaLPHA3»C9*5l0  ♦  S9*S  10*  I -ALPMAl 
1-0131)1  ♦  ie*! -aL'7KA3*S9  ♦  C9*(-ALPHAl-0(311l  1 


1  -0(21®(  (-ALPHAl-L(3)l»S9»510  ♦  ALPHA3*t9*SI0  1*C8-0(21*( <-AlP«Al- 
1  013)  1  'C9-AlPHA3'S»)»Sa 

PUI  '-.((.IP,  2)  =  *KY2»ZZ  <651'»C8s-(-ALPMA2*S10  ♦  ALPHA1*$9*C10-C9PC 
lUiM -ALPhAI-,,(  3)  11 

1  -0(2  1  «-( -1 -ALPt-Al-Ol  3  1  UC9«-C10-ALPHA2*S10  ♦  ALPMA3*59*C10»*C8 
CcfllVATlvE  FUHCTiCSS  W/R  TO  DISSIPATIVE  ENERGY  UBARl  PRRT  ONE 


0109 

POl 

1.  4.12, 

1  1 

.  S9'>S10 

019D 

POl 

Ir  5,12, 

1  1 

.  C9 

0191 

PD  1 

t  ,  8,U, 

1  1 

X  XMl  ,1  l»ZZI14l  ♦  XHII  ,1  1*ZZ(15I 

0192 

PC  1 

1  ,  9,12, 

I  1 

X  AN  ( 1  , 1  )  UZ  1 101-XL  n  ,l  1  *C9*C10  ♦  XMlltl  1*ZZ(11  1  ♦ 

1  A2>C9^S10-A3«'S9 

♦  il(4)-C9^'5l0-a(5}*S9-Al*C9*CI0 

0193 

POl 

1,10,12, 

1  1 

*  -XN( 1 , 1 US9»C 10*S6  ♦  XL  1 1 ,1  I«S99S10  ♦  XN(1,1)*S9 

1  9C10 

S'CN  ♦  t2*S90 

CIC  ♦  0(4)o59*C10  ♦  A1*S9*S10 

0194 

POl 

1  ,  4,12, 

2  I 

X  S9«510 

0195 

PC  1 

1.  5,12, 

2  1 

X  C9 

0U6 

POl 

1.  6,12. 

2  1 

*  XMll ,2|xZZ(  151  *  XN( 1  ,2}»ZZ 114) 

0197 

PD  1 

1  ,  9,12  , 

i  1 

X  -XL  (1  .2  )H9»C10  ♦  XK(l,2)oZZtll)  ♦  XN(  1  .2  )  •ZZ  1 10 

1  1  ♦ 

A2  ^Crs-Sl  j-A3 

XS9  ♦  0(4 1 »C9sS10-Cl5 1 •S9-A1»C9*C10 

0198 

PC  1 

1  ,  IC.U, 

2  1 

X  XL  (  1 ,2  1  '  S9*S10  ♦  XMl  1  ,2  1  *S9»C10*C6-XNn,Z»*S9*Cl 

lO'Se 

♦  A^.'S9*Cl': 

♦  .(4)xS9xci0  ♦  A1»S9»S10 

0199 

PO  ( 

1  ,  4,12, 

31 

X  59»S10 

0200 

PDl 

1  ,  5,12, 

3  1 

X  C9 

0201 

PD  ( 

1 .  6,12, 

3  1 

X  XH  (2  ,1  1  -ZZ  115  1  ♦  XN12  ,1  UZZ  ( 14) 

0202 

P3( 

1,  9.12, 

3  1 

=  -XL(2,l)n9K10  ♦  XM(2,1  I.ZZlll  I  *  XN(Z,n*ZZ(10 

1  1  * 

A  2  •  C  y  »  S  1 C  -  4  3 

-59  ♦  C<4  |»C9’51&-0(5)*S9-A1*C9'»C10 

02J3 

PD( 

!  .r  ,12. 

3  1 

X  Xi  (2 ,1  1  «59»  SIO  ♦  XN(2,n»S9*C10*C4-XN12,ll*S9*CI 

*  :2>S9 

1C 

♦  UI4US9-C10  ♦  A1#'S9*S10 

0204 

PDl 

1 .  4,12, 

4  1 

X  S9«S10 

0205 

POl 

1 ,  5,12, 

4  1 

X  C9 

0206 

PDl 

1  ,  6.12, 

4  1 

--  **1(2,2  1  »ZZ(15)  ♦  XNI2,2>*ZZ(14) 

0207 

p:  1 

1  ,  .  u , 

4  1 

X  -XL  (2  ,2  1 'C9‘*C10  ♦  X11(2,2I*ZZ(11  1  ♦  XNIZ.Z  )  ^ZZ  ( 10 

1 1  ♦ 

-  2  ‘  9  1 

0-A3 

:9  ♦  C(4 |x:9'S10-C(51*59-A1*C9*C10 

0208 

pj( 

1 .  r:. .  1 2 , 

4  1 

X  *l  12  .2  1 -SVX-SIO  ♦  XN12,2  )*S9*CI0*C6-XNI2,2)*S»*C1 

C 

c 

1 

*  .'.2  'S9^C  l  , 

»  ;(4)»,,9xcn  ♦  ai  959*sio 

,  ‘  I  I  k'E  FufiCT 

i:ss  m/R  id  DISSIPATIVE  ENERGY  UBAR2  PART  ONE 

0209 

c 

P3{ 

2  .  b  .12  , 

1  1 

X  ZZ  (52)M-2Z(77)*C€*2  ♦  2Z(78»*S8*2I 

0210 

PC  ( 

,  s  ,  1  2  , 

1  1 

X  ZZl  -  2IHZZ(  77l*C6*2-ZZ(78)*Si*2) 

0211 

PCI 

/■  ,  5.12, 

1  1 

X  7Z  1  82  )  X  1  ZZ  1  77  ) ->5892  ♦  ZZ(7B)*C8*2I 

0212 

POl 

2,  6,12, 

1  1 

X  ZZU2)- tZZ(77)MZZ(79)AS0  ♦  Z2(B0)*C»)*2  ♦  IZITi 

1  UIZ 

/  ( ?  9 1  «■ :  0 

cousai^i) 

0213 

PCI 

2  ,  5!  ,  1  2  , 

1 1 

X  ZZ  132 1 » (ZZl 77) *(ZZ( 76)*Ce-ZZ(751*SII*2  ♦  2ll7il* 

1  l-ZZ  (76)  ’SR-ZZ  175  1  ='C8)92  1 

OEfii.ATIvE  FUNCTICNS  N/R  TO  DISSIPATIVE  ENERGY  UBAR2  PRIT  TMO 


0214 

P:i  1 

2  . 

2,12. 

21 

0215 

POl 

2  , 

4,12, 

2  1 

0216 

PD  1 

X  , 

5,12. 

2  I 

021  7 

*  01 

2  , 

6.12, 

2  1 

ZZUl  1«(-ZZ(77)*C8*2  ♦  ZZ(78I*SI*2I 
zzifl |9|ZZ(77)»C8*2-ZZI7R)*SI*ai 
ZZ1P11»1ZZ1771*S8«2  ♦  ZZ(78>*CR*2I 
ZZ181)9(ZZ(77)*UZ(79)*S8  ♦  ZZIIOIKR)*!  ♦  IIITt 


F-12 


0^16 


0^l9 

0^i0 

QZil 

0222 

0223 


0224 

0225 

0226 

0227 

0226 


0229 

0230 

0231 

02  32 

0233 

0234 

0235 

0236 

0237 

0236 


023V 

0240 

0241 

0242 


0243 


0244 

0245 


0246 

0247 

0248 


0249 
02  50 
0251 


0252 
02  53 
02  54 


0255 

0256 

0257 


025S 

0259 

0260 


0261 

0262 

0263 

0264 


C 

c 

c 


c 

c 

c 


( 

c 

c 

c 


c 

c 

c 


c 

c 

c 

c 


c 

c 

c 


c 

c 

c 


c 

( 

c 


c 

c 

c 


ll»J2M791*C6-/ZieC»»Sp)»2) 

o:.t  2,  8,12,  21  »  ZZtai  )*1Z2(  771  •tzz  (  76)*C8-12«75»*SH42  ♦  IZOII* 

I  C -m  76)  "SB-iZ  (  751  *CB)»2J 

CcBlVATlve  FUNCTIONS  ta/R  TO  OISSIFiTiVE  EKERfcY  UBR*)  P**!  ONE 

PD)  ),  2.W,  1»  •  ZZ(7n  •tZZ)70)*S»*2-Z2169»»C8*2» 

P,)  3,  4,12,  n  •  ZZ  (  711  •  «-Z2(  70)*S8*2  ZZ1691*C8*2I 

PD(  3,  5,12,  1)  »  ZZni)*WZ<7CMC8*2  ♦  ZZ<69I*S8»21 
PDl  3,  6,12,  1)  =  ZZ(7n*lZZ17Q)M-22i72)‘S8  4  ZH73I*CI»*2  ♦  2216 
1y»*(:Z)72)^-C6  ♦  ZZ(73)<>Se  1  •Zl 

PUl  ),  8,12,  11  •  22(71 (*)ZZ)701*I-ZZ«67) •S8-ZZ<661*CI1*2  ♦  221691 
1*I  22  (  67)^'C8-Z2(66)«S0)»2I 

derivative  functions  V/R  to  dissipative  energy  U6AR3  >A*T  TNO 

PCX  3,  2,12,  2)  '  2Z(74)»(Z217C)*S8'»2'22169I*C8*2> 

PD(  3,  4,12,  2)  *  ZZ174)4l-2ZI70l*Sa*2  ♦  Z21691*C8*2I 

POl  3,  5,12,  21  •  ZZ<741*(ZI(7C)»C8»2  ♦  2ZI69J*S6*2) 

PD(  3,  6.12,  21  •  22(74)4(£Z170)*1-22I72)*S8  *  ZZI731*C8I*2  ♦  2216 

lV)»(Z2t72)»C8  ♦  22(73  1*58  l-Z) 

PDl  3,  8,12.  2)  ‘  22(741*(Z2(70)*(-Zl(67)*S8-22«66)*Ctl*2  ♦  22(691 
MUZiaTl^CS-ZZlbEMSSI'Z) 

DERIVATIVE  FUNCTIONS  V/R  TO  CENERAIUEO  FORCES  PART  Of  VECTON 

AB.ee.ce 

PG(  1,  1,12,  II  •  ZZ185)*C7  ♦  ZII9»*S7 

PCI  1,  6.12,  II  *  Z£TAeM-Z2(91*Cll»S7  ♦  ZZ(11)*C71  ♦  XIB*22(9I*S1 
11  *  (EB  4  iJI  1  )  )*(2Z(9  )»C1 1*C7  ♦  22(U)«S7) 

PCI  I,  7,12,  II  «  Z=TAB*(-Z2  Ce5l*C7-22(9l*S71  ♦  (E6  ♦  0m)4(-22(6 
15»*S7  ♦  2Z(9)*C7) 

PC(  1.  9,12,  1)  •  ZE1ABM-22  (901*S7  ♦  22U5)»C7»  *  R1 B* ( 22 (13 1  *511 
l-S9*ClO*Cin  ♦  Ite  ♦  IXl  l  )»(Z2(90I*C7  ♦  ZZ(15I*S7I 
PG(  1,10,12,  1)  «  2ETAB4I-2Z1P610S7  ♦  C9*C 1 0*S6*C7 1  ♦  Rlfi*(-C9*C10 
1  osii  «C6-C9*S10»Cin  ♦  (EB  ■*  0(  n  )  *(22  (86)  •C7  ♦  C9*C10*S6*S71 
P&(  1,11,12,  1)  •  -2£TAe*S7*1-2Z(12)''Sll-C9*C10*Clll  ♦  XIB*ZZ(691 
1*  C7<>(Ee  ♦  C(  1  )  )  *  I -221  12)  •SU-C9»C10*Cin 
PCI  1.  1.12,  2)  •  22I87)’»C7-C10«S6*S7 

POl  1,  6,12,  2)  •  2E I AB*(C10*Cll*St*S7-C10*C6»C71-XlB*ClO*Sll*S6  ♦ 
1  IE8  ♦  (X  1  )  )  *( -Ci0»Cl  1  *S6H7-C10*C6*S7» 

P0(  I,  7.12,  2)  •  ZfcTAB*! -22 (e71»C7  ♦  C10»S6*S71  ♦  (EB  ♦  fl(l»l*(-I 
lZ(e7)*S7-ClO»S6*t7l 

POl  1,10,12,  2)  •  7ET AB*I S10*S6*C7-S7*(-C10*S11-S10*CI1»C6) 1  ♦  Rlt 
i»(C10«'Cll-S10‘-Sll*C6)  ♦  (EB  *  Oil)  1»(S10*S6*57  ♦  C7*(  “C10*S11-510* 
1  Cl  1 *0  6  1  I 

POl  1,11,12,  2)  -  -ZETA9*ZZ (88 ) *57  ♦  XIE»Z2(87)  ♦  22 ( BB I *07* I EB  ♦ 

10(11) 

PCI  1.  1,12.  31  -  221901*07  ♦  221151*57 

POl  I,  6,12,  3)  •  ZH4P*(22n4l*C7-22ll5)»Cll»57J  ♦  XI6*22(15»*S11 
I  ♦  (EB  ♦  0( I >  )  *(22 ( 14)*S7  ♦  22I15)*C11*C7) 

POl  1.  7,12,  3)  «  2EIA8*I-22(901*C7-Z2I15I*S7»  ♦  lEB  ♦  Om»*(-l21 
1901*57  ♦  22I15)*C7) 

PCI  1,  9,12,  3)  •  2ETiB*IZ2llO>*C7-5T*(22ni)*Cll  ♦  C9*C10*SUH  ♦ 

I  X 1 B*l 22 (1 1 1 *51 1 -C99C  10*C1 1  )  ♦  lEB  ♦  Jl 1 > I  *  1 2 2  110 1 *57  ♦  C7*(ZZ(111 

l*C 11  ♦  C9*C10*Sl  1 )  >  ,  _ _ 

PCI  1,10,12,  3)  «  2tTAB*l-ZZI39l'>S7-S9*ClC*S6*C7»  ♦  XI  B«  (S9*C10*S1 

ll«C6  ♦  59*SlC;*tU)  ♦  lEB  *  0  U  )  )*l  Z2  I  69  1  *C  7-59  *C10*S6*571 
PCI  1,11,12,  3)  »  -ZET  AB*S7»  (-2211 3)«S11  ♦  59*ClO*Clll  ■»  Xlfl*Z2(90 

II  *  C7*(£e  ♦  Oil  n  *(-Z2ll3>»511  ♦  S9«C10*Cin 

derivative  function  m/r  tc  generalized  forces  part  of  vector 
AR.br, CR 

PGI  2,  1,12,  II  »  221851*07  ♦  22(91*57 
PCI  2,  1,12.  21  •  22 (87}*C7-C1G*S6*S7 
PCI  2,  1,12,  31  *  Z2I90)*C7  ♦  22(151*57 

derivative  functions  W/R  to  (iENERALlZEO  FORCES  PART  Of  VECTOR 
AC.BC  ,CC 

PCI  3,  1,12,  1)  *  22(05>*C7  ♦  22(91*57 
PCI  3,  1,12,  2)  *  ZZ (e7)*C7-CI0*S6*S7 

PC(  3,  1,12,  3)  -  Z2I90)*C7  *  2Zli5)*S7 

VECTOR  BUFTA,  BOFTB,  BOFTC  in  GENERALIZED  FORCES 

PCI  4,12,12,  1)  -  -BUFT*(ZZ(85)*C7  ♦  Z2(9)*S7) 

PG(  4,12.12,  2)  *  -BOFT*( ZZ(87>*C7-C10*56*57) 

POl  4,12,12,  3)  »  -BOFT*( ZZ (90)*C7  ♦  Z2115)»57l 

VECTOR  ROFTA,  ROFTB,  ROFTC  IN  GENERALIZED  FORCES 

PCI  5,12,12,  n  »  RGFT*IZZie5)*C7  ♦  2Z(9)*S7» 

P0(  5,12,12,  2)  •  R0FT*(ZZ(87)*C7-C10*St*S7) 

PCI  5,12,12,  3)  »  R0FT*(ZZ(90)»C7  ♦  ZZ(15)«571 

VECTOR  CUFT A  ,CCFTB  iCCFTC  IN  OENERALIZEO  FORCES 

PCI  6,12,12,  1)  -  C0FTXZZ(65) 'C7  ♦  ZZI9)*S7I 
PGI  6,12,12,  21  =  CUFT-  (ZZ<67)-’C7-C10*S6*S7» 

PCI  6,12,12,  31  *  COFT  ■  (  Z  Z  I  90  )  *..7  ♦  ZZ(15)'>i7) 

jELTA  IJ  in  U2  of  potential  energy 


;a(  1 ,  1 ,  1 ,  II  '  -xNi 1 ,1 )-4  3 
1  ZZ  (  1  5)  -XL  l  l  ,1)  ssvt'f.io 
_ Al  1 ,  1  ,  1 ,  2  )  *  -XNI 1  ,2 1-43 
lZm5)-XLIl,2)*S9*C10 
,1)  J  ,  1  ,  1  ,  3  1  ^  -XM2,l  J-4  3 
1(13)  *  XN(2,1  )*ZZ(15I 
r:4|  1,  1,  1,  4)  «  -XN(2  ,2  I-A3 
1(13)  ♦  XNI2,2)*ZZ(15I 


ZZ(b3)  ♦  R)<(l,l)*Zili3)  ♦  xtxil, 
ZZI63)  *  XNIl  .2»*Z2(13)  ♦  XA(ll,2l* 
ZZ  )!.' 3  ) -XU  2 , 1  1  *S9eC10  ♦  XNI2,l)*2Z 
22I83)-XL(2,21»S9*C10  ♦  XKIZ, 21*12 


DELT4  IJ  PRIME  IN  GENERALIZED  FDRCfS 


0265 


DPI  1 ,  1 ,  1 ,  1) 


-XLll.l) 


CI3)-A1  *  22(84)  4  XMll. 11*22(121  ♦  X 


026b 

0267 

026S 

0249 
02  70 


1  \(  1  .  W  *..2  <  9»  ♦  *Hl,l>*C9*C10 

.P(  1,  I.  1,  2)  «  'XL(1.2>  0O)>«}  *  IZ(S4I  *  XLIl,2t*C«*C10  «  X 

l-^n  ,2  ►•ZZC12  I  ♦  XN(1,2»»2Z«9> 

1.  J,  1.  3)  •  -XLI2,1)  ♦  6(3)-*l  ♦  ZZ(B4)  •  XL(2 . 1)  ^Cf  *€10  ♦  X 
1-12,1  »»ZZn2l  ♦  XM2,1»*ZZ191 

>.P<  1,  1,  1.  41  ■  -XL(2,2>  4^  0«3)-Al  *  ZZ(B4I  ♦  XL  12 .2 1  •C9*C10  *  X 
1'(I2,2)‘ZZ112  »  •*  XN(2,2)*ZZ(9» 

SETU^ll 

ihO 


0001 
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0007 
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0009 
0010 
0011 
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0013 
0014 
0015 
0014 
0017 
0018 
0019 
0020 
0021 
0022 
0023 
0024 
0025 
0024 
0027 
0028 
0029 
00  30 
0031 
00  32 
0033 
0034 
0035 
0034 

0037 


0038 

0039 

0040 

0041 
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C 
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SUBXOUTINE  SOLVf 

this  SU663UTIhE  DECOUPLES  ^HE  XtCTOft  DlPEEXEXTUL  EOUATICM 
**QD'.,-8  ANC  SLLWtS  fCR  THE  OCOIll  WHERE  1-1, 2,. ...11 
X  HO.IFIEO  SCUXRE  ROOT  HEThQO  IS  USED.  THIS  METHOD  TAKES  ADVANTAGE 
..f  THE  SYMMcTRV  C,  F  THE  MASS  MATRIX. 


OIME^tSlON  $1  10,1  1)  .XKlin  ,CI11  > 
COMMON  /M&HE4/ AIll ,11 ) .Fllll 
common  /OERlVt/Ollll,QOIll),X(lll 


SI  1  ,1  1 
Sll  ,21 
Sll ,31 
Sll  .4  I 
Sll , 51 
Sll .61 
Sll, 71 
SI  1  ,8  I 
Sll, 91 


Al 1  ,1) 
All, 21 
All, 3) 
All. 41 
All. 51 
All  ,4) 
All  .71 
All  .81 
All ,9) 


Sll ,10  1  «  All .  101 
Sll. Ill  >  All. II) 


Sll,n**(-1  1 

•  A(2.3)-C( l>*Sll,21*SIl,3l 

•  A{2,4)-CI1I*SI1,2)»SI1.4) 

>  A(2  ,5  )-CC  1 )  *S  II  ,2  l•SM  ,5) 

-  A(2,4  l-r.  1 1  I  •Sll  ,2)*S(I  *61 
«  AI2  .7  l-CIl  )»SI1 ,71»SI1,2) 

•  A(2,ei-cn  1*51 1 ,6  )*sii  .2) 

-  Al2,9)-cni*sil,9)*sil,2) 

A|2.10)-Cll)*S{l,10J*Sil,2l 

SI2.11)  •  A|2,U)-Clll*Sll,n)*Sll,2l 
CI2I  •  UI2,21-C  m*Sll  ,2)**2)**l-l  1 

AI3,4)-Ci2l*SI2.3l*S(2,4)-Cll)*Sll.3l«Sll,4} 
A(3,5)‘t(2)*5l2,5)«SI2,3)-Lin*Sll,3)*$ll,5l 
AI3,6>-Cf2)*SI2,4)*SI2,3l'Cn)«Sn.3l*Sll.4) 

A  I  3,7) -Cl  2  MS  I  2, 7) -SI  2,31 -cm  *5  (1,7)  *51 1,31 
A(3,6)-C(2)*SI2,S)»SI2,3}-Cmoin,8l*S(l,3l 
AI3,9J-C(2)''S(2,9)»SI2,3)-cm*S(l,91*$ll,3l 
S(3,10)  «  A(3,10)-SI2,10)*Cl2J*S(2,3»-Cm*Sll,10»*Sll,3l 
S(3,ll)  •  A  I  3,1 1  I-SI2 ,11  MC  12  )  *5  (2  ,3»-Cm  *S(  1  ,11  1*51 1 ,3  1 
cm  «  IAI3,3)-CI2)*SI2,3)»*2-CI1  1*51I,3)**2)*»I-1  I 
5  14,5  1  •  AI4,5)-C(3M5(3,5l*SI3.4  1-C(2)*S(2,S)*$l2,4)-Cm*S(l,4)* 

1511.5) 

S (4,4) 

1511.6) 

SI4,7) 

15(1,41 
514,6) 

ISIl  ,41 
5(4,91 
15(1.4} 

5(4, 10)  *  A(4,10>-S(3,10)*C(3)»S(3,4)-S(2,lO)'»C(2l*SI2,4)-Cm*Sll 
1.10)*S(1,4) 


cm  - 

512.3) 

512.4) 
5(2.5)  ' 

512.4)  < 
S(2,7) 
5(2,6) 
512,91  I 
$  12,10  ) 


5(3.41 
5(3,5) 
5(3,4) 
513,71 
513,6  > 
513,91 


A(4,4)-C(3)»S<3,4l*S(3,41-CI2)*Si2,4)*5(2.4l-Cll)»Sll,4)* 
A(4,7)-C(3)*S(3,7)*S(3,4)-C(2)*SI2,7)*5(2,4)-Cm*Sn  ,71* 
A(4,e)-C(3»*S(3,8l*Si3,4l-C(2)*S(2,8)*SI2,4)-Cm*Sll,8)* 
A(4,9)-CI3)*S(3,9)*SI3,4)-C(2  )*SI2,9)»SI2,4)-Cm*S(l,9)* 


! 


F-14 


I 


00*.? 

0043 

0044 

0045 

0046 

0047 

0048 

0049 

0050 

0051 

0052 

0053 

0054 

0055 

0056 

0057 

0058 

0059 

0060 

0061 

0062 

0063 

0064 

0065 

0066 

0067 

0066 

00  69 

00  70 

0071 

00  72 
0073 
00  74 
00  75 
00  76 

00  77 

0078 

0079 

00*^0 

00  31 

ao82 

0083 

0034 

0035 

0036 

0087 

OO-iB 

0089 


s«4,n)  •  4i4,ii)-5i3.in*ci3)*s<3.4i*su,in*c«2>»si2t4j-'C»ii*sii 
1 ,  in  >s  ( 1 .41 

CJi.)  .  (  A  i4,4>-c  t3  »*SO,4  »**2*ct2»*$<2,4l  ••2-Ctl  1  ‘Sn  ,4  ►••21**l-l 

SJ5,6>  •  AI5,6)-C  14)»S<4,5»*SI4,6»-C«3)<>5<3,5»*$13,6»-C<2>*$«2,»»* 
1S(2,6)-C(  1  ,5MS(1,6» 

515,71  .  »15,7I-C{4)*5<4,71*S<4,51-C(3)*SI3,5I*SI3,7»-CI2>*S12,5J* 
lS(2,7)-Cin»S(l,7MS(1.5> 

SI  5,3  1  «  *(5,8P-C14)*S14,e>*5l4,5>-CI3)*S0.51*$l3,8»-CI21*S(2,5)* 
l5(2.ai'C(l  |0Sn  ,8>«5(1  ,5) 

St  5, 9)  M  a(5,9  nc l4l*SI4,9)*SI4,5>’C(3l*5(3,5)*SI3,9l-C(2)*S(2tSI* 
lS(2,91-C(n*Sll,91*Stl,5l 

S(5,10>  -  At  5, 10  I -C  (4  I  »S(  4, 10)  *5 14, 5  1*5  (3, 101 -CO  I 'SO, 5  1-5  (2, 101* 
lC(2»»St2,5)-tin*Sn,10)*Stl,5) 

SIS, 111  '  A(5,n»-Ct41«St4,m*514,5)-SI3,in*CI3)*S13,5)-S12,ll)9 
lCI2»='S{2,5)-Cm‘»Sll,in*SCl,5) 

CIS)  «  I A(3,5) -C 14 )*St4,5 l•*2-C^3^•S^  3,5) ••2-CI2)*5(2,5)**2-Cll»*S 
)  II  ,5»<-»2  I**!-!  ) 

5(6,7)  •  A(6,71-CC4)*5t4,7)*5l4,6)-5t5,6l*S15,7)*C(5)-C(3)*S13,*»* 
lS(3,7)-Ct2)‘5(2,t)*5l2,7)-Ctl)*5(l,7)*Sll,6) 

S(6,5)  «  AI6,B)-CI4)*SI4,e)»St4,6)-5(5,6)*Sl5»l)*C(5)-CI31*S(3,*)* 
lS(3,a)-Ct2)‘‘5l2,b)«S(2,6)-Cll)*5(l,«)*S(l,6l 

S16,9)  -  AI6,9)-CI41*S(4,9)*5I4,6)-S(5,6)*5IS,9)*CI5)-CI3)*S13,4)* 
lS13,9l-Ct2)»SI2,6)»Sl2,9)-Ctl)*5(l,9)*S(l,6) 

St6,l0)  •  A(6,10)-C14)*S(4,lO)*5l4,6)-5l5t6)*515,lO)*C15)-$l3,IOI* 
1C13)*SI3,6)-S(2,10  l■>CI^)•5(2,6)-C(l)*S(l  ,10)*SI1,*» 

5{6,lll  •  A(6,n)-Ct4)*S(4,ll)»S(4,6)-S(5,6)*5l5,ll)*CI5)-SI3,in* 
lCI3)»S(3,6)-S(2,n)»CI2)*5(2,6)-C(l)*Sll,n)*S(l,6) 

C  16)  «  I  a(6,6)-C(4  )*SI4,6)**2-C(3  l-St  3,6)  ••2-Ct2)*S(2,6l**2-Ca  )*S 
ltl.6)*«-2-SI5,6)»‘2‘'CI5))»»(-l  ) 

517.8)  -  ft) 7,0 )-CI4)*St4,7)*St4,8)-5(5,71*SI5,«)*C15)-C(6)*S(*,f»* 
lSlb,7)-CI3)*5(3,7)«5t3,8)-Cl2)*S(2,7)*SI2,B)-Cn)*S(l,7)*51l.ll 

517.9)  *  ft(7,9)-C(4)»SI4, 71*5(4,9 1-515, 7)*S(5. 9) *CI5)-C 161*514, 
1S(6,7)-CI3)*5I3,7)»S13,9)-C(2)*S(2,7)*S12,9)-C11)*511,7)*S11»*) 

5(7,10)  •  ft( 7,10I-C(4)*5(4,7)*S(4,10)-5I5,7)*SI5,10)*C15)-C(6»*SI6 

1.10) *S(6,7)-SI3,10)*C(3)*S(3,7)-SI2,10)*C(2)*5l2,7)-Cn)*SU,7»*$l 
11  ,10) 

5(7,111  •  Ar7,ll)-C(4)»Sl4,7)»S(4,ll)-S(5,7)*S(5,ll)*C15)-CI*l*SU 

1  .11  1*5(6,7  1-5(3,11  )*CI3)*513.7I-S(2,1 1  1*01 2)*SI2,7)-Cm*Sa,1l*$( 

n  .1 1) 

C(7)  •  (A|7,7)-C(4)*St4, 7)»*2-Cl6) *5(6,7) ••2-C13)*5(3,7)**2-C12l*S 
It2,7)**2-C(l)*5n,7)**2-5(5,7)**2*CI5))**I-1  ) 

5)8,91  •  ft(e,9l-C(4l*5(4,8)*S(4,9)-Sl5,8)*S(5,9)*C15)-C17»*5l7,H* 
1S(7,9)-C(6»*5U,6)*516,9)-CI3)*S(3,»)*S(3,9I-C(2)*S(2,«)*S(2,»)-C1 
ll)*5(l,6)*S(l,9l 

'10,10)  «  A(P, 10) -C(4)*SI4,8)*SI4,l0)-5(5, 81*5(5, 10) *0(51-0171 *5 17 
1 ,«)*S(7,lO)-t (6l*5(6.8)*5(6,10)-St3,l0)*C(3l*S(3,8)-512,l0)*C(21*S 
1  12,8 )-C (1 ) *5  1 1 ,8  I »5t I  ,10) 

5(8,11)  •  A(8,n )-C (4)*S(4,8I*5(4,11 )-S(5,8)*S(5,ll )*C15)-C<7)*S(7 
l,«l*S(7,ll)-C16)f5(6,a)*5  16,n)-5(3,lll*C(3l*5(3,8)-S(2,lll*C(2)*S 
1  (2,8 )-C( 1 )*S  n  ,8  )«SI1,11 ) 

C(8)  ■  tAI8.8)-C (4)*S|4,8)**2-C(7)*SI7,6)**2-C(6)*S(4,I)**2-C(9»*S 

1I3,B)**2-C(2)*512,0)**2-Ctn  *5  (1  ,8  5,8 )**2*C 15) »**l-l  I 

5(9,10)  •  A|9,10)-C14)*S(4,9)*514,10)-CIII*5I6,9)*5(I,10I-CIT»*$17 
1  ,9  1*5(7,10  )-».  (6  )«  5(6,9  1*5(6,101-5(5, 9)*5(5,101*C(5)-S(3,10I*C(3)*5 
U3,9l-SI2,10)*C(2)*5(2,9)-C(l)*5ll,9)*SH,10) 

5(9,11)  ,  Ai9,n  )-',  (4  )*s(4,9)*s(4,in-C(8)*s(8,9)*5te,ii)-C(7)*S(7 

1,9)*,  17, 111-1(6) ‘5(6, 91*5(6, in-S(5,9)*5(5. Ill  *0  (5  )-5(3,lU*CO)*$ 
l(3,9)-S(2,ll)*C«2)*5(2,9)-Cll)*S(l.9)*S(l,ll) 

C(9)  *  ( 0 (9,9)-C (4  )  *  S 14  ,9  )**2-C (8  1*5 18,9) **2-C (71*51 7,t)**2-C(6»*i 
1 (6,9)**2-CI3 )*5(3,9)**2-C 12)*S(2,9)*»2-C 11 )*5(1,9I**2-5(5,»»**2*C( 
15)l**(-l  ) 

>(10, in  •  A(in,n)-C(4)*5(4, 101*5(4, lll-C(8)*S(8,10)*S(8,ll»-S(*, 
1  lO  >*S(9,n  )  *CI9)  -Cl  7  I  «SI  7 ,101*51  7,u  )-C(6)*S(6,10)*5(6,lll-S(5,I0» 
l*S(5,in*C(5)-S(3,10)«5(3.11)»C13)-5(2,10)*S12.11)*C12)-Cin*$(i,l 
10*5(1,11) 

CdO )  «  (A(10,10)-C(4)*S14,10)**2-C(6)*5(8,10)**2-C(7)*SI7,10I**2- 
1C(6)  5  (6,10) *«2-C in  *5(1 , 10)**2-S(9,10I ••2*C(9)-SI5,10 )**2*C15)-S( 

13.10) **2C(3l-5(2.10)»*2*C(2n**(-l  I 

con  *  (  A1  11  ,11  )-C(4)*S(4,in**2-C(10)*5ll0.111**2-0(B)*S(»,ll)** 
12-Cl 7) *5  t  7, 1 n**< -C (6)*5(t,  11  )**2-C(l l•5(l,lI)••^-5(9,ll)••2•C(♦^- 
lS(5,in**2*ClS)-5(3,in**2*CI3)-5(2,in**2*CI2>l**(-l  I 

xK(i  I  *  fin 

xK.(2i  *  f(2)-xK(  n  *C(n«sii,2) 

XK(3)  *  F  (  3) -XM  2  )*c  (2  )  *S(2,3)-XK1  n*C(l  1*5(1  ,3) 

XK(4 )  *  F|4)-XK(2 ) *C (2 ) *S(2  ,4)-XK(3  )*C 131*5(3 ,41-XKll 1*0(1 »*Sa ,4) 

XKIS)  *  F(S)-C(4)*XK(4) *5 (4,5)-XK(2)*C(2) •S(2,5)-XKI3)*C(3)*S13,5) 
1  -xx(  n  *C(  n*5(  1 .5) 

XK«6)  *  F(6)-C  14)»xk(4)*5(4,6I-S1S,6)*XK(5)*0(S)-KH12)*CI2»*S12,6) 

i-xK(n*ci3)»S(3,6)-xK(n*c<n*sii,6) 

XM7)  =  FI  n  -c  14  I*XK(4I*5(4,7)-S(5,7}*KX(5)“C(5)-XK(6I*C(6)*516,7) 
l-XM2)H(2)^5(2,7)-XK(31*r(3)*5t3,7l-XK(n*Cll)*5n.7) 
x^t8>  =  F(s)-CI4)*Xk(4)*514, 81-5(5, 81* Xft (51 *C (5)-XK<6 )*C (6) *5 1  6, 8) 
1-XK(7I*C(7)*S(7,6I-XK<2I*C(2)  *5 1 2 , 8  ) -XK ( 3 ) *0 1 3 1 *5 1 3 ,8 ) -XU  1 1 1 *0 1 1 » * 

)  >n  ,3  ) 

x-x(9)  =  F(9)-C(4)*XK(4)*5(4,9)-X»((6)»CU)*S16,9)-XK(7)*C(7)*5l7,f ) 

l-Xt'  {>,  )*C  (8  )*5(8,9)-XK(2)*C  (2  1*5(2, 9) -XK(3)*C13)*5I3,9)-XK1  5  1*515,9 
n  *c  ( 5)  -xr<  (  n  *0  in *5(  1 ,9) 

XM15)  =  F{ir)-C(4)*XK(4)*514,10)-XK|6)*C(61*5(6,101-XK(7)*C(7)*S1 
i;,iri-XK(B)*U8)*S(e,l0)-X)((9)*5t9,10)*C(9)-XK(2)*5(2  •  10  1*0(2  1-XKl 

I  ;  )*S  (3  ,1C  )  *(.  (  3  )-xM  5  nS(  5  ,  10  )»C(51-XKin*Clll*5(  1 ,10) 

x^(ln  =  F(ni-n4)*xK(4>*5(4,in-xK(6)»ci6)»5i6,i  i  i-xki  7)  *ci7»*si 
i7,ni'XK(8)n(rt)*5(8,in-xK(9)*s(9,in*C(9)-xx(ioi*ciio»*siio»ili- 
IW  (2  1  5(2,  111' C  (2  )-XK  (3  1*5(3,111*0  (31-XK(51*515,m*C(5)-XX(l)*C(I 

I I  'S(  1 , 1  n 

«.(in  •  xKdii  *c  (in 

<(10)  ‘  -xm  )  *c (10)  *5 1 10 ,11)  ♦  xK(io)*cuo) 

,(9)  =  -xni  l*S(9,ll  )*C(9)-X(10)*S19, 101*019)  ♦  XR(9)*C(9) 

<(-)  s  -X  (  1  n  ' '!(  H  *5  (  6  ,11  )-X  (10)  *C  (8)*5ie,10) -X19)*C  181*5(8,91  ♦  X 
me )  >c  (8  I 

X(7)  ,  -xn  n*C(  7)  *517,1  n-X(  101  *CI7)*S( 7,10 »-X(91*C171  *517,9 l-Xli 
n*C  (  7  )  -  S  (  7,8  1  ♦  XK  (  7)*c  (7  ) 

x(6)  -  -<(in*C(t)*5(6,in-X<10)*C  {61*5(6, 10)-X(9»  *0(61  *5 16, 91-XI9 

n^cibi-sib.Bi-xiTncibnsu,?)  ♦  xK(6)*(i6i 
X(5)  =  -5(5,6)*X(6)*C  15  )-S(5,7)*xni*C(5  1-S(S,»)*X(8)*Ct5l-Xim*S 
n5,in'-CI5)-X(l0)*S(5,10)*C(5)-XI9)*SI5,9)*ClS)  ♦  XK(SI*C15) 
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oosc 

00^1 

0092 

0093 

0094 

0095 


xn,i  «  -ct4j*xc5j<>s{4,5i*c<4j*x(m»su,in-ci4>*xno»*sj<i,io)-ci4 

|»X(9)»Sl‘*,9»-C<4)#X(«1«SI4,8J-C(4»*X(7J*S<4t7»-C14»*Xt6t*SI4,6)  ♦ 

,^*V**'‘!^i,,Cl3  1.SO.5»-H4l*CJ3)‘>S(3,4t-Xnn*5l3,ni*C(3l-l(ll0l 

ltiI'I’-M5»*c!2)M!l2!5l-MM'C(2l»S«2.4»-XI3)*CI2»*S<2.3I-K(ni*5 

t2.1l»»Ct2>-X(lOI«'S(2,10>*C<2)-X«9)*C«2>*Sl2t91-X(B)*C»2)*Si2t«l-* 

1 1 U  ( 1 , 2 )  -*  n  1 » K  <  1  >  *  s  j  1  ♦  1  n  -X 1 10 » »c  t  n  *s  J 1 , 1 0 }  -X  <  9  >  "C  j  1  j*  5  a  .9 ) 

( 8  j  * :  n )  9  s  ( 1 18 1 -X 1 7 )  sc  a  >  »sn  t7  i-x  (6 » *0 « 1 1  *5  c  1  lb )  ♦  xncn  cni 

xeTUHM 

END 
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SUBftOUTINE  XUTt* 

IHfS  SUBMOUTINE  INlEGfcbTES  THE  ELEVEN  SECCNO  D«0E9  NONLINEM 
OlFF-RrNTUL  ECOATIUHS  MHICM  OESCRIfE  ThE  NOTION  DF  THE  NllOAl 
A  fourth  order  RONCE-RUTTA  NETHOD  15  USED. 

CINENSIDN  QSAVEIll».CD5AVEni>iARCIli9> 

C'JKHON  /DATAl/TlKf*TlNEH,TIH€H2tTINEH9 

coNHjN  /DERI vi/a«m,ooan,ODDaii 

CONHHN  /DERI V9/BLf T  ,CCFT ,RQFT ,FOFC 
CnHN:N  /KUTTA2/llOfT,IROFT,l60FT 
CDI^MUN  /KUTT  A3/ IfiPT  S  ,  IRPT  S»  l&PT  $ 

CONHUN  /XREALl/BkChxaOSl  .BRCHVI 1051, ROOX  1105  l.RODVdOSU 

1  GANHAXaOSl  ,GiNPAYil05» 

DO  5  1-1,11 
OSAVE (11-0(11 
5  COSAVE  (  n-00{  1  1 

CALL  L INE  AR ( T IHt .0RCHX,BRCHV ,9GFT . tBPTS, I lOFT  I 
CALI  LlNEARdlNE  ,RCDX,ROOV  ,ROFT  ,  ItPT  S  ,  INOFT  1 

CALL^L INEAR(0(7I  ,0 ANNAX  ,&ANNAT  ,FOFt, KPT S , I60FT 1 

CALL  NANE 
CALL  SOLVE 
DO  1  L-l,ll 

AK(L . I  I -000 ( L  1  *T  INEM 

Q(L  1 -OSAVE  (L  1  ♦  TINIH2*00SAVE(L)  ♦  T  INEHB^ARU  ,  I  » 

1  CD4L 1-QDSAVE (L 1  ♦  AK(L,ll/2. 

CALL'LlNEARtTlIIe^BRCHX.eRCHV.BOFTaePTSaBOFTI 
CALL  LINEAR(TINE  ,RDDX,ROOV ,ROf T ,IiPTS,lROFT 1 

CALL^LlNF  AR(C  {  7  I  .OANNAX  ,GANNAy  ,FOF G ,  KPT  $  ,  KOFT  } 

CALL  NANE 
CALL  SOLVE 
DO  2  L-l.ll 

AK(L,2>*00D(LMT1HEH 

2  QD( L 1 -ODSAVE (L 1  ♦  AK(L,2)/2. 

CALl\inEAR(G(71  ,&ANNAX,6AMNAr,FOFG,KPTS,K0FTl 

CALL  NANE 

CALL  SOLVE 

DO  3  L-l,ll 

Ak (L . 3 1 -ODD (L  I  ♦T  (HEH 

U(L  1 'QSAVEIL  1  ♦  T  1NEH*0DSAVE  (L  I  ♦  T  INEH2*AI(  (L  ,31 

3  COIL 1-ODSAVE (L  1  ♦  AKiL,31 

CALL*LllAtAR(TIRe,6RCHX,BRCHY,eOFT,lBPTS,lBQFTl 
CALL  linear (TINE ,RCDX .ROD Y , ROf T  , IRPT5 , IROFT  1 

CALL^L 1NEAR(C(7)  ,G ANN AX , G AHH A Y ,F0FG , KPT S , I GOFT 1 
CALL  NAME 
CALL  SOLVE 
00  4  L-l.ll 


0048 

0049 

0050 

0051 

0052 


C(L  I -QSAVE  (L  KT  (MEM- (  OOSAVE  (  L  1  ♦  (AN(L,1I  ♦  ARIL.ZK*  AIHL*311/*.> 
4  tO(  L 1 -ODSAVE  (L  1  ■*  (AK(L»1»  ♦  2.-(AK(L»2)  ♦  AR{L,3))  ♦  AKIL»4ll/b. 
RETURN 
END 
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C  this  SUbROUTUE  CCHBINES  THfe  DIFFERENTIAL  EXPRESSIONS  IN  TmE 

C  PROPER  WAV  SUCH  that  it  YIELDS  THE  COEFFICIENTS  MATRIX  A  OF  THE 

,C  ACCELERATION  TERMS  TIND  THE  RIGHT  HAND  SIDES  QF  THE  DIFFERENTIAL 

C  EOUATIQNS. 

C 

DIMENSION  SUHM(5,3I 
'  DIMENSION  PKE(5.1I.12.3I 

COMMON  /DATA2/  XNN2 . XNN3 . r2 , 03 *6 ET AE .CKST . EL ANC .CU5 
COHMuN  /DATAi/flETAT 

CUMMUN  /CERi VI /Dill  I  •QDIl 1  I ,000(11  > 

CUHMUN  /DER I V  5/E  1. E2,  E3. F FI, FF2.FF3.D1.C2. 03.  X I  l.ETAltZETAl 
tOMMIlN  /OERI  Vb/Bl  .fi2,B3,BlBAR,Xl  ,2ETA 
COMMON  /DERIV7/Ptl6,12,12,3l 

CUMMUN  /OERIVe/KlB.EB.ZETAB.XIR.ER.ZETAR.XIC.ECtZETAC 
CliMHUN  /0ERIV9/6Lf  T.CUFT.ROFT  ,FOFC 

CLMHON  /NAMEl/PT(6.11tl2.3).Ph(5,12.12.3l,PU(6.11,12.6l 
1  ,PDO.n.l2.6) 

CUMMUN  /NAME2/IECS,IMASS 
CUHMUN  /NAME3/XMASS(5» 

COMMON  /NAMEA/AAAIll  .in.RHSIll  I 

COMMON  /NAME5/XUK(5I  »X |  Y Y  1 5  I  ,XI ZZ  I  5 1 , X IX V  1 5 1  ,X lYZ  1 5 1  .XI XZ  1 5 1 
.  COMMON  /NAMEb/CRAV.CBRCE.BETA 
COMMON  /NAME7/CAI l.l ,1.61 *DPI1  .1 ,1 ,6l 
common  /NANEG/CCIl »CC12,CC21 «CC22 
COMMON  /NAME9/M0I:R1 

c 

C  evaluate  PARTIAL  DERIVATIVES 

c 

CALL  DERFUC 
C 

C  DEFINE  FIRST  AND  SECOND  DERIVATIVE  TERMS  FOR  THE  KINETIC  ENERGY 
C  EXPRESSIONS  FOR  THE  FIVE  MASSES  INOT  INCLUDING  THE  ANGULAR  TERNSI 

C 

DO  2  J*l.n 
DO  2  K-1.12 
00  2  L«1.3 

PKEd.J.K.LI-PTd.J.K.L)  «  PT(2»J.K.L} 

PKE (2, J.K ,L I -PT ( 1. J.K.L)  *  PTI3.J.K.LI 
PKE(3.J.K,II-PKEII.J.K.L)  *  PTK.J.K.LI  *  PTIS.J.K.L) 
PkE(6,J,K,L)*PKE(1  .J.K.LI*PT(6.J,K»LHPT(6.J.K.L)*PT(7«J»K,L) 

2  PKE(5.JtK.LI-PKEIt.J.K.L)«PT(6,J.R,LdPT(b.JiK.L)«PTI0.J,K«LI 
C 

C  define  COEFFICIENTS  MATRIX  FOR  ODD  TERMS  RESULTING  FROM  KINETIC 

C  ENERGY  (DOES  NOT  INCLUDE  THE  ANGULAR  TERMS) 

C  ONLY  THE  UPPER  TRIANGULAR  TERMS  ARE  DEFINED  HERE 

C 

DO  5  U«I»IEOS 
DO  5  K>J.IEOS 
SUM.O. 

DU  6  1*1. INASS 
DO  4  1*1.3 


0035 

0036 


0037 

003B 


0039 

0040 

0041 

0042 

0043 

0044 

0045 

0046 

0047 

0048 


0049 

0050 

0051 

0052 

0053 


0054 


0055 


1)056 

0057 

005B 


0059 

0060 

006f 


4  SUM. SUM  ♦  KHAS5m*PKE(l,J,12.L)*PKEII.K.12.LI 

5  AAA(J,K)>SUM 
C 

C  INITIALIZE  RIGHT  HAND  SIDES  TO  ZERO 

C 

DU  6  I>1,1EQS 

6  RHSI I )*0. 

C 

C  KINETIC  ENERGY  EXPRESSIONS  FUR  TERM  3  (NOT  INCLUDING  ANGULAR  TERNS) 
C 

DO  8  I-1,1MASS 
OU  B  Ld.3 
SUMM(l.L)-0. 

DU  e  K>1,IEQS 
DO  6  J-1,1F0S 

8  SUMH( I ,L )*SUHM|  1  ,L  )  ♦  PKE ( I . J.K ,L I *QD( J) *QD(K  ) 

DO  9  J*I,IEQS 

OU  9  I «1  .IMASi 
DO  9  L*1 ,3 

9  RHS(J»*RHSIJ)  ♦  XMASS( n*PKE( 1 ,J,12,L)*SUMM|I  ,LI 
C 

C  DEFINE  THE  REST  (F  THE  COEFFICIENTS  MATRIX  FOR  ODD  TERMS  RESULTING 

C  FROM  THE  ANGULAR  VELOCITIES  IN  THE  KINETIC  ENERGY  EXPRESSIONS. 

C 

DO  11  J-l.lEOS 
DU  11  K>J,1ECS 
SUH-O. 

DO  10  1*1  ,  IMASS 

SUM. SUM  ♦  I0.5»X  IXXI  1  )«-PMI  I  ,  12  ,J,1  I  -  X  I  X  V  (  I  I  *PW  (  I  ,12  ,  J  ,2  t  )  • 

1  PU(  I  ,12,K,1 } 

2  ♦  (0.5*XIYYinePHll  ,12,J,2)  -  XI  YZll  1*PV(  1 ,12,  J,3n* 

3  Pi4(  I  ,12.K.2  I 

4  ■*  (0.5*xizzin‘’pw(i,i2,j,3i  -  xixzd  )*pw(i  ,i2.j,ni* 

5  PW( I  ,12,K,31 

10  SUM. SUM  ♦  (0.5«xuxm*pwd  ,12,K,1  1  -  X1XY(  n  «PW(  1 ,12.K.2>  I* 

1  PU(  I  ,12.  J.l 1 

2  ♦  (0.5oXlYVin‘’PWd  ,12,K,21  -  X1YZ(  noPWI  I  rl2,K,3d* 

3  PWH,12,J,21 

4  *  {0.5*XIZZ  d  )0PW  d  ,12  ,K  ,3)  -  XlXZ(n*Pkl  l,12,K,l))* 

5  PN( I  ,12, J,3I 

11  AAA( J.K) .AAAIJ ,K)  *  SUM 
C 

C  UBTAIN  ENTIRE  CUEFFIClfNTS  MATRIX 

C 

00  50  1.1,11 
DU  50  Jd  ,  1 1 
50  AAAIJ, d 'A^AI I , J1 
C 

c  RIGHT  HAND  SIDES  OF  ANGULAR  TERMS 

C 

OU  IT  J.l.lECS 
SUM.O. 

DO  II  K.l  ,  lECS 
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00b; 

00b3 


0064 


OObS 


OObb 


00b7 


0060 
0069 
00  70 
007} 


00  7i 


0073 


0076 
00  75 


00  76 
0077 


0070 

0079 


0000 
00  01 


0002 

0083 

0004 

0005 

0006 

0087 

0086 

0089 

0090 

0091 


0092 

0093 

0094 

0095 


0096 

0097 

0098 

0099 

0100 

0101 

0102 

0103 

0104 

0105 

0106 


Oil  I  )  I  *1  ,tH<SS 

SUH«50K  ♦  CO  .5*X1  KX  (  n*PW  M  ,6  ,  J.  I  )  *  X  I  XY  (  I  I  *06  1 1  ,K  i  J  ,2  I  I  • 

1  P6tl,12.12,l>*0C:<Xl 

2  *  |0.5*X|VY(M»PU(I  ,K.J.2)  -  XtVZin*PW(I.K»J,3n* 

3  PUC I .12.12t2>*OC(K) 

4  *  I  i.5*xuttn*pu<i .K,j,3i  -  xix2in*P6(i .R.j.m* 

5  PWC 1.12.t2,31»0C(KI 

SUH.5UH  *  (0.5*XiXXI  n»PWI  U12.J,1  >  -  XIXVI  l»*PW(  1 ,12,  JtU  1* 

1  PMC  I  ,K,12,U»CDtKt 

2  ♦  io.5*xiYYin*pmi  ,i2,j,2»  -  xiY2m*p«ii,i2,J,3n* 

3  PHC  t  ,K,12.2)*CD(K) 

4  *  (Q.5»XIZZ(n«PU<l,12,J,3>  <■  XIX1III*PMII.12.2»1H* 

5  PWC I ,K,12,3I«C0  CK  i 

SUH.SUH  «  10.5»XIXX<  n»PUlI.K.I2,n  -  XlXYIll*PIIII,K,12,2n* 

1  PWCI,t2,J.ll*C0CKl 

2  *  (0.5»XIYV<I  )«PWI1,K,12,2I  -  XIYZCU^PWCI, 11,12. Ill* 

3  PUI  1 ,12, J,2  )*CD(K1 

4  «  (0.5OXIZZ11  l«PUll,K,12.3)  -  X1XZII)*PU(1,K,12,1II* 

5  PM(I,12,J,3I*C0IK| 

5UM.SUH  *  10.5*XIXXCnoPKn.l2.12,ll  -  XIXYm‘Pl»n,12.12.2ll* 

1  PWC  l•K,J,ll•COI6l 

2  ♦  l0.5*X|YYCn*PU(1.12,12,2l  -  X 1 VI C 1 1  *PU  <  1  •  12 ,12 .3 11  * 

3  PWC  I  ,X,J,2I*CDCKI 

4  *  (0.5«XIZ2in*PU<  i  ,12,12,31  >  Xt  XZ  C 1 1  •PW(  I  ,  12,12 ,1 1 1  • 

5  put  I ,X,J,3I»G0CKI 

13  ftHStJ)«RHSIJI  *  SUH 

C  TfRHS  FDR  MINUS  PARTUL  K,E.  M/R  TO  THE  GENERALIZED  COORDINATES 
C  (ANGULAR  TERMS  ONLYI. 

C 

DU  693  J>1,IE0S 
SUH«9. 

DO  492  l>l,IMASS 

SUM«SUH  «  10.5«XIXX<1I*PWCI.J,12,1I  •  XIXV 1 1 1 *PU ( I , J,12,2I t* 

1  PW(l,12.12,n 

2  «  (0 .5*X I VV 1 1  I •PWC I  ,  J,12  ,2  I  -  X 1 YZ ( I )*PH ( I , J , 12 ,3 1 1 •R« 1 1 »12 ,12 ,2) 

3  ♦  C0,5«XI2Zin»PWIl  ,J,12,3)  -  X  IXZ  I  n*PU(  1 .  J,12,l  1 1  •PHI  I  ,12,12,31 

492  SUM. SUM  ♦  C0.5*XIXXC  nfPWCI,l2,12,n  -  XIXYm*PH(l,l2,l2,2n* 

1  Pwtl ,J,12.1  I 

2  ,C0.5*XlYV(  IMPbC  i  ,12,12,2>-XIVZ(I  )•PH(t,12,12,3M•PU(t•J,12»2) 

3  ♦  C0.5*X1ZZII  )*Pwn  ,12, 12, 31 -XI XZ  1 11  •PH  (1, 12, 12,  in  •?«(!,  J, 12, 3) 

493  RHSC J)-RHS( JJ  -  SUM 

C  RIGHT  HAND  SIDES  COMPLETE  AT  THIS  POINT  FOR  K.E .(TRANSLATION  AND  ROTATIONI 

C 

c 

C  RIGHT  HAND  SIDES  DUE  TO  POTENTIAL  ENERGY  OP  U1 

C 

DO  14  J*1 ,1E0S 

14  RHS( J|.RH$( J)  ♦  GRAV*(XHASS(1I*PKE(1,J, 12,31  ♦  XNA5S(2I* 

1  PKE{2,J,U,3I  ♦  XHASS(3I*PKE(3,J,12.3I  ♦  XNASSI 6 1 PPRE  16,2,12,31  ♦ 

2  XHASS(5)«PkE(5,J,12,3II 
C 

C  RIGHT  HAND  SIDES  DUE  TO  PCTENTIAL  ENERGY  OF -All 
C 

00  15  J*1,IEQS 

15  RHS(J).RHSIJnPUt2,J,12,l  l♦PU(^,J,12,2l♦PUI^,J,12,3l♦PU(^,J,12,61 

t 

c  RIGHT  HAND  SIDES  DUE  TO  PCTENTIAL  ENERGY  OF  U3 

C 

DU  16  J«1,IECS 

U  RHSC J)>RMS( Jl  •  PUC3,J, 12,11  •  PU(3,J,12,2l 

C  RIGHT  HAND  SIDES  DUE  TO  PCTENTIAL  ENERGY  OF  U6 

C 

DD  1 7  J«1 ,IECS 

17  RHS(J1*RH5CJI  ♦  PUC4,J,12,1I  ♦  PU(6,J,t2,2l 

C 

c  RMS  ^Ut  TO  POTENTIAL  ENERGY  OF  U5 

C 

S7«SINIQ(7n 
C7«CaSC0C7I I 

XLENGH  -  01«D1  ♦  {XNN2-CD2-02MC7  -  ( 03-03 1 *57 1••2  ♦  IXNN3  ♦ 

1  C02-02I‘'S7  -  C03-D3I*C7)**2 
XLENOH-SORT  CXLENGH) 

TERMl  «  F0FG»CXHR2«(-XNN3-(02-02)*57  ♦  (O3-03l*C71  -  XMN3*(-XNN2 
1  ♦  CI12-D2»®C7  ♦  ia3-D3n57n/XLENGM 
CAMST*0C7I  -  ELARG 
QBE  TAE»&Uf ToZEIAE 

IFC69fcTAE  .GT.  BETAEI  B6ETAE-BETAE 
IFC0:;{7)  ,LE.  0.)  B6ETAE-0. 

RHSC 7I=RMSC7)  ♦BBETAE  ♦  H00R1*!GKST*GAMST  ♦ 

1  OOC  /I'^CUS)  -  lEFMl 
C 

C  RHS  DUE  TO  POTENTIAL  ENERGY  OF  U6 

C 

B8ET4T  :0UfT*l-XI I 

IFCASSCiiBETATl  .GT.  9ETATI  BIETAT-DETRT 
IFCODC  in  .L£  .  U.  I  BBETAT-0. 

RHSCl  I  )*RHSm  I  ♦  BSETAT 
C 

c  RIGHT  HAND  SIDES  DUE  TO  DISSIPATIVE  ENERGY  OF  UBARl 

C 

SUMnO.  • 

SliHZ^o. 

SUM3»0. 

SUM4»0. 

DU  ?2  J=4,10 

3UM-SUHI  ♦  PDCl,J,12,n*CDC2l 
SUM2=SUM2  ♦  PDll ,J,12,2I*CD( Jl 
SUH3=5UM3  +  PCCl , J,12,3I*CD( Jl 

22  SUM4,SUH4  ♦  PD n , J , 1 2 ,4 1 »C0( Jl 
DU  23  Ki4,10 

23  RHS(K)>RHS(K)  •  (C11«SUH1«PDII,K,12,1I  •  CC12»SUtt2*PD(l ,K,12.2}  • 

1  CCZnSUHSoPCCl  .R.n.Bl  *  CC22^SUM4^P0llfR, 12,61 

C 


j 
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AIOHT  HAMO  sides  due  to  dissipative  ENEhCY  OF  UBAfl2 


sufn»o. 

suHi-o. 

DU  21  J*l.lfOS 

SUM1«SUH1  *  PD(2  I  J.12»l ••cot J) 

21  SUH2«SUM2  *  PO (2,J|12»21*CD( Jl 
DU  20  K-I.IECS 

20  <THS|K}«RhS(KI  «  CBIICE*(SUP1*P0(2»K»12,1  )  *  SUM2*P0(2  ,K ,  12,2 }  I 

KIOHT  HAND  SIDES  DUE  TO  DISSIPATIVE  ENERCV  OF  UBAR3 

SUN1«0. 

SUN2-0. 

DU  18  J«l,ll 

SUH1<SUN1  «  PDI3,J,12,1I*CD{JI 
le  SUK2«SUK2  *  P0(3,J,12,2)*CDI Jl 
DO  19  K>l,n 

19  KHS(K1>KHS<K  I  •  C ETA* < PD  I  3  ,K  ,  12 1 1 1 pSUMI ••2  *  PD( 3  ,K . 12 ,2  I •SUM2**2 1 

BRING  ALL  TERMS  FROM  LEFT  SIDE  OF  EQUATIONS  TO  RIGHT  HAND  SIDE 

DO  24  1-1,11 
24  RHSn)>-RHSin 

GENERALIZED  FORCE  DUE  TO  BREECH  FORCE 

00  too  J-1,11 
DO  100  L«1.3 

ICO  PGIl,J.12,LI>PKE(l,J,12,Ll*PT|4,J,12*Ll«PTIt.J.t2,LI«PGfl,J,12,Ll 
DO  102  J-I.IEQS 
5UN-0. 

DO  101  L*t.3 

101  SUM-SUM  *  PGI4,12,12,LI*PCI1,J,12.LI 

102  RHSIJI-RHSCJl  ♦  SUM 

GENERALIZED  FORCE  DUE  TO  ROD  PULL 
DO  103  L*1.3 

103  PGI2«1«12,L  )*PGI2,1,12,LI«PKE(1,I,12,LI-*PTU,1,12,L>*PTI6,1,12,L1 
SUM-O. 

DU  104  L*t,3 

104  SUM-SUN  *  PG(«,12,12,L)*PC(2,1,12,LI 
RHSil >-RHSIl I  ♦  SUN 

RETURN 

END 


SUBROUTINE  L 1 NE AR ( A ,X, V« VV,M , I  I 
DIMENSION  XI lOS) lYllOS) 


THIS  SUBROUTINE  CBTAINS  VALUES  BETWEEN  ADJACENT  ENTRIES  BV  LINEAR 
INTERPULATION.  LAGRANGE'S  INTERPOLATION  FORMULA  IS  USED. 


IFII  .LE.  Oli-1 

2  IFIA-Xtin3,l,l 
1  1-1*1 

IF(M-n3.3,2 

3  I-I>1 

VV*Y(  I  )-(A-XI  1*1)  I/I  XII  )-X|I*ni*Y(  1*11*1  A-Xin)/IX|l*l)>Xlin 

RETURN 

END 
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OJOl 

ooc^ 

0003 

000<> 

ooos 

0006 

0007 

0006 

0006 

0010 

0011 

0012 

OOlJ 

OOU 

0016 

0016 

0017 

ooie 

0019 

0020 

0021 

0022 

0023 

0029 

0025 

0026 

0027 

0026 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0036 

0039 

0040 

0041 

0042 

0043 

0044 

0045 

0046 

004  7 

0048 


31.0CK  data 

:  LfHmjN  /0ATA1/11ME,T  lMEH,TlHEH2,nMeH8 

C1;MHUN  /DATAi/  l»>N2tX463,C2>03.6£TlE.CKST.ELANCtCUS 

CCMH(JN  /0ATA3/e:TAT 

Ct'XMUN  /DEAtvl/g<n)tQD(tl  >.000(1!  I 

/OE^Ik2/:>l»A2,A3,ASTAR.KKTl  .AKV2 
C  JHHON  /DE4]  v3/llSUfi.  A2SUi  .A3SU9.XkKl  .A|(82,A1BAR 
<:  '  NMLiN  /0E9  Ik4./XL  12,2)  ,XH  (2,2)  .XNI  2.2  )  ,Xk(2,2  I 
C  JMHtJN  /rEftik5/EI,f2.i;3.Ffl.FF2.FF3.Dl.C2,03.Xll,ETAl.ZETAl 
/DEXU6/ai,B2.l)3.6ieAX.KI,IETA 

Ct’HHiJN  /C£eive/XtB,t(i.UTA6.x)(t,ER.lETAft,XICtEC.ZET6C 
Li'HMilN  /DEfi)  Itf/eLFl  ,CC)Ft  .KOFT  ,FUFC 
CJKHON  /NAME2/IECS.IMAS5 
CCHMDN  /NA4t3/l«A5SI5) 

tOKMDN  /NAM£*,/ltXX(5),XIVY(S),XIZZ(S),XiXY(5).XIYZ(5).XIXZ(5) 
COWHlJN  /NAKEt/C8AV,C8RCE,8EIA 
CUHMON  /NANre/CCll,CC12,CC21,Ct22 
COHHON  /kUlTA2/l6DFT .jRUFT.IGOFT 

.AT A  0(1 >.C(2), 0(3  ).Q (4), C(S).0{6),C( 7),0(8 I *019). 0(101,01 li) 

1  / 83. 675, 0.0, 0.0, 0.0, -.14 ,0.0, 0.0, 0.0, 0.0, 0.0, 0.0/ 

DATA  OD/ll»0.0/ 

DATA  Al,  A2.A3,ASTAR/0.0,19.6,-13.6,51.0/ 
data  AISUB,A2SUS,A35U6/39.5.-125.1.-7.15/ 

DATA  XKK1,XKX2/100000. ,100000./ 

LATA  A1BAR/-39.5/ 

data  XK(  1  ,1)  ,XK(1 ,2) ,XK(2 ,1 1 ,XKt 2. 2 1/4*100000./ 

data  XL(1,1),XL(1,2)  .XL (2,1 >.XL (2,2 )/53. 0,53.0, -53.0,-53.0/ 

DATA  XH(l.l)  .XCI  1,2),  XHt2,n.XH(2, 2  1/48. 0  3,-99. 97. 40.03,-99.97/ 
DATA  X0(l.ll,X6(l,2),XM2,l),XNI2,2)/-32.4,-32.4,-32.4,-32.6/ 
DATA  B1.82.B3,B18AR/39.S,-124.6.9.6.-39.S/ 

DATA  El.E2,E3/0.0,-84.6,2.0/ 

DATA  FF|  .FF2,FF3/0.0.S.2.16.65/ 
data  01,02,03/0.0,27.375,47.85/ 

DATA  XI1,ETA1.Z£TA1/0.0,23.625,-11.25/ 

data  xi,zeta/o.o,i.b/ 

DATA  XKYl ,XAY2/2*100000./ 

DATA  IE0S.IHAS$/11,S/ 

DATA  XHASStl  ),XMASS(2 ),XHASSI3) ,XMA5S(4),XNASS(5)/88.601  , 

1  4.663,20.207,10.868.36.300/ 

DATA  Xl8,Ee.ZETA8/0.0,-16.0,1.8/ 

DATA  XlR,ER,ZETAR/-3.93e,-14.0,-14.3I2/ 
data  X 1 C, EC, ZET AC/3.217,-14. 0,-12. 062/ 

DATA  80FT,C0FT,RCFT/3*0.0/ 

DATA  XIXY/5*0.0/ 

DATA  XlYZ/5-0.0/ 

DATA  XIXZ/5*0.0/ 

DATA  XIXXI 1 )  ,K|XX|2)  ,XIXX{3) ,XIXX(4) ,XIXXI5) 

1  /238342..1SC2.6,6476.7,12953.,132124./ 
data  XIYV(1),x1YY|2)  ,XIYY(3) ,XlVV(4| ,XlVVI5! 

1  /93 782.  .6  1 14.,  3 68 6. ,22  79.0,673.6/ 

DATA  XIZZd  )  ,XIZZt2),XlZZI3)*XUZ(4),XIZZ(S) 

1  /260622..15&2.6,38e6.,2279.B,132124./ 

DATA  GRAY, CORCE, PET A/386. ,0.0, 0.0/ 


0099  DATA  CCll,CCI2,CC21,CC22/4*0.0/ 

0050  DATA  TIME ,T|hEM/0. 0,0. 0005/ 

0051  data  iboft,iruft,icoft/i,i,i/ 

0052  DATA  XNN2,XS).3,02 ,03 ,80  T  AE  ,GKST  ,EL  ANG.CUS 

1  /73. 0,-10. 5, 27. 375, 65. 5,  100000. ,3*0./ 
0053  data  eiTAT/O.O/ 

0054  END 


! 
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r 


IlM; 


€TA 

o,ooos 

V 

X 

Y 

Z 

PHI 

GAH 

NU 

THET 

PSI 

TAU 

S3. 6747 

-0.0004 

-0.0000 

-0.0000 

-0.1400 

-0.0000 

0.0000 

-0.0000 

-0.0000 

0.0000 

-o.ooco 

-1.3103 

-1.4131 

0.0000 

-O.OflSt 

-0.0015 

-0.0001 

0.0001 

-0.0155 

0.0000 

-0.0000 

0.0000 

-3173.0664 

-aei4.727i 

0.0040 

-173.9220 

2.1001 

-0.1507 

-0.0376 

-30.7734 

0.0005 

-0.0000 

0.0002 

0.0010 

83.6736 

-0.0014 

0.0000 

-0.0001 

-0.1400 

-0.0000 

0.0000 

-0.0000 

0.0000 

-0.0000 

0.0000 

-3.0575 

-2.6090 

0.0000 

-0.1744 

-0.0004 

-0.0002 

0.0001 

-0.0307 

0.0000 

-o.cooo 

0.0000 

-3815.5720 

-2764.5732 

0.0031 

-176.7103 

2.3345 

-0.1157 

-0.0559 

-29.6741 

0.0004 

-0.0000 

0.0001 

0.0U15 

83.6716 

-0.0032 

0.0000 

-0.0002 

-0.1400 

-0.0000 

0.0000 

-0.0000 

0.0000 

-0.0000 

0.0000 

-5.0464 

-4.1709 

0.0000 

-0.2637 

0.0008 

-0.0002 

0.0000 

-0.0451 

0.0000 

-0.0000 

0.0000 

-4147.6633 

-2660.6604 

0.0013 

-180.2143 

2.3917 

-0.0728 

-0.0765 

-27.8573 

0.0(302 

-0.0000 

0.0001 

0.0020 

83.6685 

-0,0056 

0,0000 

-0.0004 

-0.1400 

-0.0000 

0.0000 

-0.0001 

0.0000 

-0.0000 

0.0000 

-7.1326 

-5.4830 

0.0000 

-0.3549 

0.0011 

-0.0003 

0.0000 

-0.0584 

0.0000 

-0.0000 

0.0000 

-4193.2461 

-2563.6677 

0.0043 

-186.0071 

-7.6169 

-0.1948 

0.3169 

-25.3544 

0.0005 

-0.0000 

0.0002 

0.002S 

83.6643 

-0.0066 

0.0000 

;0.0006 

-0.1400 

-0.0000 

0.0000 

-0.0001 

0.0000 

-0.0000 

0.0000 

-9.8062 

-6.7294 

0.0000 

-0.4487 

-0.0003 

-0.0003 

0.0001 

-0.0  703 

0.0000 

-0.0000 

o.oooc 

-6827.7148 

-2417.9802 

0.0021 

-190.8412 

-7.8755 

-0.1386 

0.2923 

-22.2401 

0.0003 

-0 .0000 

0.0001 

0.0030 

83.6583 

-0.0123 

0.0000 

-0.0008 

-0.1400 

-0.0000 

0.0000 

-0.0001 

0.0000 

-0.0000 

0.0000 

-14.4217 

-7.6962 

0.0000 

-0.5447 

0.0007 

-0.0003 

0.0000 

-0.0805 

0.0000 

-0.0000 

o.oooc 

-11320.4141 

-2243.2422 

-0.0004 

-194.7598 

1.6270 

0.0962 

-0.1506 

-18.5563 

-0.0000 

0.0000 

-0.6009 

0.0035 

e3.6«.94 

-0.0165 

0.0000 

-0.0011 

-0.1400 

-0.0000 

0.0000 

-0.0002 

0.0000 

-0.0000 

0.0000 

-21.4621 

-8.9698 

0.0000 

-0.6476 

-0.0001 

-0.0002 

-0.0000 

-0.0688 

0.0000 

-0.0000 

0.0000 

-16840.9570 

-2047.8862 

0.0026 

-216.8883 

-4.7696 

0.1973 

0.1109 

-14.4264 

0.0003 

-0.0000 

0.0001 

0.0040 

83.6364 

-0.0212 

0.0000 

-0.0014 

-0.1400 

-0.0000 

0.0000 

-0.0002 

0.0000 

-0.0000 

0.0000 

-31 .3036 

-9.9379 

0.0000 

-0,7610 

0.0008 

-0.0000 

-0.0001 

-0.0949 

0.0000 

-0.0000 

o.cooo 

-22525.3^55 

-1825.4407 

0.0014  - 

-236.5583 

8.4888 

0.6485 

-0.4617 

-9.6626 

0.0002 

-0.0000 

0.0001 

0.0045 
83.61 77 

-0.0264 

0.0000 

-0.0019 

-0.1400 

-0.0000 

0.0000 

-0.0003 

0.0000 

-0.0000 

o.cooo 

-43.7253 

-10.7916 

0.0000 

-0.8854 

0.0009 

0.0003 

-0.0002 

-0.0986 

0.0000 

-0.0000 

0.0000 

-27154.6914 

-1587.4905 

0.0026 

-259.7427 

1.8229 

0.7590 

-0.2020 

-5.0231 

0.0003 

-0.0000 

0.0001 

o.onso 

83.5922 

-0.0320 

0.0000 

-0.0023 

-0.1400 

-0,0000 

-0.0000 

-0.0003 

0.0000 

-0.0000 

0.0000 

-50,5205 

-11.5215 

0.0000 

-1.0208 

0.0026 

0.0008 

-0.0003 

-0,0998 

0.0000 

-0  .0000 

0.0000 

-32026.2969 

-1334.2512 

0.0017 

-281.7900 

4.9631 

1  .0443 

-0.3588 

0.0175 

0.0002 

-0.0000 

0.0001 

0.0055 

63.558B 

-0.0379 

0.0000 

-0.0029 

-0.1400 

0.0000 

-0.0000 

-0.0004 

0.0000 

-0.0000 

0.0000 

-75.7521 

-12.1219 

0.0000 

-1.1681 

0.0060 

0.0014 

-0.0005 

-0.0985 

0.0000 

-0.0000 

O.OCCQ 

-36900.1641 

-1068.9670 

-O.OOOl 

-307.1282 

8.7887 

1.3725 

-0.5404 

5.1234 

-0.0000 

0.0000 

-0.0000 
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